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The  Federal  Aviation  Administration  and  the  National 
Severe  Storms  Laboratory  are  cooperating  in  search  of 
improved  methods  for  severe  storm  prediction  and  warning 
for  aviation.  Here  NSSL  Operations  staff  reports  on  tests 
involving  transmission  of  contour-mapped  WSR-57  weather 
radar  from  NSSL  headquarters  to  a- display  unit  at  the  Okla- 
homa City  Flight  Service  Station. 

Our  study  follows  other  investigations  of  the  com- 
parative value  of  various  radar  systems  for  severe  storm 
surveillance.  Improved  signal  processing  and  communication 
techniques  and  equipment  now  permit  rapid  dissemination  of 
Information  concerning  storm  location,  intensity,  and  move- 
ment and  offer  a new  dimension  in  weather  display  for 
general  aviation. 
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Computer  Applications  and  Techniques  for  Storm  Tracking  and  Warning 

W.  David  Zlttel 


1.  INTRODUCTION 

This  report  extends  tests  of  the  remote  radar  display  described  in 
Volume  1 and  examines  the  feasibility  of  the  display  as  a graphics  termi- 
nal. A storm  tracking  program  has  been  combined  with  an  echo  contouring 
scheme  to  produce  graphic  warning  areas  based  on  size  and  motion  of  storm 
echo  areas. 

A detailed  description  is  provided  of  the  mathematical  techniques 
employed  and  a quasi-real  time  software  program  is  outlined.  In  addi- 
tion, three  case  studies  utilizing  the  above  logic  are  presented.  i 

Finally,  a summary  of  results  and  suggestions  for  Improvements  and  future  j 

work  are  discussed.  /i 

I 

2.  BACKGROUND 

Information  regarding  storm  motion  and  growth  tendencies  are  now  pre- 
sented to  the  FSS  pilot-brief er  in  two  forms.  First,  a numerical  coding 
in  the  hourly  radar  report  transmitted  by  teletype  (RAREP)  Indicates  the 
past  tendency  of  storm  pattern  growth  or  decay.  Motion,  in  polar  coordi- 
nates, of  both  the  pattern  and  Individual  cells  are  Included. 

Secondly,  a plain  language  summary  provides  a "layman's"  geometric 
description  of  the  storms  with  geographical  references  to  outline  the 
present  and  projected  coverage.  Severe  Storm  Warnings  and  special  advi- 
sories to  airmen  (SIGMETS  and  AIRMETS)  carry  Information  on  hazardous 
flight  conditions.  At  a few  locations,  these  messages  are  augmented  by 
a facsimile  machine  replica  of  the  Plan-Position  Indicator  (PPI)  with 
appropriate  annotations  (Bigler,  1969). 

The  reliability  of  both  types  of  advisories  varies  directly  with  the 
spatial  and  temporal  variance  of  radar  echo  patterns.  Information  con- 
tained In  the  RAREP  Is  usually  a sterile  summary  of  the  radar  scope  dis- 
play, condensing  details  observed  and  coded  during  a specific  15  minute 
period.  Plain  language  summaries  and  advisories  may  Include  Information 
on  the  position  and  movement  of  fronts  and  squall  lines,  and  observations 
of  recent  severe  weather  events. 

During  periods  when  echo  coverage  and/or  Intensity  change  rapidly 
special  observations  supplement  hourly  reports.  During  periods  of  severe 
weather,  the  National  Weather  Service  radar  scopes  are  monitored  con- 
stantly, but  because  the  flow  of  Information  is  restricted  by  communication 

1 


j 


facilities  and  the  heavy  work  load  required  to  meet  various  local,  state 
and  national  needs,  messages  to  the  FSS  are  periodic. 

The  Volume  I tests  have  shown  that  if  calibrated  contoured  data  are 
available  routinely  at  the  FSS,  pilot  briefers  can  interpolate  between 
National  Weather  Service  advisories  and  maintain  a "user's  watch"  of 
storm  locations  and  Intensities.  However,  neither  time  nor  expertise 
is  available  at  the  FSS  to  relate  echo  patterns  to  s3moptlc  scale  dis- 
turbances (wind,  pressure,  and  moisture  fields)  and  severe  weather 
reports.  Even  with  such  data,  it  is  difficult  for  meteorologists  to 
predict  changes  in  gross  features  of  precipitation  areas. 

Fortunately,  large  severe  storms  tend  to  be  steady-state  and  lend 
themselves  to  tracking  and  extrapolation.  The  principal  objective  of 
this  study  is  to  apply  semi-automated  computerized  logic  to  identify, 
track,  and  extrapolate  those  storms  of  sufficient  intensity  and  size 
to  produce  hazardous  weather  conditions  and  to  map  out  a warning  area 
for  the  extrapolated  storm  positions. 


3.  RATIONALE  FOR  SELECTING  AND  TRACKING  ECHO  CENTROID 
(OR  WHY  LEAVE  A PERSON  IN  THE  PICTURE) 

Several  years  experience  in  field  operations  at  NSSL  have  stressed 
the  value  of  retaining  the  meteorologist  for  real  time  decision  making. 
It  seems  difficult,  if  not  impossible,  to  anticipate  all  the  compli- 
cated elements  occasionally  present  in  real  weather  situations,  in  a 
computer  program. 

Several  objective  techniques  have  been  suggested  by  Kessler  and 
and  Russo  (1963),  Wilson  (1966),  and  Blackmer  and  Duda  (1972),  which 
rely  on  spatially  correlating  PPI  information  to  derive  echo  motion  and 
speed.  As  shown  in  Figure  1,  significant  storms  on  the  radar  scope  may 
have  quite  different  motions.  Also  under  certain  conditions  severe 
storms  split  with  one  portion  often  moving  to  the  left  of  the  average 
ambient  wind  direction  and  faster  than  the  wind  speed,  while  the  other 
portion  moves  to  the  right  and  slower  than  the  ambient  wind  (Newton  and 
Fankhauser,  1964).  Under  these  conditions  a single  speed  and  direction 
of  motion  for  the  whole  scope  would  be  misleading. 

Even  if  individual  echoes  are  first  isolated  (Wllk,  1966),  there  are 
several  drawbacks  to  using  this  approach.  In  a matrix  analysis,  a mini- 
mum of  two  PPI's  must  be  stored  at  the  same  time  requiring  a large  com- 
puter core.  Generally,  the  two  fields  should  have  uniform  grid  density 
which  an  R,  9 system  doesn't  have.  In  an  R,  9 system,  data  must  be 
scan-converted  to  rectilinear  coordinates  before  correlating  the  data. 
Both  scan  conversion  and  correlation  techniques  are  time  consuming. 

Care  also  must  be  taken  when  scan  converting  to  assure  that  spatial 
averaging  has  not  changed  the  distribution  of  intensity  integers. 
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By  contrast,  the  echo  centroid  extrapolation  technique  requires 
only  a small  amount  of  core  and  is  very  fast.  One  can  operate  the 
program  with  radar  data  extracted  manually  from  the  PPI  scope  display. 
One  may  also,  in  a relatively  short  period  of  time,  use  limited  auto- 
mation to  scan  a PPI  for  centroid  information,  and  display  it  regularly 
without  full  time  monitoring.  (Such  a technique  is  presented  in  sec- 
tion 5.)  Tests  during  the  NSSL  Spring  Program  (Wllk  and  Gray,  1970) 
indicate  an  operator  can  easily  filter  extraneous  or  unwanted  data. 

Some  storms  may  be  moving  beyond  the  radar  scope's  range  while  others 
may  be  part  of  a broad  band  of  non-severe  stratiform  rain  whose  overall 
movement  is  slower  and  more  persistent  (fig.  2). 

An  operator  may  recognize  splitting  or  merging  storms  which  need 
to  be  treated  as  new  echoes.  (Computer  programs  to  date  have  not 
proved  reliable  in  echo  matching  and  we  make  no  attempt  to  do  this 
here . ) 

One  final  reason  for  leaving  an  operator  in  the  picture  is  to 
Insure  detection  of  system  failures  and  to  recognize  spurious,  non- 
meteorologlcal  results  when  computerized  objective  analysis  software 
systems  are  in  operation.  The  following  sections  are  devoted  to  the 
operation  of  a man-machine  mix  using  examples  of  real  data  sets. 


Figure  1.  WSR-57  radar  PPI,  100  n ml. 
range,  20  n ml.  range  marks, 

1454  CST,  April  3,  1964.  Indi- 
vidual echoes  are  numbered  one 
to  five. 


Figure  2.  WSR-57  radar  PPI,  200  km 
range,  40  km  range  marks,  1132  CST, 
September  24,  1974.  Light  to 
moderate  stratiform  showers 
Indicated  over  most  of  the  radar 
scope. 
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4.  MATHEMATICAL  FORMULATION 
4.1  Introduction 


Three  sources  of  information  are  used  to  construct  a graphic  pres- 
entation of  a warning  area.  In  order  of  calculation,  they  are  (a)  echo 
centroid  and  shape,  (b)  echo  motion,  and  (c)  a measure  of  the  variance 
of  the  echo  motion.  The  method  used  to  calculate  echo  centroid  and 
shape  basically  requires  fitting  an  arc  length  function  to  the  echo's 
perimeter  and  was  suggested  by  Blackmer  and  Duda  (1972),  later  developed 
by  Ostlund  (1974).  Calculation  of  echo  motion  and  variance  using  cen- 
troid positions  was  developed  by  Barclay  and  Wilk  (1970)  and  run  opera- 
tionally during  the  1970  Spring  Data  Collection  Program. 

4. 2 Echo  Shape  and  Centroid  Calculation 

Use  of  the  arc-length  function  to  describe  echo  shape  requires  that 
one  first  determine  echo  perimeter.  In  the  computer  logic  developed  for 
this  report,  data  are  entered  into  core  and  all  bins  with  intensity  less 
than  a specified  level  are  first  set  to  zero.  Then,  beginning  with  zero 
degrees  azimuth,  the  PPI  is  searched  until  an  echo  is  found.  Then  the 
program  isolates  the  echo,  moving  around  the  perimeter  in  a counterclock- 
wise direction  until  it  comes  upon  the  starting  point.  This  logic  dif- 
fers from  Ostlund's  in  at  least  two  respects.  First,  the  echo's  perim- 
eter is  defined  in  an  R,0  coordinate  system;  and  second,  the  program 
minimizes  echo  area.  The  following  two  examples  in  B scan  format  illus- 
trate these  points. 


In  Figure  3 the  arrows  indicate 
the  path  followed  in  the  boundary 
search.  S is  the  starting  gatt , X's 
represent  echo,  dots — no  echo. 

Echo  1 is  joined  to  echo  2 by  a 
single  gate  along  a common  radial. 
The  program  ignores  that  gate  since 
it  would  have  to  be  used  twice  in 
order  to  close  the  boundary  and  com- 
bine echo  1 and  2.  Likewise,  between 
echo  2 and  echo  3 there  is  a coinnon 
comei . But  because  the  comer  gates 
would  have  to  be  used  twice,  the 
echoes  are  separated. 


No  gate  is  used  twice  except  the 
starting  one  and  no  gate  is  accepted 
unless  there  is  an  acceptable  gate 
beyond  it.  The  subroutine  BNDRY 
(called  from  CONTUR)  contains  this 
logic . 


RANGE  (KM) 

21  22  23  2^  25  26  27  28  29  30 


Figure  3.  Representative  echo 
samples  show  the  path  and  order 
of  points  to  describe  echo 
perimeter. 
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The  obvious  result  of  the  above,  is  that  cores  tend  to  be  discrete 
with  a more  regular  shape. 

After  a closed  boundary  Is  found,  the  area  of  echo  Is  calculated  by 
summing  up  all  bins  within  and  Including  the  perimeter.  The  area  of  a 
bin  is  given  by  (SirRi^  _ 0irR22)/36O®,  which  can  be  factored  into 
0tv(R2  ” R2)(Ri  + R2)/360®.  Since  the  difference  between  R^  and  R2  is  the 
gatelength  of  the  radar,  G,  and  0 Is  the  angular  difference  between 
radlals,  the  area  of  an  Individual  bin  Is 


07tG(2R3^  - G) 
360® 


(1) 


If  an  echo's  area  Is  less  than  a specified  threshold,  it  Is  Ignored 
and  the  program  searches  for  a new  echo.  If  an  echo  exceeds  the  speci- 
fied area,  Fourier  analysis  of  Its  shape  Is  performed.  If  the  area  Is 
more  than  five  times  the  specified  area,  the  lowest  intensity  is  purged 
and  the  remaining  core  treated  as  a new  echo.  This  process  is  iterated 
until  the  echo  is  less  than  the  specified  area. 

Once  an  echo  meets  the  size  criterion,  the  program  enters  subroutine 
OSTLND.  Here  the  paired  azimuth  and  reuige  perimeter  data  are  converted 
to  Cartesian  points.  Fourier  analysis  of  X(s)  and  Y(s)  Is  performed 
where  s is  the  arc  length  function.  Mathematically  these  functions 
are: 


X(s) 


00 

■ y a cos 
n*0 


/2n7TS\  1 . . /2n7TSv 

(— > + ^n 


(2) 


V 


.. y \ c /2n7TS\  ( . . ^zn*n^S\ 

Y(s)  = I c cos  — ) + d^  sin  (— r — ) 
n 1.  n L 

n=0 


,2mrs. 


(3) 


The  coefficients  a , b , c , and  d may  be  expressed  as 
n’  n’  n’  n 


n L 


Y / X(s)  cos  ds 


(4) 


Y / X(s)  sin  (-?^)  ds 


n L 


(5) 


: =■  |-  / Y(s)  cos 

n L ^ L 


(6) 


5 


(7) 


I 

I 


1 

(i 


I 

t 

r 


dn  - f A (8)  sin  (^)  ds 

where  n,  an  Integer,  Is  the  nth  harmonic. 

Each  coefficient  may  be  rewritten: 

, M ®1+1  , 

®n  “ L ^ ^ (-^^)  ds 

“ ^ i-1  ^ 


(8) 


where  M is  the  number  of  discrete  points  in  the  echo's  boundary.  Also, 
since  X(s)  may  be  considered  as  consisting  of  a series  of  discrete  line 
segments,  one  may  set  X(sj^)  “ Pi  + qis^.  Setting  this  expression  into 
Sq.  (8)  and  Integrating  yields: 


7 “ 

r ^ 

^ 1-1 


Pi^  2n7TSi  q^L 

■z sin  ( ;; )+ — 

^ (2mr)2 


cos  (■ 


2nirsi 


2n7rs. 


’1+1 


(9) 


Eq.  (9)  and  the  corresponding  equation  for  each  of  the  other  coefficients 
are  calculated  in  the  computer.  The  0th  harmonic  yields  the  mean  of  each 
series  and  thus  the  centroid  of  the  echo's  shape.  Eight  harmonics  in 
addition  to  the  mean  are  calculated. 

Because  the  values  are  derived  initially  from  a polar  scan,  the- 
density  of  points  about  the  perimeter  Is  not  constant,  biasing  the  cen- 
troid location  towards  the  radar.  In  a h3rpothetical  case,  using  a circle 
of  10  Von  radius,  the  resulting  centroid  error  varies  as  a function  of 
range  (fig.  4)  and  Is  a maximum  at  10  km.  However,  for  an  echo  with 
stable  motion  there  is  little  error  because  all  centroids  have  the  same 
bias. 


Figure  4.  Graphic  presentation 
of  error  in  determining  the 
echo  centroid  due  to  non- 
uniform  perimeter  density  for 
a hypothetical  circular  echo 
with  a 10  km  radius. 
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4.3  Echo  Motion  Calculation 


The  calculation  of  echo  motion  uses  linear  least  squares  (LSS)  equa- 
tions fitted  through  an  echo's  past  centroid  positions  expressed  para- 
metrically as  a function  of  time  as 


X(t)  * A t + B 

X X 

(10) 

Y(t)  = A t + B 

y y 

(11) 

1 where  A^  and  A^  are  found  by  solving 

. NZXt  - EXEt 

(12) 

j 

^x  * 2 2 

NEt  - (Et)^ 

I 

^ NEYt  - EYEt 

(13) 

i 

^ NEt^,-  (Et)^ 

1 The  ordinate  axis 

Intercepts,  usually  found  by  solving 

■ 

‘ 

B = X - A t 
X X 

(14) 

\ and 

B « Y - A "t  . 

y y ’ 

(15) 

: are  here  given  as 

• 

1 : 

j 

B - X.  - A t. 

X f,  X £. 

(16) 

/ and 

B - Y.  - A t 

y y 

(17) 

where  t£,  Xjj  and  Yjj,  are  the  echo's  last  position  In  time  and  space.  This 
condition  forces  the  LSS  equations  through  the  last  point. 

From  Eqs.  (10)  and  (11)  echo  speed  Is  simply 

SPD-  (A^^  + (18) 

and  the  direction  of  motion  Is 

DIR  - TAN~^(A  /A  ) + 180  . (19) 

X y 
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A measure  of  an  echo’s  predictability  in  time  and  distance  is 
determined  by  comparing  the  time  of  the  latest  centroid  position  to  the 
predicted  tiro  of  closest  passage  to  that  centroid  position.  The  dif- 
ference, At,  is  defined  as  the  error  in  time,  e^.  The  distance  between 
the  predicted  point  of  closest  passage  and  the  actual  centroid  location 
is  defined  as 

^d  “ ^^y’'  ®y  ■ (20) 

where  t is  the  time  of  closest  passage  and  X and  Y are  the  Cartesian 
coordinates  of  the  latest  centroid  position.  We  can  solve  for  the 
unknown  time,  t,  by  first  squaring  terms  in  Eq.  (20)  and  then  differen- 
tiating them  with  respect  to  t yielding: 

— 3J—  - 2A^(A^t  + - X)  + 2Ay(Aj,t  + - Y)  . C21) 


Setting  the  expression  on  the  right  equal  to  zero  and  solving  for  t yields 

(22) 


t = 


A(X-B)+A(Y-B) 
X X y y' 


2 2 
A + A ^ 
X y 


Therefore  is  t - t£  where  is  the  time  of  the  latest  centroid. 
Given  t,  e^j  can  be  calculated  directly  from  Eq.  (20). 

Finally,  and  are  normalized  to  one  hour  and  root  mean  square 


errors  (RMSE)  computed  from 


2“ll/2 


and 


RMSEj.  = 


(23) 


(24) 


i 

i 


for  n ^ 3,  where  n is  the  nimber  of  points  in  a discrete  track. 


4.4  Warning  Area  Calculations 

We  can  now  combine  the  results  of  sections  4.2  and  4.3  to  determine 
a warning  area.  Given  a beginning  time,  t|,,  and  ending  time,  tg,  we 
first  solve  Eqs.  (10)  and  (11)  for  the  starting  and  ending  points  P],  and 
Pg,  of  the  warning  area  specified  in  the  time  domain,  which  lie  on  the 
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echo  path  (fig.  5).  A measure  of  a storm's  predictability  in  time  is 
included  in  determining  P|j  and  as  follows: 

For  Pjj 

^ - ^t>  + ®x  (25) 

Y = - O^)  + By  (26) 

and  for  P 

X = A^(t^  + a^)  + (27) 

Y = Ay(tg  + o^)  + By  . (28) 


Next,  two  line  segments  are  found  which  are  parallel  to  and  the 
same  length  as  the  major  axis  of  the  echo  and  which  pass  through  Pjj  and 
Pg,  respectively.  The  length  and  orientation  of  the  line  segments  are 
determined  by  finding  an  ellipse  which  approximates  the  echo  at  hand. 


A parametric  form  of  an  ellipse  is  given  by  the  zeroth  and  first 
harmonics  of  the  echo.  For  convenience  the  ellipse  is  translated  to  the 
origin  eliminating  the  zeroth  harmonic  from  further  calculations.  From 
Eqs.  (2)  and  (3)  Cartesian  coordinates  expressed  as  a function  of  the 
arc  length,  s,  for  any  point  on  the  ellipse  are  given  by 


Figure  5.  Illustration  of  use  of  an  approximating  ellipse,  linear 

least  squares  predicted  echo  trajectory  and  the  RMSE  values 
to  calculate  graphic  warning  area  for  oi  > 10  degrees. 
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X(s)  - a^cos(^)  + b^8ln(-^) 


Y(s)  - c^cos(^)  + d^8ln(^) 


Al8o,  an  elllpae  haa  the  aecond  degree  form 

AX^  + BXY  + CY^  - K (31) 

where  A,  B,  C and  K are  conatants.  By  combining  Eqa.  (29)  and  (30)  with 
Eq.  (31)  we  can  find  three  equations  with  which  to  solve  for  A,  B and  C 
from  which  the  orientation  of  the  ellipse  is  found.  A value  for  K is 
specified  below.  Setting  Eqs.  (29)  and  (30)  into  (31)  yields: 

K - A(a^cos(^)  + b^sln(^))^  + c(c^cos(^)  + d^sin(-?p))^ 

(32) 

+ B(a^cos(-^)+  b^sln(^))(c^cos(^)  + d^sin(^))  • 

Expanding  and  combining  like  terms  gives  us: 

K - (Aa^  + Bac  + Cc^)  cos^(-^)  + (Ab^  + Bbd  + Cd^)sin^(-^) 

(33) 

+ (2Aab  + B(ad  + be)  + 2Ccd)sin(^)cos(^)  • 

L L 

When  2Trs/L  - 0,  sin(0*)  » 0 and  cos(0°)  = 1;  when  2tts/L  = 7r/2, 
co8(tt/2)  ” 0 and  sin(Tr/2)  * 1. 

Under  these  two  conditions  Eq.  (33)  reduces  to  the  following  two 
identities 

2 2 

K = Aa  + Bac  + Cc  (34) 

2 2 

K = Ab  + Bbd  + Cd  • (35) 

Since  they  contain  only  constants  they  are  valid  for  all  s and  K can  be 
set  into  Eq.  (33)  as: 


K - K cos^(^)  + K sln^(^) 
u L 


+ (2Aab  + B(ad  + be)  + 2cd)  sin  (^)  cos  (-^). 

L L 
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Prom  trigonometry  sin^  + cos^»l  and  similarly  K sln^u)  + K cos^O)  = K. 
Hence,  for  Eq.  (36)  to  be  valid  when  sin(2lTs/L)  and  cos(2Trs/L)  both  ^ 0 

2Aab  + B(ad  + be  + 2Ccd)  = 0 (37] 


must  be  true.  Since  a,  b,  c,  and  d are  simply  the  coefficients  of  the  first 
harmonic,  only  A,  B,  C and  K are  unknown. 


From  Eq.  (29)  a maximized  value  of  X(s)  is 


+ b^) 


2x1/2 


and  from 


Eq.  (30)  a maximized  value  of  Y(s)  is  (c^  + d^)^/^.  Therefore  a maximum 
value  for  (X(s)2  + Y(s)2)^'2  is  a2  + + d2)l/2  ^hlch  is  set  equal  to 

K.  This  gives  us  a fairly  accurate  measure  of  the  ellipse's  semimajor 
axis.  K,  as  computed  above,  will  always  be  slightly  greater  than  the 
true  length  of  the  semimajor  axis.  However,  this  is  quite  satisfactory 
since  the  length  of  the  echo's  axis  is  the  sum  of  several  harmonics  and 
not  just  of  the  first  harmonic  alone. 


Using  determinants.  A,  B,  and  C may  be  found  by  solving  Eqs.  (34), 
(35)  and  (37)  simultaneously.  Specifically 


0 j(ad  + be) 


|ab  -jCad  + be) 


Expanding  the  determinants  and  combining  like  terms  yields: 

K(bcd2  + i(a^  - d^)(ad  + be)  + adc^) 

A — — - _ ■ ■ 

2 2 2 2 7 2 2 2 2 

2(a  bed^  - ab^c^d)  + y(b  c - ad) (be  + ad) 


B and  C are  computed  likewise. 
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Once  the  ellipse's  coefficients  have  been  found.  Its  orientation  can 
be  determined  using  the  relationship  (Morris  and  Brown,  1937) 


tan (2(d)  * (40) 

or 


There  is  an  ambiguity  as  to  which  axis  of  the  ellipse  from  which  co  is 
measured.  This  ambiguity  may  be  resolved  by  considering  the  sign  and 
relative  size  of  A,  B and  C.  However,  if  the  ellipse  is  first  rotated 
through  angle  (d  eliminating  B,  its  equation  becomes 

a'x2  + c'y2  = K • (42) 

Then  if  A*  is  less  than  c',  (D  is  measured  with  respect  to  the  major  axis; 
if  C is  greater  than  A , (d  is  measured  with  respect  to  the  minor  axis. 

Now  if  the  slope  of  a line,  the  length  between  two  points  on  that 
line,  and  the  coordinate  of  one  of  the  points  are  all  known,  we  may 
solve  for  the  coordinates  of  the  unknown  point  by  combining  the  equation 
of  a straight  line  with  the  equation  for  the  distance  between  two  points. 
Further ^assume  that  the  given  line  intersects  the  echo  path  at  (x’,Y*) 
where  X and  Y are  known  and  the  unknown  coordinates  are  X and  Y.  Then 

Y - Y'  = M(X  - X')  (43) 

where  M = tan  (d  is  the  slope  of  the  line  and 

f?  - Cl  - + (X  - X')^  (44) 


where  R is  the  distance  between  (X,Y)  and  (X’,Y').  We  square  Eq.  (43)  and 
set  it  into  Eq.  (44)  yielding 


= M^(X  - X')^+  (X  - X')^ 
Factoring  and  transferring  terms  yields 


K 2 

(X  - X')^ 

M^+1 


(A5) 


(46) 


Finally,  taking  the  square  root  of  each  side  and  solving  for  X gives  us 

(47) 


X = X’  ± (R^/(M^  + 1))^^^ 
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and  Y for  each  X is  given  by 


Y - Y'  + M(X  - X')  . (48) 

Thus,  we  have  found  two  points~one  above  and  one  below  the  echo  path 
which  define  the  initial  warning  boundary.  There  are  also  two  points 
which  define  the  final  warning  boundary.  These  four  points  define  the 
warning  area.  R in  the  above  equations  is  equivalent  to  K in  the  general 
equation  of  the  ellipse.  However,  in  solving  for  the  above  points,  R is 
modified  to  Include  a^j.  is  scaled  linearly  such  that  the  total  length 
for  R is  given  by 


*^  * *^  "^  ^d  Is  ‘ si  S (^9) 

and  R . K + |t^  - t^|  at  . (50) 


tjj  is  the  time  of  the  last  echo  observation. 


From  Figure  6 it  will  be  recognized  that  the  warning  area  is  a 
modified  parallelogram.  However,  as  the  echo's  major  axis  becomes  more 
closely  aligned  with  the  echo  motion,  the  parallelogram  closes.  There- 
fore, whenever  the  echo's  motion  and  the  echo's  major  axis  are  within 
10  degrees  of  each  other,  the  minor  axis  of  the  best  fit  ellipse  is  used 
for  K and  the  slope  of  the  lines  passing  through  Pw  and  P_,  respectively, 
are  given  as  -1/XM  where  XM  is  the  slope  of  the  echo  motion  line  passing 
through  both  Pj^  and  P^  (fig.  6). 


CASE  FOR  uxlO* 

sample  echo 

ECHO 
MOTION 


Figure  6.  Illustration  of  use  of  an  approximating  ellipse,  linear  least 
squares  predicted  echo  trajectory  and  the  RMSE  values  to  cal- 
culate graphic  warning  area  for  OJ  S 10  degrees. 
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5.  REAL  TIME  SYSTEM  AND  PROGRAM  OPERATION 


f' 


I. 
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Since  the  remote  radar  display  system  has  the  bullt-ln  capability 
to  be  Interfaced  to  a computer,  we  adopted  the  programming  philosophy 
of  duplicating,  as  nearly  as  possible,  a real-time  operation.  In  this 
section  we  shall  first  describe  a model  system  and  Its  components; 
second,  describe  the  decision  making  and  choices  within  the  software 
available  to  an  operator;  and  third,  offer  some  guidelines  for  using  the 
program  logic. 


5.1  Hypothetical  Systems  Configuration 

In  addition  to  the  electronic  components  already  described  In  Vol- 
ume I,  the  system  requires  a central  processor  with  a 50K  decimal  word 
memory  core.  In  order  to  operate  In  a pseudo-real-time  manner,  memory 
cycle  time  should  be  about  one  ysec.  (The  Systems  Engineering  Laboratory's 
model  8600,  on  which  the  software  was  developed,  has  a memory  cycle  time 
of  600  nanosec.) 

Secondly,  some  sort  of  mass  storage  unit  Is  needed.  When  not  being 
used,  the  prediction  and  display  logic  resides  there.  Otherwise,  the 
Fortran  program  would  have  to  be  entered  each  time  the  system  Is  used. 
Also  stored  on  disk  are  three  data  files:  a)  coordinates  for  graphically 
displaying  the  State  of  Oklahoma,  b)  coordinates  for  graphically  dis- 
playing the  Victor  Airways,  and  c)  a list  of  Oklahoma  airports.  (The 
use  of  these  files  Is  explained  below.)  Last,  an  I/O  device  such  as  a 
teletype  or  alphameric  CRT  with  keyboard  entry  is  needed.  The  operator 
must  manually  Insert  comnands  and  echo  Information  Into  the  software  and. 
In  turn,  receives  back  numerical  values  of  echo  speed  and  direction  pf 
motion  and  a measure  of  the  predictability  of  echo  motion.  Information  • 
flow  is  shown  in  the  systems  flow  chart  (fig.  7). 

5.2  Hypothetical  Software  Logic 

First,  let  us  assume  that  the  program  already  resides  on  disk  and 
has  been  given  the  name  ECHOPRED.  The  operator  then  merely  enters 
ECHOPRED  to  bring  the  program  to  an  operational  status.  The  operator's 
first  decision  Is  whether  or  not  to  Initialize  the  program.  (Figure  8 
illustrates  the  command  structure  which  Is  presented  in  this  section.) 
This  depends  upon  whether  or  not  the  operator  Is  working  a new  storm 
day.  For  a new  day  or  a long  break  in  operation,  the  operator's  response 
will  be  'YES',  otherwise  we  presume  he  Is  still  working  the  same  storms 
and  the  response  Is  'NO',  to  the  question,  'INITIALIZE'.  When  the  answer 
Is  'NO',  echo  Information  is  retrieved  from  disk.  The  program  should  be 
left  operational  during  storm  conditions;  only  If  a power  failure  dis- 
rupts operation  should  the  operator  not  Initialize  the  program. 

Next,  the  program  will  ask  for  'COMMAND'.  Assuming  this  Is  the 
start  of  operation  or  recovery  after  a failure,  a systems  check  should 
be  made.  After  the  operator  responds  with  'RQC',  Radar  Quality  £ontrol. 
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Figure  7.  Schematic  of  systems  flow  chart. 


the  computer  will  type  first  'PPI  CHECK'  with  appropriate  response  being 
'YES'  or  'NO',  and  then  'TEST  PATTERN'.  Again  the  operator  responds 
'YES'  or  'NO'. 

In  'PPI  CHECK'  the  program  checks  the  housekeeping  Information  (date, 
time,  azimuth)  In  detail  and  also  counts  the  number  of  bins  of  each  Inten- 
sity In  the  PPI  and  presents  this  Information  to  the  operator.  A system- 
atic decrease  In  the  frequency  of  higher  Intensities  should  occur  when 
only  ground  clutter  returns  are  present.  A low  count  at  especially  the 
first,  second  or  fourth  Intensity  levels  may  Indicate  hardware  failure. 

A few  random  housekeeping  errors  will  occur  due  to  telephone  line  noise 
and  are  not  serious. 

The  'TEST  PATTERN'  Is  a computer  generated  field  of  seven  concentric 
rings  20  degrees  In  width  corresponding  to  each  Intensity  switch  sur- 
rounded by  seven  radial  stripes  of  10  degrees  width  starting  with  0 deg- 
rees AZM,  and  repeated  through  360  degrees.  The  purpose  Is  to  check  the 
fidelity  of  the  receiver  memory.  For  least  confusion,  we  recommend  set- 
ting the  Intensity  switches  to  the  following  gray  shade  pattern: 

1 2 3 1 2 3 2. 

When  the  computer  has  finished  with  one  or  both  of  the  above  tests. 
It  will  again  type  'COMMAND'.  At  this  point  the  operator  may  wish  to 
enter  a value  for  the  Time  Weight  £onstant,  TWC.  After  TWC  Is  entered, 
the  computer  will  type  'HHMM'  for  which  the  operator  enters  a number  such 
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^ ^ T1 lustration  showing  various 

* a (Here  and  on  adjoining  pag  ) 

Figures.  i„  ECBOPRED. 
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as  2400.  If  no  entry  Is  made  for  TWO  the  program  uses  a 30  minute  default 
value.  The  TUG  exponentially  weights  the  Influence  that  past  centroids 
have  when  predicting  echo  motion,  giving  greatest  weight  to  the  most  recent 
point  (see  section  5.3.1). 

Another  parameter  the  operator  may  wish  to  change  Is  the  £round  flutter 
Mstance,  GCD.  The  default  value  Is  normally  set  at  20  km  for  the  NSSL 
radar,  to  omit  all  of  the  ground  targets  from  the  analysis.  When  the 
ground  clutter  Is  extended  by  abnormal  propagation,  spurious  echoes  may 
be  processed.  By  setting  the  size  and  Intensity  criteria  high  enough, 
these  echoes  will  be  Ignored.  However,  time  will  be  lost  determining 
this  fact.  As  echoes  move  Into  the  ground  clutter,  spurious  echoes  will 
complicate  the  shape,  but  not  seriously  affect  total  echo  area  and  cen- 
troid position.  Here  a simple  and  expedient  method  Is  to  tilt  the  radar 
antenna  at  two  degrees  which  will  effectively  remove  the  ground  targets 
from  the  scope  while  leaving  the  echo  pattern  mostly  unchanged.  When 
anomalous  echoes  are  extensive,  some  program  speed-up  can  be  realized  by 
setting  the  GCD  value  artificially  large,  say  100  km.  'However,  the  risk 
here  Is  that  the  operator  will  fall  to  reset  the  value  as  echoes  approach 
that  range. 


The  next  command  by  the  operator  Instructs  the  computer  to  accept  the 
remote  radar  data  and  locate  echoes.  Before  doing  this  the  computer  will 
reply  'AREA/INTENSITY' . The  operator  must  then  respond  with  two  values, 
for  example,  100  km^  and  the  4th  code  switch.  This  means  that  only  echoes 
whose  areas  are  greater  than  100  km^  and  whose  Intensities  are  greater  than 
or  equal  to  the  dBZ  value  corresponding  to  the  fourth  Intensity  switch  are 
contoured.  For  the  data  used  In  this  report  the  dBZ  value  Is  about  40 — a 
rainfall  rate  of  12  mm  hr“^  (0.5  In  hr“l).  The  program  also  checks  for 
echoes  whose  area  Is  five  times  that  given.  Whenever  this  occurs  the  lowest 
Intensity  In  the  echo  Is  purged  and  the  next  Intensity  level  checked  to  see 
If  It  meets  100  km^  criterion.  If  It  does,  the  Information  for  the  higher 
Intensity  core  Is  saved  as  well  as  the  lower  Intensity. 


At  the  completion  of  the  PPI,  the  computer  writes  out  each  echo's 
area  and  centroid  location  (azimuth  and  range) ; the  echoes  being  sorted 
by  Intensity.  The  program  id.ll  continue  to  accept  new  PPI  Information 
unless  Interrupted  by  a 'STOP'  command  typed  In  by  the  operator.  The 
program  does  no  matching  of  echoes  between  two  PPI's  and  no  Information 
from  a previous  PPI  Is  saved  In  the  computer  except  that  manually  entered 
by  the  operator.  Since  two  centroid  positions  are  required  to  make  a 
prediction,  the  operator  should  wait  two  scans  before  entering  time,  azi- 
muth and  range  data. 


When  the  program  has  received  a 'STOP'  command.  It  will  ask  for  a 
new  command.  The  operator  will  then  want  to  enter  the  echo  Information. 
This  Is  done  by  first  typing  'ENT'.  The  computer  will  then  ask  for  an 
echo  number,  time  of  observation  and  echo  azimuth  and  range;  It  will  type 
the  following: 

'N/  HHMM/  AZM/  RNG/  DIRECTION/  SPEED/  STDDS/  STDTM' 
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The  echo  may  have  any  number  assigned  between  one  and  99.  The  time  Is 
entered  as  a four  digit  number.  The  first  two  digits  are  for  the  hour 
(a  24-hour  day  is  used);  the  second  two  digits  are  for  the  nearest  whole 
minute.  Azimuth  is  entered  to  the  nearest  whole  degree;  range  to  the 
neeurest  kilometer.  After  the  first  two  entries  and  each  subsequent  entry 
for  the  same  storm,  immediately  following,  the  program  will  respond  with 
a speed,  direction  of  motion  and  the  standard  deviation  in  distance  and 
time  of  the  echo's  motion  (cf  section  4).  The  program  checks  the  manu- 
ally entered  data  for  simple  entry  errors.  Also,  if  the  last  time 
entered  is  the  same  as  the  time  of  the  last  PPI  scanned  by  the  program, 
the  program  will  match  the  manually  entered  centroid  position  to  the  com- 
puter derived  centroid  position.  The  echo  selected  is  the  one  whose  cen- 
troid distance  is  a minimum  from  that  manually  entered.  If  no  match  is 

found,  the  observation  is  deleted  and  an  error  message  generated.  If 
the  computed  echo  speed  is  greater  than  120  km/hr  and  error  message  is 
also  typed  out  but  the  observation  is  not  purged.  (Echoes  moving  at 

that  speed  are  rare.)  When  there  are  no  more  observations  to  enter,  the 

operator  enters  a zero  for  the  echo  number;  the  program  will  respond  by 
asking  for  a new  command.  One  last  operation  which  can  be  performed  is 
to  delete  an  observation  by  entering  a minus  sign  in  front  of  the  echo 
number  followed  by  time  and  centroid  positions.  Up  to  ten  different 
echoes  can  be  stored  at  one  time. 

There  are  four  different  commands  the  operator  may  wish  to  give  the 
program  now.  They  are  'DIS',  'WHE',  'AIR',  and  'POS'. 

When  'DIS'  is  entered,  the  operator  has  asked  for  a graphic  display 
of  projected  storm  motion  on  the  remote  terminal.  The  progran  will  ask 
for  a prediction  time  interval,  a specified  range,  a graphic  overlay, 

and  the  operator  assigned  numbers  of  echoes  that  the  operator  wishes  to  ' 

see.  The  computer  will  type: 

'BTIME/ETIME/RANGE/OVERLAY/ECHO  NUMBERS' 

BTIME  and  ETIME  are  the  beginning  and  ending  times  for  the  prediction 
period  entered  in  the  same  four  digit  format  as  for  entering  echo  time 
data.  A useful  beginning  time  might  be  the  time  of  the  last  PPI  and  the 

ending  time  might  be  one  hour  later.  RANGE  lets  the  operator  choose 

between  a 200  or  a 400  km  range.  If  radar  data  were  previously  being 
transmitted  at  one  range,  the  operator  would  probably  also  want  the 
graphic  display  to  be  scaled  the  same.  There  are  two  choices  for  'OVERLAY'. 
They  are  the  Oklahoma  State  outline  entered  as  'STA'  and  the  low  level 
Victor  Airways  entered  as  'VIC'.  Esthetlcally , the  former  is  more  suited 
to  400  km  range,  while  the  latter  to  a 200  km  range.  If  that  entry  is 
left  blank,  no  overlay  is  produced  on  the  remote  radar  scope.  The  ECHO 
NUMBERS  are  those  assigned  by  the  operator. 

Another  option  is  to  ask  WHEn  the  centroid  of  a storm  will  be  near- 
est a given  point.  After  the  operator  has  entered  'WHE'  the  program 
will  type: 
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'ECHO  NO/AZM/RNG' 


The  operator  then  enters  the  appropriate  information,  where  AZM  and  RNG 
give  the  position  of  the  point  in  question,  not  the  centroid  of  the  echo. 
The  computer  will  return  the  distance  from  the  point  normal  to  the 
extrapolated  echo  path  and  the  standard  deviation  of  that  distance.  Time 
of  arrival  and  two  other  times,  one  before  and  one  after  the  predicted 
time,  are  also  computed.  These  times  are  t ± a^. 

Another  option  available  to  the  operator  is  to  ask  for  the  AIRports 
which  lie  in  the  echo  path.  The  echo  path  Is  considered  to  be  a cone 
whose  apex  is  the  centroid  position  at  the  time  of  the  last  observation. 
The  cone  expands  downstream  as  a function  of  time  (fig.  9)*  After  'AIR' 
has  been  entered  as  a conmand,  the  computer  responds  by  typing 

'BTIME/ETIME/STD  DEV/ECHO  NO. S' 


BTIME  and  ETIME  are  enterable  for.  DIS.  STD  DEV  is  entered  as  a whole 
number  (e.g.,  1,  2,  or  3).  The  expansion  rate  of  the  cone  is  determined 
by  Oj  (STD  DEV).  ECHO  NO.S  are  the  user  assigned  storm  numbers. 

The  computer  will  return  with  the  predicted  encounters  listed  by 
storm  in  the  order  of  arrival  time,  t^.  Also  given  is  the  distance  to 
the  echo  path  at  the  arrival  time  and  t^^  + and  t^  - 0^. 

One  last  question  the  operator  may  address  is,  "What  will  the  echo's 
position  be  at  a later  time?"  After  'POS'  has  been  entered,  the  computer 
will  type: 

'ECHO  NO/  HHMM' 


The  operator  enters  the  required  data  in  the  same  manner  as  under 
the  conmand  'ENT'.  After  the  above  Information  is  entered,  the  compjuter 


Figure  9.  Illustration  of 
cone  used  to  locate  air- 
ports in  an  echo's  path. 
ECHO  Integer  3 was  entered 

RATH  for  'STDDEV.' 


^Instead  of  a cone  for  determining  what  airports  may  be  affected,  the 
modified  parallelogram  of  the  preceding  section  could  be  used. 
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will  return  the  azimuth  and  range  of  the  centroid  for  the  time  given  and 
two  values  which  are  Oj  and  3o<j. 
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Because  computer  storage  restrictions  permit  only  ten  echoes  to  be 
tracked  simultaneously,  the  operator  will  need  to  DELete  echoes  from  time 
to  time.  The  operator  should  consider  first  those  echoes  which  are  no 
longer  being  actively  tracked. 

After  entry  of  'DEL',  the  computer  asks: 

'WHICH  ECHOES' 

To  this  the  operator  responds  by  entering  the  assigned  echo  numbers.  Before 
asking  for  a new  command,  the  computer  lists  the  active  echoes. 

On  occasion  the  operator  may  wish  to  enter  information  of  a textual 
nature.  This  might  include  severe  weather  events  associated  with  a partic- 
ular storm,  or  storm  tendencies  for  a new  operator  coming  on  duty.  The 
command  'IGN'  for  IGNnore  is  entered.  After  this,  a text  of  any  number  of 
lines  may  be  entered.  To  restore  the  program  to  an  operational  mode,  the 
operator  types  'KEY'  at  the  beginning  of  a new  line. 

To  terminate  the  program,  the  operator  types  'BYE'  for  a command. 

5.3  Some  Practical  Guidelines  for  Using  Echo  Prediction  Software 
5.3.1  Time  Weight  Constant  (TWC) 

At  the  time  the  echo  prediction  logic  was  developed,  it  was  recognized 
that  severe  storms,  especially  tornado  producers,  "turn  right"  as  they 
become  severe  (Newton  and  Fankhauser,  1964).  In  order  to  take  the  path 
curvature  into  account  in  predicting  future  echo  positions,  an  exponen- 
tially assigned  weighting  function  was  incorporated  into  the  computer  soft- 
ware. Mathematically,  the  function  is  W = where  k is  In  2,  TWC, 

the  value  entered  by  an  operator  and  At  the  time  interval  between  two 
observations.  The  rate  of  decay  is  a function  of  time  (fig.  10)  such  that 
when  the  time  elapsed  from  the  last  point  entered  is  equal  to  TWC,  the 
value  of  the  previous  point  is  decreased  by  half.  Although  the  prediction 
logic  was  operationally  tested  in  1969  and  1970,  TWC  was  always  made  large 
enough  (24  hours)  such  that  W*=l.  In  other  words,  all  centroid  positions 
carried  the  same  weight  when  fitting  LLS  equations. 

As  a first  step  in  testing  the  utility  of  changing  the  weighting 
function,  points  were  arbitrarily  entered  at  10  degree  intervals  at  a 
range  of  10  km  through  a 90  degree  sector.  The  sampling  interval  was 
entered  as  five  minutes.  The  results  for  k = 24  hrs,  30,  15  and  5 min- 
utes are  shown  in  Figure  lla-d.  An  examination  of  the  figures  shows  that 
the  greatest  improvement  in  following  the  data  trend  is  between  the  15  and 
5 minute  weighting  function.  Of  the  four  parameters,  direction,  speed, 
standard  deviation  of  distance  and  standard  deviation  of  time,  the  latter 
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showed  the  greatest  overall  Improve- 
ment— nearly  a factor  of  three  from 
.121  hours  to  .047  hours.  Least 
affected  was  the  speed. 

In  another  test,  real  data  from 
two  storms  were  used  to  determine  the 
effects  of  varying  the  weighting 
function.  One  storm  was  tracked 
between  1225Z  and  1310Z;  the  other 
between  1247Z  and  1310Z.  Two  time 
weight  constants  (TWO) — 30  and  5 min- 
utes— were  used.  Also,  because  the 
data  density  was  greater  than  other 
cases — 2 to  3 minute  Intervals — two 
passes  at  each  TWC  were  made,  one  at 
a 2 to  3 minute  spacing  and  the  other 
at  a 4 to  6 minute  spacing.  The 
results  are  shown  In  Table  1.  Con- 
trary to  the  results  In  the  first 
test,  there  is  little  if  any  improve- 
ment In  0(]  and  0(  between  a 30  and  a 
5 minute  TWC  for  2-3  minute  spacing. 
Using  5 minute  data  spacing  and  com- 
paring the  and  between  a 30  and 
5 minute  TWC  shows  some  overall 

Improvement.  By  far  the  greatest  improvement  Is  to  use  5 minute  data 
Instead  of  2-3  minute  data.  Several  sources  of  error  suggest  why  this 
Is  so.  One  Is  the  natural  variability  of  the  echo;  that  Is  the  random 
gain  or  loss  of  echo  due  to  small  scale  echo  changes.  Second,  radar 
system  fluctuations  of  1-2  dB  would  cause  small  scale  changes.  Third, 
and  perhaps  the  most  important.  Is  the  system  resolution.  Typically,  an 
echo  might  move  1-2  km  In  a 2 minute  period.  At  a range  of  100  km  this 
might  be  a change  of  centroid  location  of  1 degree  azimuth  or  1 km  range 
or  both.  Such  a motion  results  in  a very  noisy  path.  Barclay  and  Wilk 
(1970)  also  noted  erratic  echo  movement  from  centroid  positions  when 
using  data  of  similar  density.  In  the  remotlng  system,  one  must  also  con- 
sider the  time  element.  If  one  were  to  use  radar  data  directly  and  had 
the  time  of  each  radial,  the  centroid  time  would  be  precisely  known.  How- 
ever, the  time  of  the  data  displayed  on  the  remote  scope  is  truncated  to 
the  nearest  minute  so  that  the  time  the  echo  was  sampled  could  conceivably 
be  one  minute  off.  When  the  truncation  error  Is  a large  fraction  of  the 
AT,  the  projected  error  from  this  cause  tends  to  be  large. 


Figure  10.  Illustration  of  response 
curve  for  weighting  function,  W. 


Let  us  consider  now  the  speeds  and  directions  actually  computed  under 
the  different  conditions.  In  general,  the  TWC  had  more  effect  than  data 
density.  Echo  1,  with  TWC  equal  to  5 minutes,  accelerated  sharply  after 
1245  nearly  doubling  Its  speed  by  1306.  With  a TWC  of  30  minutes,  the 
acceleration  is  smoothed  considerably.  The  direction  also  shows  consider- 
ably more  variance  with  a 5 minute  TWC  than  with  a 30  minute  TWC.  Echo  2 
shows  basically  the  same  features  as  Echo  1. 


ELAPSED  TIME  (min) 


ELAPSED  TIME  (min) 


(c)  (d) 

Figure  11.  Illustration  of  t)ie  effect  of  varying  the  TWC  for  an  echo 
whose  track  is  curved  in  calculating  (a)  direction. 


(b)  speed,  (c)  a^,  and  (d)  a 
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Table  1.  Comparison  of  results  of  using  a TWC  of  both  30  and 
5 minutes  on  data  sampled  at  both  2-3  minutes  and 
4-6  minutes  for  tracking  two  echoes. 


WC  - 30  MIN  TWC  - 5 MIN 


ECHO  1 
2 MIN  DATA 


TIHE 

DIR 

SPD 

STDDS 

STD  TH 

DIR 

SPD 

STDDS 

STD  TM 

1230 

254.6 

23.8 

27.41 

.803 

247.3 

24.3 

27.41 

.803 

1236 

273.8 

37.6 

39.05 

.662 

272.1 

39.2 

40.64 

.614 

1240 

258.7 

34.6 

37.54 

.624 

252.2 

34.6 

38.67 

.583 

1245 

253.0 

36.1 

33.66 

.609 

248.4 

37.9 

34.59 

.579 

1251 

238.5 

39.7 

40.65 

.570 

227.9 

44.7 

39.87 

.546 

1255 

228.7 

43.2 

41.17 

.620 

216.0 

52.0 

38.63 

.594 

1300 

224.8 

50.5 

42.34 

.670 

219.2 

62.3 

41.19 

.612 

1306 

218.8 

55.4 

47.77 

.661 

210.3 

67.6 

45.60 

.602 

1310 

219.3 

55.1 

49.99 

.695 

217.0 

58.4 

47.67 

.655 

S HIN  data 

1230 

251.0 

24.0 

0.00 

.000 

251.0 

24.0 

0.00 

.000 

1236 

273.9 

36.4 

14.10 

.394 

276.2 

38.8 

14.10 

.394 

1240 

263.4 

34.1 

18.90 

.410 

259.4 

34.0 

19.20 

.429 

1245 

257.1 

36.7 

17.63 

.403 

252.2 

38.8 

17.52 

.422 

1251 

238.4 

39.3 

28.34 

.368 

223.6 

45.7 

27.20 

.386 

1255 

228.2 

43.6 

31.01 

.405 

214.7 

53.6 

27.09 

.415 

1300 

227.0 

48.7 

29.75 

.448 

222.3 

60.1 

28.25 

.397 

1306 

220.9 

53.0 

32.47 

.436 

211.2 

65.6 

30.34 

.377 

1310 

220.7 

53.6 

32.62 

.437 

216.6 

58.7 

31.25 

.412 

ECHO  2 
2 MIN  DATA 


1230 

- 

- . 

- 

- 

- 

- 

- 

1236 

- 

- 

- 

- 

- 

- 

. 

1240 

- 

- 

- 

- 

- 

- 

- 

_ 

1247 

346.1 

74.1 

0.00 

0.000 

346.1 

74.2 

0.00 

0.000 

1251 

270.7 

38.9 

62.34 

1.099 

268.6 

36.3 

62.06 

0.944 

1255 

247.0 

41.7 

57.18 

1.020 

241.1 

42.1 

56.50 

1.049 

1300 

233.5 

66.4 

56.79 

1.362 

230.9 

73.9 

53.82 

1.372 

1306 

228.3 

74.7 

52.68 

1.236 

226.3 

79.8 

48.81 

1.201 

1310 

229.6 

70.0 

48.35 

1.146 

230.4 

66.3 

45.93 

1.155 

5 MIN  DATA 


1230 

- 

- 

- 

- 

- 

- 

- 

- 

1236 

- 

- 

- 

- 

- 

- 

. 

1240 

- 

- 

- 

- 

- 

- 

- 

_ 

1247 

- 

- 

- 

- 

- 

- 

- 

- 

1251 

286.1 

33.2 

0.00 

0,000 

286.1 

33.2 

0.00 

0.000 

1255 

251.8 

38.2 

32.21 

0.075 

245.0 

41.3 

32.21 

0.075 

1300 

235.8 

65.0 

36.27 

0.912 

231.8 

78.6 

32.46 

0.863 

1306 

231.6 

72.0 

32.85 

0.817 

228.9 

79.1 

29.19 

0.774 

1310 

232.6 

66.9 

30.79 

0.803 

232.5 

64.5 

27.63 

0.772 
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Although  two  storms  are  admittedly  a small  sample,  our  experience  In 
a large  number  of  cases  makes  us  believe  that  these  results  apply  to  other 
storms  as  well.  Essentially,  what  Is  Indicated  Is  that  for  a five  minute 
forecast  of  echo  motion,  one  should  use  a five  minute  TWC,  with  five  minute 
data  resolution.  However,  for  a 30  to  60  minute  forecast,  a larger  TWC, 
such  as  30  minutes,  should  be  used.  Over  this  length  of  time,  trends  In 
the  overall  echo  motion  are  more  important  than  short  term  fluctuations 
which  should  be  smoothed  out.  Also,  since  warning  areas  are  mapped  out 
based  on  echo  speed  and  motion,  large  variance  from  one  time  to  the  next 
would  only  confuse  the  user  and  cause  the  warning  areas  to  shift  consid- 
erably from  one  prediction  PPI  to  another. 


5.3.2  Selection  of  Area  and  Intensity  Criteria 


In  choosing  the  threshold  area  and  intensity  the  user  is  hampered  by 
being  limited  to  ten  echoes.  However,  the  ability  to  assimilate  and  follow 
even  ten  Is  questionable.  On  a scope  containing  many  echoes,  therefore, 
the  user's  attention  should  be  directed  to  the  largest  and/or  strongest 
echoes. 


Barclay  and  Wilk  (1970)  determined  that  for  echo  extrapolation  using 
centroid  data,  threshold  values  of  10^  to  10^  mm^m“^  for  Isolated  storms 
and  10^  to  10^mm^“^  for  squall  lines  worked  best.  Based  on  his  own  expe- 
rience, the  author  believes  these  are  good  criteria.  With  higher  inten- 
sities, tracking  is  difficult  because  the  lifetimes  of  the  intense  cores 
are  short  and  in  squall  lines  the  probability  of  mergers  and  splits  also 
increases.  If  one  uses  too  low  an  intensity,  information  concerning  those 
areas  most  hazardous  to  aircraft  is  lost.  A general  rule-of-thumb  is  to 
use  the  lowest  Intensity  for  which  a discrete  echo  can  still  be  defined. 


In  a master's  thesis  in  1969,  R.  A.  Houze,  Jr.  defined  three  areal' 
sizes  associated  with  New  England  precipitation:  "synoptic  scale  areas" 

on  the  order  of  10^  - 10^  mi^  and  a duration  of  about  10  hours  to  pass  a 
point;  "mesoscale  areas"  on  the  order  of  lO^  - 10^  mi^  and  a duration  of 
about  one  hour;  and  "cells"  with  a 1-10  mi^  area  and  lasting  about  one  minute. 
The  scale  size  which  is  associated  with  severe  weather  and  with  which  we 
have  concerned  ourselves  in  this  report,  is  the  mesoscale.  The  synoptic 
scale  pattern  is  dependent  on  larger  atmospheric  disturbances  and  its 
movement  is  better  forecast  by  existing  NWS  software.  Also,  the  average 
radar  Intensity  does  not  constitute  a hazard  to  aviation  from  strong  shear 
or  turbulence.  Small,  Intense  cells,  on  the  other  hand,  are  of  too  short 
a duration  to  be  tracked  and  where  they  do  occur — imbedded  in  mesoscale 
systems — the  entire  area  should  be  avoided. 


The  average  area  of  eight  extremely  severe  isolated  storms  which  have 
occurred  in  Oklahoma  over  the  past  six  years  was  208  n mi^  (713  km^)  with  a 
range  of  122  to  427  n mi^.  Squall  lines  generally  are  about  three  times  the 
size  of  isolated  echoes,  although  cells  within  squall  lines  are  usually 
slightly  smaller  than  severe,  isolated  storms.  Since  the  areas  given 
above  are  for  total  storms,  the  actual  area  threshold  used  is  reduced  by  a 
factor  of  2 or  3 depending  on  the  intensity  threshold.  In  summary,  the 
area  threshold  ranges  between  100  and  1000  km^,  while  the  Intensity  varies 
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from  lO^mm^m'  ■3  to  10^mm^m~^.  Experience  guides  determination  of  the  best 
combination  for  any  given  storm  system. 

Once  the  criteria  are  established,  the  user  Is  well  advised  to  resist 
frequent  change.  The  purpose  in  et..:abllshing  criteria  is  to  give  the  user 
a history  for  summarizing  events  on  the  display.  Frequent  changes  will 
cause  loss  of  continuity  of  pattern. 

If  it  becomes  difficult  to  match  echoes  due  to  frequent  splitting  or 
merging,  however,  then  a lower  intensity  and  a larger  area  should  be 
selected.  Conversely,  if  persistent  significant  core  elements  are  being 
omitted,  then  higher  intensity  and  smaller  area  threshold  are  indicated. 

5.3.3  Splits  and  Mergers 

One  of  the  biggest  problems  in  using  echo  centroids  to  track  and 
extrapolate  future  positions  is  how  to  handle  splitting  or  merging  storms. 
For  example,  if  a storm  splits  into  two  distinct  cores,  should  you  regard 
one  of  the  cores  as  a continuation  of  the  old  echo  and  tag  the  other  core 
as  a new  echo  or  should  both  be  treated  as  new  echoes?  A similar  problem 
exists  with  mergers. 

The  best  procedure  the  author 
has  found  is  to  follow  trends  in 
area  and  centroid  position  for  each 
core  in  question.  Typically,  in  a 
five  minute  period,  an  echo  will 
move  from  3 to  6 km.  If  motion  is 
along  a radial,  then  the  centroid 
position  is  simply  a function  of 
range.  If,  however,  the  azimuth 
changes,  then  the  incremental  change 
is  range  dependent.  Figure  12 
shows  the  change  in  azimuth  as  a 
function  of  range  for  motion  tan- 
gent to  a circle  at  that  range. 

The  important  consideration  remains 
to  look  for  discontinuities  in 
either  range  or  azimuth.  A change 
in  range  10  km  and/or  10  degrees 
greater  than  expected  coupled  with 
a 20  percent  or  greater  change  in 
area,  should  be  considered  a new 
echo. 

6.  TEST  CASES 

Three  days  were  selected  for 
program  testing.  On  two  days, 

June  6 and  June  16,  1975,  squall 
lines  producing  damaging  surface 
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Figure  12.  Illustration  of  the  angular 
change  in  centroid  position  as  a 
function  of  range.  X is  the  dis- 
placement distance  of  the  centroid. 


winds  moved  across  the  State;  the  third  case  on  November  2,  1974,  produced 
localized  flooding  in  the  northwest  Oklahoma  City  with  cumulative  rainfall 
amounts  in  excess  of  five  inches. 

The  three  cases  above  were  chosen  because  of  their  danger  to  the  avia- 
tion community.  In  the  case  of  fast  moving  squalls,  the  Inherent  danger 
is  from  the  turbulence  and  high  winds  in  and  around  them.  In  the  case  of 
flooding,  aside  from  the  Intensity  of  storms,  the  sheer  persistence  of  heavy 
rain  in  one  location  implies  time  delays  for  air  traffic  leaving  and/or 
arriving,  or  for  circumnavigating  those  features.  Certainly,  it  is  of  the 
utmost  Importance  for  a pilot,  en  route,  to  know  of  adverse  conditions  which 
will  not  clear  his  destination  by  his  ETA. 

The  case  studies  were  made  to  develop  the  logic  to  generate  graphic 
forecasts  and  to  determine  the  feasibility  of  using  the  polar  coordinate 
display  as  a graphics  terminal.  The  analyses  excluded  quantitative  veri- 
fication of  the  predicted  echo  motion.  The  cases  were  analyzed  using 
variations  of  the  echo  tracking  logic,  involving  different  cutoff  limits 
for  magnitudes  of  a^j  and  Oj.^in  expanding  the  warning  areas. 

In  simulating  real  time  operations,  radar  data  archived  on  magnetic 
tape  were  reformatted  to  "look"  like  data  received  by  the  remote  radar 
display  program.  A second  tape  was  generated  by  the  echo  prediction  logic 
which  contained  the  graphics  information,  also  in  the  format  for  the 
remoting  system.  The  graphics  tape  was  then  played  through  a tape  drive* 
and  transmitted  to  the  receiver  where  the  graphics  were  photographed. 

There  are  three  components  to  the  graphics:  individual  warning  areas 

for  one  or  more  storms;  a contour  at  a specified  intensity  threshold  for 
each  echo;  and  a background  reference  field,  showing  the  Victor  Airways  or 
the  state  outline.  Each  component  is  coded  at  one  of  three 'intensity' 
values  to  produce  differential  gray  shading.  Warning  areas  were  coded  in 
the  seventh  intensity  level  (not  normally  used  in  the  real  time  echo  dis- 
play) and  displayed  at  the  brightest  gray-shade  level;  echo  contours  were 
coded  in  the  second  intensity  level  and  displayed  at  the  medium  gray-shade 
level;  and  the  Victor  Airways  were  coded  in  the  first  intensity  level  and 
displayed  at  the  lightest  gray  shade  level. 

After  reviewing  the  radar  data  on  the  remote  scope,  threshold  criteria 
were  selected  and  an  initial  run  was  made  using  the  echo  prediction  logic. 
The  area  and  centroid  data  from  this  run  were  matched  and  a second  run  was 
made  using  these  results  to  generate  the  warning  areas.  The  following 
cases  Illustrate  the  logic  and  generation  of  the  graphics. 

Case  1,  June  6,  1975 

A stationary  front  was  indicated  on  the  12Z  surface  map  (fig.  13)  as 
extending  from  the  Oklahoma  Texas  panhandles  northeastward  across  the 


*As  part  of  the  specifications,  a half-inch  magnetic  tape  Interface  at 
the  transmitter  was  provided. 
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southern  third  of  Kansas;  a cold  front  was  entering  Nebraska  to  the  north. 
The  upper  level  flow  at  500  mb  (fig.  14)  was  from  the  west-northwest. 

Storms  broke  out  along  the  stationary  front  in  Kansas  in  the  early  after- 
noon and  generally  moved  southeast  in  the  same  direction  as  the  flow 
aloft.  Abundant  moisture  was  available  in  central  Oklahoma  (18-20  g kg”"^) 
and  the  air  was  potentially  unstable.  Surface  winds  were  from  the  south- 
southeast  at  15-20  knots.  By  2000  CST,  a well-developed  dry  line  oriented 
north-south  existed  in  the  Texas  and  Oklahoma  panhandles  as  shown  by  the 
surface  analysis  with  streamlines  in  Figure  15.  Cold  air  produced  by  the 
storms  in  western  Kansas  had  produced  a pseudo  cold  front  as  well.  There- 
fore, an  area  of  strong  convergence  formed  in  northwest  Oklahoma  (fig.  16).. 
As  the  Kansas  storms  were  following  the  same  track,  growth  was  favored  on 
the  southern  flank  of  these  storms  and  a line  formed  propagating  southward. 

As  shown  in  Figure  17,  the  closest  echoes  at  1600  were  about  200  km 
away  along  the  Kansas-Oklahoma  border.  Because  the  motion  at  that  time 
was  east-southeast,  the  echoes  were  not  expected  to  move  within  Doppler 
radar  range.  Hence,  operations  ceased  and,  between  1625  and  1830  CST, 
no  radar  data  were  collected.  The  radar,  monitored  at  1830  CST,  indicated 
a large  echo  was  still  beyond  Doppler  radar  range.  Within  two  hours, 
significant  changes  occurred.  A line  of  storms  developed  in  the  northwest 
quadrant  and  proceeded  to  move  southward  into  central  Oklahoma. 


For  experimental  development  of  the  echo  prediction  and  display  logic, 
analysis  was  begun  at  2040  CST  when  the  echoes  were  within  120  km.  A con- 
densation*^ echo  positions  and  motions  is  presented  in  Table  2.  The  com- 
plete program  command  structure  with  annotation  is  Included  as  Appendix  1. 
For  this  storm  period,  the  size  and  intensity  criteria  selected  were 
250  km^  and  code  switch  4 (corresponding  to  40  dBZ) , respectively.  Data 
at  0 degree  elevation  were  available  at  five  minute  intervals.  The  time, 
weight  constant  was  set  at  2400  (24  hr)  and  the  ground  clutter  distance 
at  20  km. 


Figure  13.  June  6,  1975,  12Z,  sur-  Figure  14.  June  6,  1975,  12Z, 500  mb 
face  analysis  adapted  from  DOC,  analysis  adapted  from  DOC,  NOAA, 

NOAA,  EDS  daily  weather  maps,  EDS  daily  weather  maps.  Weekly 

Weekly  Series.  Series. 


28 


SBbct 

m 

Inn 

IBfi 

w^pBB 

]|^B&| 

I^H 

V^BHKr 

IBrrlffi 

h^^SHI 

nSB^kmI 

Sjm^^H 

||KnH|l 

I 


Figure  17.  WSR-57  radar  PPI, 
400  km  range,  100  km  range 
marks,  1600  CST,  June  6, 
1975. 


Gray  shades  for  each  intensity  level  were  set  at  1220333  and  resulted 
in  echoes  being  contoured  at  cancellation  level  on  the  remote  scope. 

At  2040  CST,  as  shown  in  Figure  18,  two  large  cells  were  located  by 
the  computer  search  at  azimuth  29®,  range  136  km  and  at  azimuth  335®,  range 
157  km.  After  the  second  PPI  at  2045  CST,  the  same  two  cells  were  matched 
manually  to  the  previous  cells.  They  were  labeled  echo  1 and  echo  2 and 
their  respective  times  and  positions  were  entered  for  tracking.  Signifi- 
cantly different  motions  were  obtained:  echo  1 moved  from  west-northwest 

(293®)  at  93  km  hr~l  while  echo  2 moved  from  north-northwest  (335®)  at 
52  km  hr“l.  Warning  graphics  were  then  displayed  for  these  echoes  (fig.  19). 
The  starting  position  for  each  warning  area  is  15  minutes  after  the  last 
observation,  the  ending  position  is  one  hour  and  15  minutes  later.  All 
subsequent  warning  areas  were  also  for  one  hour  duration  starting  15  min- 
utes downstream.  A circle  was  drawn  on  the  graphics  display  after  the  PPI 
at  2055  (fig.  20)  to  represent  the  echo  shape.  The  program  does  not  save 
echo  shape  information  for  more  than  one  PPI.  Whenever  a graphic  display 
of  a warning  area  is  requested  for  an  echo  for  which  no  entry  was  made  for 
the  last  PPI,  a circle  with  a 10  km  radius  automatically  replaces  the 
initial  shape  function. 

A third  echo  located  at  azimuth  at  60®  range  195  km,  at  2050  showed 
no  movement  by  2055  CST.  The  resulting  warning  area  was  a line  drawn  along 
the  major  axis  of  the  echo  (fig.  20).  With  the  addition  of  a third  point 
at  2100  CST,  an  exorbitantly  large  Of.  was  calculated.  When  the  prediction 
field  was  first  displayed,  the  large  value  of  a|-  distorted  the  warning  area. 
Therefore,  a "cutoff"  limit  was  introduced  such  that  whenever  at  exceeded 
the  prediction  interval,  tg  - t^j,  0t  was  set  to  (tg  - t|,)/2.  Later,  a 
"cutoff"  limit  was  also  established  for  a^  restricting  it  to  the  length  of 
an  echo's  axis  used  in  defining  the  warning  area.  The  choice  of  these 
limits  was  purely  arbitrary  and  may  need  further  adjustment  based  on  a 
large  sample  of  data.  Obviously,  many  tracks,  based  on  the  first  few  obser- 
vations, will  show  considerable  scatter.  As  additional  locations  are  added 
over  a larger  period  of  time,  a better  estimate  of  the  mean  velocity  will 
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Table  2.  Echo  centroid  positions,  direction  and  speed  predicted, 
and  the  RMSE  values  for  June  6,  1975. 

ECHO  TIME  AZIMUTH  RANGE  DIRECTION  SPEED  ( 

NO,  (CST)  (DEGREES)  (KM)  (DEGREES)  (KM  HR"1)  (KM)  (H( 


I 

292.5 

92.9 

- 

- 

335.0 

51.5 

- 

- 

283.1 

78.7 

14.36 

0,1 

319.0 

58.2 

17.90 

0.0 

- 

— 

— 

— 

276.9 

75.9 

16.73 

0.1 

360.0 

0.0 

- 

- 

277.1 

73.1 

15.58 

0.1 

334.2 

64.3 

23.26 

0.1 

241.0 

6.1 

6.06 

* 

- 

- 

- 

- 

- 

- 

- 

- 

331.8 

62.8 

24.05 

0.1 

138.9 

12.2 

20.27 

* 

356.8 

49.0 

- 

- 

9.8 

76.7 

- 

- 

339.1 

66.1 

40.99 

0.2 

139.7 

20.2 

18.42 

* 

345.0 

55.0 

12.69 

0.1 

272.1 

79.2 

15.85 

0.1 

341.2 

68.1 

38.06 

0.2 

137.7 

8.1 

16.84 

* 

340.6 

57.4 

12.34 

0.1 

277.9 

80.7 

44.36 

0.1 

340.5 

66.6 

37.05 

0.2 

153.9 

7.1 

17.39 

* 

329.4 

57.6 

25.53 

0.1 

- 

- 

— 

- 

279.9 

82.1 

41.50 

0.1 

341.5 

69.1 

36.33 

0.3 

178.8 

7.0 

18.34 

- 

312.9 

66.7 

50.69 

0.3 

305.0 

61.0 

- 

- 

- 

— 

— 

- 

278.3 

81.8 

46.15 

0.1 

341.7 

69.5 

34.56 

0.3 

310.0 

68.5 

47.67 

0.3 

318.3 

55.8 

14.75 

0.1 

71.1 

46.1 

- 

- 

ting 

the  graphics  the 

values  for 

at  excee 

e observation 
rs. 

for  cell 

3 at  2055 

CST,  at 

be  made.  At  2105  CST  three  cells  were  associated  with  the  squall  line 
(fig.  21).  Generally,  the  core  in  the  middle  cell  in  the  line  was  too 
small  to  be  tracked. 

No  major  changes  occurred  in  the  pattern  but  several  more  PPI's  and 
graphics  (figs.  22-27)  on  this  data  are  presented  as  examples.  The  last 
PPI  at  2130  CST  and  graphic  warning  areas  which  were  displayed  at  a 400  km 
range  with  the  State  outline.  The  last  radar  PPI  at  2249  CST  (fig.  28) 
(about  the  ending  time  of  the  last  warning  area)  shows  that  the  line  had 
moved  as  expected. 
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Figure  18.  (left,  above)  Remote 
radar  display  PPI,  200  km  range, 
2040  CST,  June  6,  1975. 

Figure  19.  (above)  Remote  radar 
display  with  computer  generated 
warning  areas,  echo  contours  and 
Victor  airways,  June  6,  1975. 
Prediction  period  is  2100  to 
2200  CST. 

Figure  20.  (left)  Remote  radar 
display  with  computer  generated 
warning  areas,  echo  contours  and 
Victor  airways,  June  6,  1975. 
Prediction  period  is  2110  to 
2210  CST. 
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Figure  21.  Remote  radar  display 
PPI,  200  km  range,  2105  CST, 
June  6,  1975. 


Figure  22.  Remote  radar  display 
PPI,  200  km  range,  2115  CST, 
June  6,  1975. 


Figure  23.  Remote  radar  display 
with  computer  generated  warning 
areas,  echo  contours  and  Victor 
airways,  June  6,  1975.  Predic- 
tion period  is  2130  to  2230  CST 


Figure  24.  Remote  radar  display 
PPI,  200  km  range,  2120  CST. 
June  6,  1975. 


Figure  25.  Remote  radar  display 
computer  generated  warning 
areas,  echo  contours  and  Victor 
airways,  June  6,  1975.  Predic- 
tion period  is  2135  to  2235  CST 


Figure  26.  Remote  radar  display 
PPl,  200  km  range,  2130  CST, 
June  6,  1975. 
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Figure  27.  Remote  radar  display 
with  computer  generated  warning 
areas,  echo  contours  and  Okla- 
homa State  outline,  400  km  range 
June  6,  1975.  Prediction  time 
is  2145  to  2245. 


Figure  28.  WSR-57  radar  PPI, 

200  km  range,  40  km  range  marks 
2249  CST,  June  6,  1975. 


Case  2,  June  16,  1975  | 

j 

On  this  day,  moisture  from  the  Gulf  of  Mexico  and  drier  continental 
air  were  separated  by  a stationary  surface  front  extending  from  the  Okla- 
homa panhandle  across  southern  Oklahoma  and  northern  Arkansas  (fig.  29). 

During  the  next  24  hours  this  boundary  moved  northeastward  as  a warm  front 
ahead  of  a cold  front  approaching  from  the  northwest. 

The  flow  at  500  mb  was  essentially  zonal  at  20-30  kt  (fig.  20).  Mix- 
ing ratio  values  during  the  afternoon  of  the  16th  were  16-18  g kg“^,  and 

the  air  mass  was  potentially  unstable.  When  convective  temperature  was 
reached  after  1500  CST,  numerous  showers  and  thunderstorms  developed  over 
Kansas,  Oklahoma,  north  Texas  and  New  Mexico.  The  sequence  of  satellite 
photographs  and  WSR-57  PPI  displays  (figs.  32a,  b,  c,  and  33a,  b,  c)  trace 
the  growth  and  movement  of  these  storms  dtiring  the  afternoon. 

] 

Until  1830  CST,  activity  was  beyond  200  km,  and  the  radar  integrator  1 

was  range  delayed  (fig.  32)  to  provide  1 km  resolution  data  between  200  ! 

and  400  km  (Sirmans  and  Doviak,  1973).  After  1830  CST,  data  collection  ■ 

was  normal  with  data  recorded  at  five  minute  intervals.  ! 

The  area  of  thunderstorms  visible  over  northern  Texas  on  the  GOES 
satellite  photos  appeared  as  a line  of  moderate  to  intense  echoes  on  the 
WSR-57  display.  These  storms  intensified  and  moved  rapidly  across  central 
Oklahoma  as  an  organized  squall  line,  preceded  by  a strong  damaging  gust 
front.  The  NSSFC  issued  a tornado  watch  at  1815  CST  valid  from  2000  CST 
p to  0200  CST  to  cover  western  Oklahoma,  but  it  lagged  spatially  behind  the 

storms  because  of  their  unexpected  acceleration.  The  steering  level  wind 
|i  velocity  was  210  degrees  at  25  kts  (46  km  hr~^) . However,  line  motion  was 

much  faster  from  270  degrees  at  45  kts  (80  km  hr~l) . Figures  34  and  35, 

30  minutes  before  tracking  began,  shows  the  storm  influenced  surface  stream- 
i;  line  and  convergence  patterns,  respectively. 


Figure  29.  June  16,  1975,  12Z  sur- 
face analysis  adapted  from  DOC, 
NOAA,  EDS  dally  weather  maps. 
Weekly  Series. 


Figure  30.  June  16,  1975,  12Z  500  mb 
analysis  adapted  from  DOC,  NOAA, 
EDS  dally  weather  maps.  Weekly 
Series. 
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Figure  31.  SMS-3  satellite  pic 
tures,  June  16,  1975: 

(a)  1545  CST  (2145Z) 

(b)  1615  CST  (2215Z) 

(c)  1745  CST  (2345Z) 


(a) 


Figure  32.  WSR-57  radar  PPI, 
400  km  range  with  200  km 
range  delay  to  first  gate, 
40  km  range  marks,  June  16, 
1975: 

(a)  1545  CSX 

(b)  1615  CSX 

(c)  1745  CSX 
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Figure  33.  Subsynoptic  surface  data  provided  by  NWS,  FAA,  and 


NSSL  stations,  1800  CST,  June  16,  1975,  with  stream- 
lines superimposed. 


Testing  of  the  echo  prediction 


06-16-7S  1800CST 


Figure  34.  Divergence  field  derived 
from  surface  data,  1800  CST, 

June  16,  1975. 


logic  began  with  the  1830  CST  observa- 
tion. Using  intensity  and  area  thresh- 
olds of  40  dBZ  and  250  km^,  respectively, 
the  storm  centroids  were  determined  at 
5 minute  Intervals.  The  four-level 
intensity  code  for  the  remote  display 
was  set  to  light,  medium,  cancel  and 
bright  (1220333).  The  time  weight  con-  , 

stant  (TWC)  was  2400,  and  the  minimum 
range  (GCD)  was  20  km.  All  warning  areas  i 

were  for  one  hour  Interval  starting  with  " ‘ 

the  time  of  the  last  observation.  i 


After  two  PPI's  (figs.  35  and  36)  ' C 

centroid  positions  for  cell  1 were 
entered  and  an  estimate  of  the  echo's 
motion  (from  337  degrees  at  44  km  hr"^-)  ^ 

was  obtained  (Table  3) . Figure  37  shows 
the  subsequent  warning  area  based  on 
two  centroid  positions.  t-'*' 


i...  ■ 
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Figure  35.  Remote  radar  display 
PPI,  200  km  range,  1830  CST, 
June  16,  1975. 
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Figure  37.  Remote  radar  display 
with  computer  generated  warning 
area,  echo  contour  and  Victor  air- 
ways, June  16,  1975.  Prediction 
period  is  1835  to  1935  CST. 


Figure  36.  Remote  radar  display 
PPI,  200  km  range,  1835  CST, 
June  16,  1975. 


It  is  important  to  note  here  that 
as  an  echo  enters  the  scope,  its  cen- 
troid position  will  be  Influenced  by 
the  area  change.  As  a result,  using 
centroid  positions  to  calculate  echo 
motion  will  underestimate  true  echo 
speed.  Likewise,  direction  of  motion 
will  be  affected.  If  the  Increase  in 
each  area  occurs  only  at  one  end  of  a 
line,  the  calculated  direction  of 
motion  will  be  pulled  towards  that  end. 
Until  such  time  as  an  echo  or  squall 
line  has  fully  entered  the  scope,  it 
is  probably  better  to  track  a point 
along  its  leading  edge.  (That  was  not 
done  for  this  case.) 

From  Table  3,  an  increase  in  line 
speed  is  evident  through  the  first  half 
hour  of  tracking.  Because  of  the  bound- 
ary problem  just  described, large  values 
of  and  aj  resulted.  An  example  of 
a PPI  and  its  associated  warning  area 
at  1900  CST  are  shown  in  Figures  38  and 
39. 


Table  3.  Echo  centroid  positions,  direction  and  speed  predicted, 
and  the  RMSE  values  for  June  16,  1975. 


ECHO 

NO. 

TIME 

(CST) 

AZIMUTH 

(DEGREES) 

RANGE 

(KM) 

DIRECTION 

(DEGREES) 

SPEED 
(KM  HR“1) 

(KM) 

(HOUR) 

1 

1830 

280. 

188. 

_ 

1 

1835 

279. 

186. 

338.0 

45.9 

- 

- 

1 

1840 

278. 

183. 

331.3 

49.3 

6.59 

0.070 

1 

1845 

279. 

179. 

302.5 

39.2 

31.11 

0.516 

1 

1850 

276. 

174. 

314.0 

52.1 

43.65 

0.843 

1 

1855 

273. 

171. 

324.0 

63.3 

46.94 

0.845 

1 

1900 

271. 

165. 

325.4 

71.5 

43.74 

0.812 

2 

1905 

262. 

170. 

- 

- 

- 

— 

3 

1905 

287. 

151. 

- 

- 

- 

- 

2 

1910 

261. 

167. 

305.9 

55.2 

- 

- 

3 

1910 

285. 

142. 

315.6 

135.9 

- 

- 

2 

1915 

261. 

159. 

276.1 

69.7 

39.15 

0.134 

3 

1915 

287. 

133. 

287.0 

110.4 

59.72 

0.299 

1 

1920 

264. 

146. 

321.3 

77.7 

41.41 

0.761 

1 

1925 

263. 

140. 

318.7 

78.4 

41.81 

0.720 

1 

1930 

263. 

132. 

315.3 

77.8 

46.93 

0.690 

1 

1935 

260. 

126. 

313.3 

78.7 

44.85 

0.668 

1 

1940 

258. 

120. 

311.7 

79.3 

43.72 

0.640 

1 

1945 

259. 

111. 

309.0 

79.1 

50.48 

0.626 

4 

1950 

263. 

no. 

- 

- 

- 

- 

4 

1955 

264. 

108. 

219.9 

36.2 

- 

- 

4 

2000 

267. 

97. 

235.9 

91.2 

28.38 

1.627 

5 

2005 

263. 

82. 

- 

- 

- 

- 

5 

2010 

263. 

74. 

263.0 

105.0 

- 

- 

5 

2015 

263. 

64. 

263.0 

110.4 

0.01 

0.082 

6 

2020 

285. 

51. 

- 

- 

- 

- 

6 

2025 

284. 

46. 

294.1 

66.6 

- 

- 

6 

2030 

287. 

41. 

276.9 

62.1 

20.79 

0.099 

% 

I 
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Figure  38.  Remote  radar  display 
PPI,  200  km  range,  1900  CST, 
June  16,  1975. 


Figure  39.  Remote  radar  display 
with  computer  generated  warning 
area,  echo  contour  and  Victor  air- 
ways, June  16,  1975.  Prediction 
period  is  1900  to  2000  CST. 


At  1905  CST,  the  program  logic  isolated  two  discrete  cells  within  the 
line  at  intensity  40  dBZ  where  there  had  been  only  one  five  minutes  before. 
The  split-off  of  the  smaller  of  the  two  cells  caused  the  centroid  of  the 
larger  cell  to  be  shifted  nine  degrees  in  azimuth  and  moved  back  5 km.  After 
the  split,  the  two  cells  were  identified  as  cell  2 and  cell  3 and  tracked 
for  the  next  ten  minutes.  The  PPI  at  1915  CST  (fig.  40)  and  the  graphics 
for  that  time  (fig.  41)  indicate  quite  different  speeds  for  each  cell  (see 
also  Table  3) . 

At  1920  CST,  cells  2 and  3 merged.  Since  the  resulting  centroid  position 
was  consistent  with  the  extrapolated  position  for  cell  1,  the  merge  was 
reassigned  as  cell  1 (fig.  42).  It  was  traced  until  1945  CST.  Two  radar 
PPIs  and  their  respective  warning  areas  (figs.  43-46,  at  1930  and  1945  CST) 
show  the  line's  movement  during  this  period. 

After  1945  CST,  in  the  initial  pass  through  the  data,  the  large  core 
which  had  been  cell  1,  fragmented  and  attempts  to  match  the  new  cells  proved 
futile.  Therefore,  a second  pass  was  made  with  new  threshold  criteria  of 
1000  km2  and  intensity  level  3.  The  gray  shading  was  recoded  at  1203333  to 
reflect  the  decrease  in  the  intensity  threshold. 

Cell  4,  tracked  from  1950  to  2000  CST,  indicated  motion  in  a north- 
easterly direction  due,  primarily,  to  growth  on  the  north  end  of  the  line. 
Motion  had  been  southeasterly  previously.  Shown  in  Figures  47  and  48  is  a 
radar  PPI  and  graphic  warning  for  2000  CST. 
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Figure  40.  Remote  radar  display 
PPI,  200  km  range,  1915  CST, 
June  16,  1975. 


Figure  41.  Remote  radar  display 
with  computer  generated  warning 
areas,  echo  contours  and  Victor 
airways,  June  16,  1975.  Predic- 
tion period  is  1915  to  2015  CST. 


Figure  42.  Remote  radar  display  Figure  43.  Remote  radar  display 

PPI,  200  km  range,  1920  CST,  PPI,  200  km  range,  1930  CST, 

June  16,  1975.  June  16,  1975. 
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Figure  A6.  Remote  radar  display 
with  computer  generated  warning 
area,  echo  contour  and  Victor 
airways,  June  16,  1975.  Predic- 
tion period  Is  1945  to  2045  CST. 


Figure  47.  Remote  radar  display 
PPI,  200  km  range,  2000  CST, 
June  16,  1975. 


At  2005  CST,  due  to  a substantial 
change  in  area  and  range,  the  largest 
echo  was  relabeled  as  cell  5.  Like 
cell  4,  cell  5 was  also  tracked  for  only 
15  minutes.  As  shown  in  Figures  49  and 
50,  and  Table  3,  the  warning  area  and 
echo  motion  after  10  minutes  reflect 
the  line's  overall  speed  and  direction 
of  motion  better  than  any  earlier  time. 

A substantial  change  in  area  and 
centroid  position  again  occurred  about 
2020  CST  with  the  additional  growth  on 
the  north  end  of  the  line.  Cell  5 was 
dropped.  The  larger  echo,  renumbered 
as  cell  6,  moved  slower,  although  the 
direction  was  the  same.  The  PPI  and 
warning  area  (fig.  51  and  52)  are  shown 
for  2030  CST  when  tracking  ceased. 

Probably  a longer  sampling  inter- 
val should  have  been  used  for  this  case 
(e.g.,  15  minutes).  However,  several 
factors  have  to  be  considered.  The 
optimum  sampling  Interval  determined 
by  Wllk  and  Gray  (1970)  was  45  minutes. 
Obviously,  one  can't  wait  that  long  before  making  a prediction.  Also,  the 
lifetime  of  the  storm  may  be  less  than  an  hour.  Conversely,  if  one  samples 
too  frequently,  lack  of  spatial  and  temporal  resolution  will  produce  ficti- 
tiously large  errors. 

Another  dilemma  arises  when  matching  echo  manually.  The  more  frequently 
one  samples,  the  easier  it  is  to  follow  echo  motion  and  to  account  for  splits 
and  merges.  (On  some  occasions  it  was  difficult  to  match  five  minute  data.) 
However,  an  unwarranted  amount  of  time  may  be  spent  tracking  snail  cells  of 
short  duration  which  neither  produce  severe  weather  nor  have  sufficient  pre- 
dictability to  provide  meaningful  extrapolations. 

Case  3,  November  2,  1974 

The  storms  developed  early  in  the  morning  of  November  2.  On  the  previous 
day  a stationary  front  extended  across  southern  Texas  and  northern  Louisiana 
(fig.  53).  During  the  day,  it  evolved  into  a warm  front  which  moved  into 
southern  Oklahoma  (fig.  54).  The  position  of  the  front  changed  little  during 
the  afternoon  and  evening  of  November  2 (fig.  55)  as  an  upper  level  low 
developed  over  California,  and  maintained  a stationary  pattern  of  southwest- 
erly flow  aloft  over  Oklahoma. 

At  500  mb  (fig.  56)  a broad  trough  covered  most  of  the  United  States 
with  two  low  pressure  centers,  one  centered  over  northeast  Wyoming,  the  other 
over  northern  California  and  Nevada.  Between  November  1 and  3 the  northern 
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Figure  48.  Remote  radar  display 
with  computer  generated  warning 
area,  echo  contour  and  Victor 
airways,  June  16,  1975.  Predic- 
tion period  is  2000  to  2100  CST. 
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Figure  49.  Remote  radar  display 
PPI,  200  km  range,  2015  CST, 
June  16,  1975. 


Figure  50.  Remote  radar  display 
with  computer  generated  warning 
area,  echo  contour  and  Victor 
airways,  June  16,  1975.  Predic- 
tion period  is  2015  to  2115  CST 


Figure  51.  Remote  radar  display 
PPI,  200  km  range,  2030  CST, 
June  16,  1975. 


Figure  52.  Remote  radar  display 
with  computer  generated  warning 
area,  echo  contour  and  Victor 
airways,  June  16,  1975.  Predic- 
tion period  is  2030  to  2130  CST. 


i|  Figure  53.  November  1,  1974,  12Z  Figure  54.  November  2,  1974,  12Z 

tj  surface  analysis  adapted  from  surface  analysis  adapted  from 

' ' DOC,  NOAA,  EDS  daily  weather  DOC,  NOAA,  EDS  daily  weather 

maps.  Weekly  Series.  maps.  Weekly  Series. 


low  migrated  northeastward  into  Canada; 
the  other  moved  southward  as  a separate 
closed  low.  This  resulted  in  a shift 
of  the  primary  trough  to  a northeast- 
southeast  orientation  and  caused  the  jet 
stream  to  retrograde  westward  (figs.  57 
and  58).  Since  there  was  insufficient 
frontal  lift  and  no  surface  heating, 
the  triggering  mechanism  for  the  early 
morning  storms  on  November  2 was  prob- 
ably a short  wave,  produced  in  the  lee 
of  the  Rocky  Mountains.  Over  central 
Oklahoma,  mixing  ratio  values  were 
12-14  g kg~l,  and  low  clouds  and  high 
relative  humidity  were  widespread 
(fig.  59).  At  0600  CST,  Oklahoma  City 
reported  a ceiling  of  200  ft,  Hobart 
and  Cllnton-Sherman,  150  km  to  the  west, 
200-300  ft  ceilings,  and  Ardmore,  600  ft. 
At  stations  north  of  Oklahoma  City, 
rain  and  fog  were  reported. 

Surface  winds  were  generally  light  5-10  kt  (fig.  59)  which  was  unrepre- 
sentative of  the  mean  flow.  Surface  mixing  was  not  occurring,  as  indicated 
by  the  winds  recorded  at  444  m on  the  WKY  tower  (Goff  and  Zittel,  1974), 
which  were  20-30  kts  form  the  south-southeast.  This  flow  more  accurately 
reflects  the  true  inflow  into  the  storms. 

Testing  of  the  echo  prediction  logic  began  with  the  1225Z  observation 
, and  continued  until  1210Z.  In  this  case,  digital  radar  data  were  available 

I every  two  to  three  minutes.  However,  graphic  displays  were  produced  as 
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Figure  55.  November  3,  1974,  12Z 
surface  analysis  adapted  from 
DOC,  NOAA,  EDS  dally  weather 
maps.  Weekly  Series. 


Figure  56.  November  1,  1974,  12Z  Figure  57.  November  2,  1974,  12Z 

500  mb  analysis  adapted  from  500  mb  analysis  adapted  from 

DOC,  NOAA,  EDS  daily  weather  DOC,  NOAA,  EDS  dally  weather 

maps.  Weekly  Series.  maps.  Weekly  Series. 


before  at  approximately  five  minute 
Intervals.  The  thresholds  for  inten- 
sity and  area  were  40  dBZ  (intensity 
switch  4)  and  150  km^,  respectively. 
The  coding  of  gray  shades  for  photog- 
raphy was  1220333.  The  time  weight 
constant  was  2400  and  the  minimum 
range  (to  omit  ground  clutter)  was 
20  km. 

Because  the  radar  pattern  changed 
little  during  the  45  minutes  of  track- 
ing, only  selected  radar  PPIs  and 
graphic  PPIs  are  shown  for  this  case 
(figs.  60  through  69). 

Warning  areas  were  drawn  for  a 
one-hour  prediction  interval  starting 
from  the  time  of  the  last  observation. 


Figure  58.  November  3,  1974,  12Z 
500  mb  analysis  adapted  from 
DOC,  NOAA,  EDS  dally  weather 
maps.  Weekly  Series. 


In  the  first  PPI  at  1225Z,  two 
cells  were  Isolated  and  labeled  cell  1 

and  cell  2,  respectively.  Cell  1,  at  250  degrees  azimuth  and  76  km  range, 
was  tracked  for  the  entire  45  minute  period.  As  shown  in  Table  4,  excluding 
the  first  prediction,  cell  motion  was  southeasterly  gradually  shifting  to  the 
northeast  and  accelerating.  By  1310Z,  cell  motion  was  from  220  degrees  azi- 
muth at  53  km  hr~l. 


Cell  2 split  after  the  first  PPI  becoming  two  cells.  These  merged  later 
at  1234Z  and  tracking  of  cell  2 was  resumed  (table  4).  This  cell,  imbedded 
in  the  line  northwest  of  Oklahoma  City,  moved  from  240  degrees  azimuth  about 
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Figure  59.  Subsynoptic  surface  data  provided  by  NWS,  FAA  and  NSSL 
stations,  1200Z.  November  2,  1974,  with  streamlines 
superimposed. 

60  km  hr  At  124SZ  when  a sudden  Increase  in  the  line's  area  occurred, 
cell  2 was  dropped. 

After  cell  2 split  at  1228Z,  one  of  the  new  cells  was  labelled  cell  3 
and  tracked  until  1234Z.  The  other  cell  was  never  tracked  because  little 
motion  was  shown. 

Cell  5 was  isolated  at  the  north  end  of  the  line  at  1230Z  and  tracked 
for  15  minutes.  During  that  time  it  moved  from  a more  southerly  direction 
(214  degrees)  than  the  other  cells. 

At  1245Z,  a much  larger  cell  was  isolated  due  to  the  increase  in  inten- 
sity and  size.  This  cell,  cell  6,  was  tracked  until  testing  ceased  at  1310Z. 

As  has  been  pointed  out  in  previous  sections,  the  large  values  for 
and  Ojj  (table  4)  are  in  part  the  result  of  too  frequent  sampling.  In  some 
cases,  the  graphics  examples  do  not  show  the  high  values  of  and  because 
they  exceeded  the  predetermined  limits,  as  mentioned  in  case  1. 
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Table  4.  Echo  centroid  positions,  direction  and  speed  predicted, 
and  the  RMSE  values  for  November  2,  1974. 
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ECHO 

NO. 

TIME 

(Z) 

AZIMUTH 

(DEGREES) 

RANGE 

(KM) 

DIRECTION 

(DEGREES) 

SPEED 
(KM  HR-1) 

(KM) 

<^t 

(HOUR) 

1 

1225 

251. 

76. 

2 

1225 

311. 

95. 

- 

- 

- 

- 

1 

1228 

250. 

75. 

303.3 

33.1 

- 

- 

3 

1228 

304. 

95. 

- 

« 

- 

- 

1 

1230 

251. 

74. 

256.0 

23.8 

27.41 

0.803 

3 

1230 

306. 

95. 

215.0 

98.8 

- 

- 

5 

1230 

349. 

126. 

- 

- 

- 

- 

1 

1232 

249. 

73. 

283.1 

30.8 

43.72 

0.755 

5 

1232 

350. 

127. 

235.1 

72.5 

- 

- 

1 

1234 

249. 

71. 

278.3 

36.6 

41.90 

0.729 

2 

1234 

315. 

92. 

247.7 

47.9 

- 

- 

5 

1234 

349. 

127. 

168.9 

15.0 

15.84 

1.060 

1 

1236 

249 

70. 

274.1 

37.3 

38.72 

0.675 

2 

1236 

317. 

92. 

242.8 

53.3 

20.54 

0.502 

1 

1238 

250. 

69. 

265.2 

35.6 

39.77 

0.675 

2 

1238 

319. 

91. 

242.1 

59.6 

17.93 

0.619 

5 

1238 

351. 

130. 

214.7 

40.9 

36.72 

1.329 

1 

1240 

250. 

68. 

259.9 

34.6 

37.31 

0.635 

2 

1240 

320. 

90. 

243.2 

61.7 

18.14 

0.556 

1 

1242 

250. 

67. 

256.7 

34.0 

35.21 

0.601 

2 

1242 

322. 

90. 

242.5 

64.9 

18.45 

0.546 

1 

1245 

250 

64. 

254.1 

35.8 

33.48 

0.618 

6 

1245 

311. 

83. 

- 

- 

- 

- 

1 

1247 

251. 

63. 

250.5 

36.4 

33.63 

0.593 

5 

1247 

310. 

81. 

346.1 

74.1 

- 

- 

1 

1249 

252. 

61. 

246.8 

37.8 

33.45 

0.599 

6 

1249 

313. 

80. 

268.5 

62.2 

69.93 

0.944 

7 

1249 

353. 

132. 

- 

- 

- 

- 

1 

1251 

255. 

60. 

240.6 

39.0 

40.84 

0.576 

6 

1251 

312. 

80. 

271.2 

39.4 

62.41 

1.099 

1 

1253 

256. 

59. 

236.1 

39.9 

39.81 

0.555 

6 

1253 

314. 

79. 

260.7 

43.3 

58.90 

1.076 

1 

1255 

259. 

57. 

231.4 

41.9 

41.75 

0.624 

6 

1255 

315. 

80. 

248.2 

41.6 

57.34 

1.015 

1 

1257 

260. 

54. 

228.9 

44.4 

40.65 

0.677 

6 

1259 

319. 

80. 

237.8 

54.6 

57.82 

1.402 

1 

1259 

262. 

52. 

227.1 

46.9 

39.60 

0.685 

6 

1259 

320. 

79. 

236.7 

59.5 

54.73 

1.313 

1 

1300 

261. 

51. 

226.7 

48.3 

42.70 

0.680 

6 

1300 

322. 

80. 

234.2 

65.0 

57.53 

1.359 

1 

1302 

267. 

51. 

224.3 

50.2 

50.80 

0.701 

6 

1302 

324. 

81. 

231.5 

69.2 

55.49 

1.292 

1 

1304 

269. 

50. 

222.3 

51.8 

49.67 

0.685 

6 

1304 

326. 

82. 

229.2 

72.6 

53.43 

1.235 

1 

1306 

270. 

49. 

221.0 

52.8 

48.48 

0.670 

6 

1306 

327. 

81. 

229.0 

73.3 

52.15 

1.190 

1 

1308 

267 

47. 

221.1 

52.9 

51.53 

0.710 

6 

1308 

328. 

81. 

229.0 

72.6 

50.14 

1.149 

1 

1310 

270. 

46* 

220.9 

53.2 

50.72 

0.699 

6 

1310 

328. 

80. 

229.8 

69.8 

48.96 

1.143 

49 


A 


S ' 


Figure  60.  Remote  radar  display 
PPI,  200  km  range,  1225Z, 
November  2,  1974. 


Figure  61.  Remote  radar  display 
PPI,  200  km  range,  1230Z, 
November  2,  1974. 


Figure  63.  Remote  radar  display 
with  computer  generated  warning 
areas,  echo  contours  and  Victor 
airways,  November  2,  1974.  Pre- 
diction period  is  1240  to  1340Z. 
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Figure  02.  Remote  radar  display 
PPI,  200  km  range,  1240Z, 
November  2,  1974. 
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Figure  64.  Remote  radar  display 
PPI,  200  km  range,  1251Z, 
November  2,  1974. 


Figure  65.  Remote  radar  display 
with  computer  generated  warning 
areas,  echo  contours  and  Victor 
airways,  November  2,  1974.  Pre- 
diction period  is  1251  to  1351Z. 
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Figure  66.  Remote  radar  display 
PPI,  200  km  range,  1300Z, 
November  2,  1974. 


Figure  67.  Remote  radar  display 
with  computer  generated  warning 
areas,  echo  contours  and  Victor 
airways,  November  2,  1974.  Pre- 
diction period  is  1300  to  1400Z. 


r* 


51 


Figure  68.  Remote  radar  display 
PPI,  200  km  range,  1310Z, 
November  2,  1974. 


Figure  69.  Remote  radar  display 
with  computer  generated  warning 
areas,  echo  contours  and  Victor 
airways,  November  2,  1974.  Pre- 
diction period  is  1310  to  1410Z 


7.  SUMMARY,  CONCLUSIONS  AND  RECOMMENDATIONS 
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The  conclusions  and  recommendations  Included  here  are  divided  Into  three 
areas:  A)  Phase  II  objectives,  B)  software  refinements,  and  C)  future  tests 
and  suggestions  for  Implementation. 

A.  Phase  II  Objectives 

After  consideration  of  research  on  storm  motion  and  predictability  In 
progress  at  NSSL,  Lincoln  Laboratories,  Massachusetts  Institute  of  Technology, 
and  the  National  Weather  Service  Techniques  Development  Laboratory,  we  con- 
clude that  statistical  techniques  for  automation  of  tracking  and  warning  pro- 
cedures are  Incomplete  and  probably  will  require  substantial  improvements  in 
hardware  (e.g.,  more  sophisticated  radars  and  signal  processing  equipment) 
and  significant  changes  in  the  operational  configurations  of  manpower.  We 
believe  that  until  such  time  as  our  ability  to  measure  and  understand  severe 
storm  dynamics,  a simplified  man-machine  mix  using  echo  centroids  for  track- 
ing and  extrapolating  echo  motion  represents  the  best  technique  for  issuing 
advisories  at  the  Flight  Service  Station. 

We  also  believe  that  the  remote  terminal  described  in  Phase  I is  an 
adequate  method  for  remotely  displaying  radar  Imagery;  and  the  Phase  II  study 
has  demonstrated  the  feasibility  of  using  the  R,  9,  coordinate  system  as  a 
graphics  terminal.  To  this  end  we  have  developed  logic  which: 

1.  Shows  current  echo  positions  and  coverage  by  displaying  contoured 
echoes  at  user  defined  area  and  intensity  threshold  and  shows  that 
echoes  from  severe  storms  can  be  Isolated  routinely  using  an  inten- 
sity threshold  of  30-40  dBZ  and  an  area  threshold  of  150-1000  km^. 
This  threshold  Is  sensible  because  ninety-nine  percent  of  the’  storms 
with  hail  have  a maximum  intensity  above  30  dBZ  and  moderate  turbu- 
lence can  be  expected  in  storms  of  that  intensity. 

2.  Displays  a graphic  warning  area  which  is  derived  from  echo  size  and 
motion  and  expanded  downstream  for  a user  specified  prediction  inter- 
val to  show  a measure  of  the  storm's  predictability. 

3.  Provides  the  user  a choice  of  two  computer-generated  background 
reference  maps  — - the  State  of  Oklahoma  outline,  suitable  for  a 
400  km  range  display,  and  the  low  level  Victor  airways,  suitable 
for  a 200  km  range  display.  Other  maps  can  be  added  with  only  a 
slight  Impact  on  the  existing  logic. 

4.  Achieves  high  visibility  and  easy  interpretation  by  using  different 
grey  shades  for  the  graphics  elements.  The  best  results  were 
achieved  when  the  warning  areas  were  coded  "bright",  the  echo  con- 
tours coded  "medium",  and  the  reference  maps  coded  "dim". 

What  we  have  not  done  is  to: 

1.  Generate  alphameric  messages  on  the  remote  radar  display  system. 
Although  possible,  this  was  rejected  for  three  reasons: 
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a)  degradation  of  the  display's  resolution  as  a function  of  range 
would  require  excessively  large  characters  to  maintain  readability; 

b)  the  logic  would  be  complex,  time  consuming  and  require  additional 
core,  and  c)  low  cost,  high  speed  alphameric  displays  are  available 
for  use  as  a satellite  display. 

2.  Transmit  simulated  real  time  test  cases  of  the  graphics  products  to 
the  Wiley  Post  FSS.  However,  their  personnel  were  kept  Informed  of 
our  work  and  on  occasion  did  have  an  opportunity  to  see  a simulated 
test  case  at  NSSL.  Reactions  were  quite  favorable  to  the  graphics 
described  above. 

3.  Incorporate  growth  and  decay  trends  within  the  graphics.  They  were 
found  to  be,  on  a storm-by-storm  basis,  too  short-lived  to  extrapo- 
late. General  statistics,  such  as  echo  coverage,  average  Intensity, 
etc.,  were  not  acceptable  to  FSS  briefers  as  guidance  Information. 

B.  Software  Refinements 

The  computer  techniques  used  for  this  report  have  proved  reliable  and 
stable  and  no  major  modifications  are  suggested.  Small  changes  will  be  nec- 
essary to  facilitate  real  time  operation  and  improve  aesthetic  appearances. 
Although  developed  for  an  R,  9 coordinate  system,  the  logic  can  be  adapted 
with  a considerable  modification  effort  to  an  X,Y  type  display.  Additional 
modifications  of  existing  logic  might  include: 


1.  Creating  a permanent  file  for  the  echo  shape  function  so  that  a 
contoured  echo,  regardless  of  the  time  of  the  last  entry,  will 
retain  Its  shape. 

2.  Converting  the  logic  to  determine  the  airports  in  the  paths  of 
echoes  from  a cone  to  the  rectangular  area  used  for  generating 
warning  areas.  (It  may  be  desirable  to  perform  the  airport  search 
within  the  program  area  which  generates  graphics.) 

3.  Routing  radar  PPI  information  Into  the  computer  for  blending  with 
background  reference  fields  before  displaying  on  the  remote  terminal. 
This  would  require  a hardware  change  and  should  be  switch  control- 
lable to  avoid  computer  dependency. 

4.  Including  the  echo's  depth  when  drawing  a warning  box.  The  contour 
of  large,  slow  moving  echoes  may  actually  circumscribe  a box  under 
the  existing  logic. 

5.  Introducing  the  weighting  function  into  the  calculations  of  and 
0(j.  Currently,  only  the  centroid  positions  are  weighted  to  have 
decreasing  influence  with  time.  However,  the  errors  and 
should  also  have  decreasing  influence  with  time.  We  know  that  with 
few  points  to  estimate  echo  motion,  the  errors  will  be  large.  How- 
ever, as  more  points  are  added  we  should  not  be  penalized  by  the 
large  earlier  errors.  Conversely,  if  having  tracked  a storm  for 
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awhile,  its  predictability  decreases,  we  should  be  made  aware  of 
this  fact. 

6.  Dividing  the  scope  into  sectors.  Some  program  speed-up  could 
probably  be  realized  if  certain  quadrants  did  not  have  to  be 
scanned  for  echo. 

C.  Future  Tests  and  Suggestions  for  Implementation 

Based  on  experience  with  the  three  case  studies,  the  author  believes 
some  skill  is  needed  to  use  the  echo  centroid  tracking  logic  presented  in 
this  report.  Basic  to  a successful  operation  is  the  selection  of  area  and 
intensity  thresholds  and  sampling  Intervals  for  the  different  types  of  storm 
situations  (e.g.,  isolated  severe  storms,  squall  lines,  and  slow  moving 
flood  producers).  We  believe  that  some  training  of  FSS  personnel  will  be 
needed . 

A second  consideration  is  the  configuration  between  man  and  machine. 

One  possibility  is  to  have  one  person  responsible  for  identifying  the  haz- 
ardous storms  and  operating  the  echo  tracking  logic.  Then,  all  user  (pilot 
briefer)  requests  are  channeled  to  the  "hazards"  briefer  who  is  cognizant 
of  the  complete  echo  pattern.  A second  possibility  is  to  allow  each  pilot 
briefer  access  to  a computer  terminal  which  lists  input  and  output  of  echo 
locations  and  velocities  (e.g.,  echoes  which  have  been  assigned  user  num- 
bers for  tracking  could  be  so  labeled) . 

One  very  Important  question  which  needs  to  be  answered  is  how  much  time 
is  needed  for  decision  making.  Not  only  Initially  when  establishing  thresh- 
olds, but  during  operation  when  matching  echoes  manually.  This,  of  course, 
will  vary  with  radar  patterns  of  differing  complexity.  (The  author  had  the 
advantage  of  being  able  to  study  the  numerical  results  at  leisure,  but  was 
handicapped  by  not  having  a simultaneous  radar  display  for  comparison.) 

For  the  above  reasons,  a limited  operational  test  is  needed.  This  test 
should  be  divided  into  two  sections  of  one  month's  duration  each  during  a 
storm  season.  The  first  half  of  the  test  would  be  to  prepare  software  logic 
for  operational  use  and  to  test  (a)  the  feasibility  of  operationally  using 
the  existing  logic,  and  (b)  determining  the  optimum  man-machine  configuration. 
It  should  Include  sufficient  hardware  (e.g.,  remote  alphameric  terminal  con- 
nected to  a time  share  larger  computer)  to  allow  efficient  operation  of  the 
ECHOPRED  logic  without  major  modifications. 

If  the  outcome  of  the  first  test  is  favorable,  the  system  should  then  be 
Installed  at  a Flight  Service  Station,  such  as  Wiley  Post,  for  operational 
test  and  evaluation  by  FSS  staff. 
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APPENDIX  A 


Annotated  ConDand  Structure,  June  6,  1975 

C3MMAN3 

TWC 

HHMN 

2400 

6C0 

20 

C3MMAMQ 

ACC 

A^EA/nrE^SITY 
4 250  4 

DAY/  THE/  TILT 
157  2040  0 


2 VST  ECHOES  FOUND  WITH  AREA  GREATER  THAN  250  SQ  KN 
AREA/A2IHUTH/RANGE 

986  29  i36C.eti  1 ~ wcM»  iioUZ  tfuf  tAAcfcotg  Only  \JtAjy  i^tAJong 

1907  335  157  2 - new  ( (l/ST)  ceZti  thti  Aun 


2 SEV  ECHOES  FOUND  WITH  AREA  GREATER  THAN  250  SG  KN 

area/azimuth/range 

501  327  161  ^MAe.  ^&v^A&  czZti  ^bedded,  in.  tne 

359  340  155  ceJiti  dbo\ft 

DAY/  time/  tilt 
157  2045  0 

STOP 


2 VST  ECHOES  FOUND  WITH  AREA  GREATER  THAN  250  SC  KM 
AREA/AZIMUTH/RANGE 

1127  32  137cetfc  1 - Aome  gAjouoth  evident 

1974  335  153eett  2 


2 SEV  ECHOES  FOUND  WITH  AREA  GREATER  THAN  250  SC  KM 


AREA/AZIMUTH/RANGE 
407  328  154 

496  338  149 


57 


J 


i 


C3N^A*I3 

E^T 


N/ 

HHMi/  AZM/  RN6/ 

DIRECTION/ 

SPEED/ 

STDOS/ 

STOTM 

1 

20A0 

29.  136. 

1 

20A5 

32,  137. 

292.5 

92.9 

0.  00 

0.000 

2 

20A0 

335.  157. 

2 

20A5 

335.  153. 

335.0 

51.5 

0.00 

. 0.000 

0 

0 

0.  0. 

CONNANO 


BTIHE/ETIME/RANGE/OVERLAY/ECHO  NUMBERS 

21B0  2200  200  VIC  1 2 0 0 • 0 0 • • 

COMMAND 

ACC 

AREA/ INTENSITY 
250  A 

DAY/  TIME/  TILT 
1ST  2050  0 

STOP 


A VST  ECHOES  FOUND  MITH  AREA  GREATER  THAN  250  SC  KM 
AREA/AZIMUTH/RANGE 


1180 

3A 

lA0cee£  / 

353 

61 

188  ctU  3 - new 

30A 

310 

163 ^ 

188A 

336 

IA7  cM  1 

1 SEV 

ECHOES 

FOUND  WITH  AREA 

AREA/AZIMUTH/RANGE 

9A3 

335 

lAA 

(wcdZ  6e  KejtcUmd  In  next  PPJ 


COMMAND 

ENT 


N/ 

1 

HHMM/  AZM/  RN6/  DIRECTION/ 
2050  3A.  1A0. 

SPEED/ 

STDCS/ 

STOTM 

2 

283.1 

2050  336.  1A8. 

78.7 

IA.36 

0.182 

0 

319.0 

0 0.  0. 

58.2 

17.90 

0.098 

' \:t 


f-' 


( UL. 

I • ’ 


» r-fi*, 
f 


I : 

»• 


58 


COMMAND 

DIS 

.BTIME/E TIME /RANGE /QVERLAY/ECHO  NUMBERS 
21B5  2205  zee  VIC  1 2 0 0 

COMMAND 

ACC 

AREA/INTENSITY 
250  A 

DAY/  TIME/  TILT 
157  2055  0 

STOP 


3 VST 

ECHOES 

FOUND  WITH 

AREA/AZIMUTH/RANGE 

1225 

36 

145  zeJUi  1 

388 

61 

188  CZtt  3 

2509 

324 

147 

1 SEV 

ECHOES 

FOUND  WITH 

AREA/AZIMUTH/RANGE 

895 

336 

140 

(HO  oAA-cgrvnenC  AzoeAxU.  imaJU,  cjeUU  in  tinz 
{meAgzd,  may  bz  A.z&utt  mdoA  powzn  ituctuation 


COMMAND 

ENT 


N/ 

HHMM/ 

AZM/  RNG/  DIRECTION/ 

SPEED/ 

STDCS/ 

STDTM 

1 

2055 

36. 

144. 

276.9 

75.9 

16.73 

0.162 

3 

2050 

61. 

188. 

3 

2055 

61. 

188. 

360.0 

0.0 

0.00 

0.000 

0 

0 

0. 

0. 

COMMAND 

DIS 


BTIME/ETIME /RANGE /OVERLAY/ECHO  NUMBERS 
2110  2210  210  VIC  1 3 2 0 


59 


COMMAND 

ACC 

AREA/INTENSITY 
25C  A 

DAY/  TIME/  TILT 
157  2100  0 

STOP 


5 VST  ECHOES  FOUND  WITH  AREA  GREATER  THAN  250  SC  KM 


I 


AREA/AZIMUTH/RANGE 


1077 

38 

I46ce^e  1 

340 

61 

189cee£  3 

357 

298 

194cete  4 

561 

310 

155ce^e  6 

1656 

335 

136catt  2 - cent/LO^id  po^^OUon  comZitaU  loUh  hut  ceJLL  2 

position 


1 SEV  ECHOES  FOUND  WITH  AREA  GREATER  THAN  250  SC  KM 

AREA/AZIMUTH/RANGE 
95A  335  133 

COMMAND 

ENI 


N/  HHMi/  AZF/  RNG/ 

DIRECTION/ 

SPEED/ 

STDCS/ 

STOTM 

1 2100  38. 

146. 

277.1 

73.1 

15.58 

0.159 

2 2100  335. 

136. 

334.2 

64.3 

23.26 

0.102 

-3  2055  61. 

188. 

3 2100  61. 

189. 

241.0 

6.1 

0.00 

0.000 

0 0 0. 

0. 

COMMAND 

OIS 

BTIME/ETIME/RANGE /OVERLAY/ECHO 

NUMBERS 

2115  2215  200  VIC 

1 2 

3 0 

0 

0 0 

60 


I 


B 

r' 


• s 


COMMAND 

ACC 

AREA/INTENSI1Y 

2se  A 

DAY/  TIME/  TILT 
137  2105  § 

STOP 


6 VST  ECHOES  FOUND  tolTH  AREA  GREATER  THAN  250  SC  KM 


AREA/AZIMUTH/RANGE 

762  33  1561  Ic&tt  1 ipt/itp  MiLt  to  ACC  6pLLt  A.ejncuju  b&io^e.  itiakinQ 

252  A5  122  f \nm  asugment  ^ 

A32  60  187CCCZ  3 

600  2 97  19 2ccZZ  4 - gAouith  but  ceMtAotd  po&itLon.  con6t6tent 

586  308  152CCU  6 

19A6  3 36  132cctC  2 - gfvovoth  cv-ufcjtt,  but  czntfioid  po&ltJLon  ti 

conAlttu/it 

1 SEV  ECHOES  FOUND  WITH  AREA  GREATER  THAN  250  SC  KM 

AREA/AZIMUTH/RANGE 
1087  333  127 

COMMAND 

ENT 


N/ 

’ HHMM/  AZM/  RNG/ 

DIRECTION/ 

SPEED/ 

STOCS/ 

STDTM 

2 

2105 

336. 

132. 

331.8 

62.8 

2A.05 

0.1A2 

3 

2105 

60. 

187. 

137.5 

11.7 

22.61 

2.883 

A 

2100 

298. 

19A. 

6 

2100 

310. 

155. 

A 

2105 

297. 

192. 

3 56.8 

A9.0 

0.00 

0.000 

6 

2105 

308. 

152. 

9.8 

76.7 

0.00 

0.000 

0 

0 

0. 

0. 

COMMAND 

DIS 


BTIME/ETIME /RANGE /OVERLAY/ECHO  NUMBERS 
2120  2220  200  VIC  2 3 A 6 


61 


COMMAND 

ACC 


AREA/INTENS11Y 
250  4 

DAY/  TIME/  TILT 
157  2LI0  0 
STOP 


4 VST  ECHOES  FOUND  WITH  AREA  GREATER  IHAN  250  SC  KM 
AREA/AZIMUTH/RANGE 

879  40  162  no  (L&6<gninejit,  htXXX.  not  convtnc^xi  tcho  1 t/oon't  Aeappe/Vi 

531  59  107  cete  3 

772  296  188cett4 

1854  332  124  cate  2 


I SEV  ECHOES  FOUND  WITH  AREA  GREATER  THAN  250  SC  KM 

area/azimuth/range 

726  337  123 

COMMAND 

ENT 


N/ 

HHMM/  AZM/  RNG/  DIRECTION/ 

SPEED/ 

ST  DCS/ 

STDT 

3 

2110 

59. 

167. 

138.9 

19.7 

20.00 

2.746 

4 

2110 

296. 

188. 

345.0 

55.0 

12.69 

0.110 

2 

2110 

333. 

124. 

339.1 

66.1 

40.99 

0.230 

0 

0 

0. 

0. 

COMMAND 

OIS 


btime/etime/range/overlav/echo  numbers 

2125  2225  200  VIC  123400000  j 

COMMAND  ’ I 

ACC  i 


AREA/INTENSITY 
250  4 


L 

j 

i 

! 


I 


4 VST  ECHOES  FOUND  WITH  AREA  GREATER  THAN  250  SC  KM 

AREA/AZ IHUTH/RANGE  t oA.ea  -cA  coifU'Utejfit  ulcth  tlUt  dzZJi  1 cJita  and  -6s 

1060  AA  i62ce^  lionZy  one  -in  the  a/iea,  naptd  movement  -6s  tndlcAted 

A97  61  lascete  3 ' 

B15  295  l8Acet6  4 

1868  333  Il8ce66  2 {cett  6 and  maZteA  paM.  celt  1 iptit 

dfLopped] 

1 SEV  ECHOES  FOUND  WITH  AREA  GREATER  THAN  250  SC  KM 

AREA/AZIMUTH/RANGE 
797  337  117 

COMMAND 

ENT 

N/  HHMM/  AZf'/  RNG/  DIRECTION/  SPEED/  STDCS/  STDTM 


1 

2115 

AA> 

162. 

272.1 

79.2 

15.85 

0.171 

3 

2115 

61. 

188. 

135.5 

7.1 

17.95 

3.333 

A 

2115 

295. 

18A. 

3A0.6 

57.A 

12.  3A 

0.109 

2 

2115 

333. 

118. 

3A1  .2 

68.1 

38.06 

0.215 

0 

0 

0. 

0. 

command 
DIS 

BTIME/ETIME/RANGE/OVERLAY/ECHO  NUMBERS  - 

2130  2230  200  VIC  13A2000B0  ] 

COMMAND 

ACC  j 

i 

AREA/INIENSIIY  | 

250  A \ 

DAY/  TIME/  TILT 
157  2120  0 

STOP 


I. 

63  ; 


5 VST  ECHOES  FOUND  WITH  AREA  GREATER  THAN  250  SC  KH 


AREA/AZIMUTH/RANGE 

1236 

48 

160  cete  I , 

452 

60 

189  celt  3 

674 

295 

178  celt  4 

323 

305 

131  ceU  5 - 

169S 

334 

115  celt  2 - 

1 SEV 

ECHOES 

FOUND  W1 TH  1 

AREA/AZIMUTH/RANGE 

608 

334 

107 

poshly  old  cdJU  6,  but  jait  OA  tou,y  to  fieaiM^n 
&onz.  cuinja.  dccAcoAC,  but  cojntAold  poi-ctlon 
conitstent 


COMMAND 

ENT 

N/  HHMM/  AZM/  RNG/  DIRECTION/  SPEED/  STOCS/  STDTM 


1 

2120 

48. 

160. 

277.9 

80.7 

44.36 

0.159 

3 

2120 

60. 

189. 

155.2 

6.2 

18.65 

3.435 

4 

2120 

295. 

178. 

329.4 

57.6 

25.53 

0.112 

2 

2120 

334. 

115. 

340.5 

66.6 

37.05 

0.274 

0 

0 

0. 

0. 

COMMAND 

CIS 

BTIME/E TIME/RANGE /OVERLAY/ECHO  NUMBERS 

21TS  2235  zee  VIC  134200000 


COMMAND 

ACC 

AREA/INTENSITV 
250  4 


DAY/  TIME/  TILT 
157  2125  0 
STOP 


64 


6 VST  ECHOES  FOUND  WITH  AREA  GREATER  THAN  250  SC  KH 


AREA/AZIMUTH/RANGE 

I89ce££  7 - nm 

I65ce££  t - gAjowtk  t/Lond  JjndJjMtsd  cxjn&^tojnt 

l9lceZt  3 - ItVLge  oA&a  dtcAmAz,  bat  cmtAoid  position/' 
Ibbczti  4\  iftapid  gfiowth  ob6e/ived  In  both  cett6,  but 
I26c&tt  5/  (cent'totd  poitttom  cjon&t&tz^ 
l05cetZ  2 


377 

13 

1421 

50 

267 

60 

906 

2 96 

668 

305 

1639 

332 

I SEV  ECHOES  FOUND  WITH  AREA  GREATER  THAN  250  SC  KM 

AREA/AZIHUTH/RANGE 
637  336  10A 

COMMAND 

ENT 


N/ 

HHMM/  AZM/  RNG/ 

DIRECTION/ 

SPEED/ 

STOCS/ 

STOTM 

5 

2 1 20 

305. 

13I. 

• 

I 

2125 

50. 

165. 

2 79.9 

GO 

• 

41.50 

0.155 

3 

2125 

60. 

191. 

184.5 

6.5 

CD 

• 

3.190 

4 

2125 

296. 

166. 

312.9 

66.7 

50.69 

0.313 

5 

2I25 

305. 

126. 

305.0 

61.0 

0.00 

0.000 

2 

2125 

332. 

105. 

341.5 

69.  1 

36.33 

0.3S9 

0 

0 

0. 

0. 

COMMAND 

CIS 

BT IME/E TIME /RANGE /OVERLAY/ECHO  NUMBERS 
2140  2240  200  VIC  I 3 4 5 

COMMAND 

ACC 

AREA/INTENSITY 
250  4 

DAY/  TIME/  fiLT 
157  2130  0 

STOP 


65 


r 


1 


1 

i 


1 

' I 


1 

i 


[i 


5 VST  ECHOES  FOUND  WITH  AREA  GREATER  THAN  250  SC  KM 


ARcA/AZIMU TH/RANGE 

515 

12 

l87ceZX  7 - fiapld  gfumth  inam  pA.ev^uA  PPI 

1376 

50 

173ce££  1 - gfuMth  txmi  KnveJU^lng 

935 

295 

163  cete  4 

6AA 

3 0A 

122  cJiJUL  5 

1732 

3 32 

101  C.M  2 ic&tt  3 cUwpped  by  p^ogfum) 

2 SEV 

ECHOES 

FOUND  WITH  AREA  GREATER  THAN  250  SC  KM 

ARcA/AZIMUTH/RANGE 

260 

57 

156 

62A 

333 

93 

COMMAND 

ENT 


N/ 

HHMM/  AZM/  RNG/ 

DIRECTION/ 

SPEED/ 

STDCS/ 

STDTM 

7 

2125 

13. 

189. 

7 

2130 

12. 

187. 

71.1 

A6.  I 

0.00 

0.000 

1 

2130 

50. 

173. 

278.3 

81.8 

46.15 

0.167 

A 

2130 

295, 

163. 

310.0 

68.5 

47.67 

0.31fr 

5 

2130 

30A. 

122. 

318.3 

55.8 

14.75 

0.121 

2 

2130 

332. 

101. 

3A1.7 

69.5 

34.56 

0.383 

0 

0 

0. 

0. 

COMMAND 

OIS 

BTIME/E TIME /RANGE /OVERLAV/ECHO  NUMBERS 

21A5  22A5  A20  STA  7 1 A 5 2 0 0 

command 

bVc 


66 


i 


1 

t 

I 


1 


0 


f- 


I 


APPENDIX  B 


Annotated  Command  Structure,  June  16,  1975 


ClJMMA.'Ji) 

1 rtC 

^iCD 

2H 

4”) 

ACC 

AK  f A / I j 1 r S I r Y 
252  ^ 

DAY/  TIML/  TILT 

167  183C1  / 


1 VST  LCHJf'i  FQU'JO  MTH  '*RcA  JrlLATFR  THAN  2‘>kl  SC  KM 

■\kC  A/AZMLTH/RANGF. 

17  75  2FZ  ii]QC.M  1 - new  {wWL  tnack  VST  ceZti  tftu  mm 


1 SrV  CCHCCS  FJUMO  wITH  AREA  GREATER  THA\  2^^  SC  KM 
iurh/RANGL 

1^32  2^1  IH9  Ignore  SEl/  ceZU  -imbedded  -in  the  VST  eeUA  above 

n\Y/  TlMf-  / ri  LT 
167  1R35  2 

STOP 


1 VST  ^CHOrS  FOUND  wlTh  \KcA  .iRtATEK  THAN  25;)  SC  KM 

A<tA/AZ  li-lUlH/R.\NoE 

222  1 275  lP6Cete  1 


I SLV  hCH.;ES  FOUND  klTH  AKLA  vJREAlER  THAN  250  SC  KM 

AREA/AZIMUTH/RANUE 
17AA  279  IP6 


67 


COMMA  ■40 
fc  4T 

4/  HHMM/  ROG/  DIRtCTIGN/  SPEL-D/  STCOS/  STDIM 

1 1G30  kbii,  188. 

I 1835  270,  186, 

33S,  U 45.9  a.KiJ 

a i)  «, 

CUM  '1 AOO 

0 1 s 

eriME/!'  I lM'l/RA'4Gt-/UVF.RL&Y/ECHU  i'iUMb£R'-> 

18  35  1 93  5 2Hi6  VIC  I H <i  13  yi  M y 

CU.MM  IMH 
ACC 

AR^  A/I  JTl  NSIT Y 
25<A  4 

OAY/  TIMF/  riLT 
16  7 18413  12 

STOP 


i VST  i.CHUES  POUOO  MTH  vR£.\  GPLAIGR  THAN  25ii  SC  K*^ 
AK£A/A/I,MUTH/MA;M.jE 

264  8 278  I R3  cetd  1 - growth  due  to  AquuxJU  line  enteJUng  PPI 


2 SFV  ECHUilS  FOUND  wlTri  ..RtA  UKlATLR  ThAi\  25M  :jC  KM 

AR3A/AZ IMUTh/kANGE 
1429  26S  188 

4 6iJ  29  3 185 

CU-M.MANI) 

E 4 1 

■4/  hH'3,3/  \ZM/  RNG/  OIRlUION/  SPEED/  STCLS/  STCi.m 

1 1H4D  2 78,  18  3. 

331,3  4 >,  3 6,59  (1,l37J 

2 •'  2 » i3 , 

. »'<  *A  4 > 

<1  > 

•II  * / II-  /•<  * 4Gf /.JVC  RL-V/EtrlO  NUMbERS 

• . I I'  , 4 / V I ' I ,1  13  2 3 W 


A^EA/IiTiNSirv 
25i;  A 


.J4Y/  TIME/  TILT 
lfc7  HAb  0 

srup 


l VST  JCT<iJCS  FOUVO  WITH  'u<£A  G.-EATl'<  25i?  SC 

>XtA/M  IHL  IF7'5>!^NGt 

306  7 2 7S  179  cM  1 - line,  itui  gfLouiing 


,,  1 SEV  tCHGES  found  MTH  UktlTEK  THAN  25k)  SC  KM 

I i 

.Wr  A/A/I.^iLrH/7),N05 
1661  2H  l?l 

Cu;-ri;,  kU' 

E TT 

J/  HH1M/  AZfl/  iirJG/  OICE:,  rlUN/  SPFtD/  STCDS/  STOn'^’ 

I IHA5  279.  170, 

)k’?.  S 37,2  31.11  0.516 

0 2 0 ..  0 . 

OuMM  ANi' 

01  S 

3 r I :U  / ■:  r I -1  - a ng  f / ov  ^ -a  a y /k  c h o . jum  k t s 

ldA5  17Ad  230  VIC  I 0 0 0 0 0 ,)  g 

CuMiNAN.) 

ACC 

WI  jTcNSITY 
250  A 

DAY/  TIMr/  tilt 
16  7 1‘)5k-  0 

STOP 


69 


1 


i 

I 


I 


2 VST  tCHGCS  FOLIVO  «ITh  GkcATEK  TH^N  ^5'^  SQ  KM 

A/^Z  IMUfH/RAfJGL- 

2 36  3 2 76  I 74  ceZl  1 {cM  &pJUt  on  nofvth  tnd  tint,  viiJUi 

26  5 jt:i  . 186  no  diiignrmnt  unit  to  4ee  mate  celt  -66  aeXained) 


1 Stv  cCHlLS  FQUNU  WITH  -RtC  Gi'.LATcH  VHA\'  25^  SC 

A.<l-A/AZ1  •IGTH/RANGE 
1636  268  179 

CuMMA.JC* 

L- xiT 

HH3M/  M\/  RMo/  OlRESlIuN/  SPE'.Li/  STCDS/  SICIM 
I 18S-1  276,  174, 

314.0  S2.  I 43.65  J',«4  i 

ti  .)  >1,  e. 


CJMMAxli) 

015 

BT  I M /u  T I /R  A'JG F / UVr  RL  A Y /ESHU 

135^  r»prf  2A2  VIC  I 0 0 0 0 H 0 K 

Cl. -1,3  \M  1 
ACC 

AriCA/l  (Ff-  VSIf  Y 
252  4 

0\y/  TI‘1_/  KLI 

16  7 I 35 3 ? 

STOP 


1 VST  ECHiTS  FOLNO  WITH  .,PEA  Gl<eATtP  THAN  25  3 sC  K'' 

A<l  A/A/  1 'UTH/RAvjGf. 

2t75  2 73  I licett  1 - ajiea.  decAzoi-ing,  6matt  ceZt  dAjopped 


I SlV  CCHH'^S  FGUVD 

WITH  -9  6 A 

GKcA  T6  f 

THA,\  250 

SC  KM 

A 3cA/A/ I :UTH/ range 
1362  266  172 

CUMMA VO 

1 < r 

V/  m'^^/  *.  2,1/  RNG/ 

DlktC 1 lUN/ 

SPEtO/ 

sreos/ 

STCl  M 

1 1855  273,  171, 

324,  0 

t>3.  3 

46,94 

0.845 

0 .3  K.  0. 

70 

j 


CUMMi,'*n 

DIS 


df IML/iTiM"/«AvlGr-/CVri<LAY/FCHU  :\IL'WbFRS 

1855  I J55  VIC  1 0 H « 0 ^ e 

CUMMA;>n 

ACC 


4RE4/I  'J  IF-MSirV 
25e  A 

DAY/  TIME/  TILT 
16  1’  1^16  0 V. 

STOP 


I VST  cCHCIS  FdbNL)  wITH  -'TEA  oRlAIF^  THA\  250  SC  KM 
AKFA/A/rUlH/RANGE 

2505  271  l65ce££  1 - oAtcL  itctt  dtcAejCi&^ng , contMld  position. 

coM^ient 

l SCV  6CH._ES  F0U>JU  '«ITH  A<tA  GRtATEk  THAN  250  SC  KM 


ARE  i/A2I.-lUTH/RANGr 
12 7A  266  167 

CuMMAV'* 

ERT 

Hh-U/  RNG/  OIRLCTIUM/  SPELO/  GUOS/  STCI  0 

1 1900  271.  165. 

326.6  71.5  63.76  0.812 

o 0 0 . 0 . 

DIS 

UT  IMt/FT  lKt;/R  6MGL/UV;RLAY/ECHU  NUMfERS 

1700  ?O0O  200  VIC  10000000 


C0M/-U-7!) 

ACC 

ARFA/nrhGSITY 
250  6 

DAY/  tim:/  tilt 

167  1^05  V 


2 VST  eCHCES  FOUND  WHH  ^RcA  GkEMEK  THAN  2i>V^  SC  KM 
AREA/A2ifiUTH/«ANGE 

2238  262  I72ce££  2 - new)  -in  czntAo-id  poi-ition  i-igni^-icAnt 

J6t{  287  Iblce^  3 - new  (w^/i  thh>  &ptit,  findtiim 


I SEV  ECHCtS  F0U\)1  wITH  AREA  SRLATEk  THAN  2S«  SC  KM 

ARE4/A2I.-1L7H/RANGC 
l<iA9  263  166 

DaX/  time/  TILT 
167  191i!  H 

STOP 


2 VST  ECHOES  FOUND  VnITH  AREA  SKCATER  THA.-i  252  SC  KM 
ArltA/AZiMUTH/RANGE 

?38K  261  I67cet£  2 - oAejCL  and  centJw-id  po&ltion  com-utent 

262  28b  lA2ce£^  3 - lobv-Loiu  ^ftam  PPI  Aome  ced£  04  above 

Aap^d  mo-tion  indlccuttd) 

1 SEV  lCH’.'ES  FOUND  wiTH  ARtA  GkEATcP  THAN  25H  SC  KM 

AKEA/AZ IMLTH/RANSe 
1626  261  163 

COMM  \N"I 
E IT 


X 

X 

DIRlCTIUN/ 

SPEED/ 

STCDS/ 

srci  H 

2 1925  262.  170. 
2 1910  261.  167. 

305,9 

55.2 

•3.00 

J 1 J05  297.  151. 
3 in.)  2 85.  I'i2. 

31  5.6 

1 35.  9 

0.  00 

ECHO  SPEt-D  GlfiATtR  THAN  120  KM/HK,  CHECK 

THIS  U 

W • 

a 2 2.  0. 

COMMAND 

OIS 

3T IML/ETI M-  /M  ANGC/CV 

CkLAY/ECHC 

NUMbERS 

ivie  2010  2ae  vie 

2 3 

0 0 

0 

k) 

cummvho 

ACC 


I 


\ReA/I  JIENSITV 
252  A 

OAy/  TIME/  TILT 
16  7 l‘)15  2 

srup 


2 VST  ECHOES  FOUND  WITH  AREA  GREATER  THAN  SC  KM 

AREA/AZ IMLTF/RANGE 
2A31  261  159  dtU.  2 

36A  2E7  133  c.eXZ  3 - oAca  toht  >ut  pAev^oaA  PPI 


I SGV  tCHCFS  FOUND  WITH  AREA  GREATEP  THAN  252  SC 

AREA/AZ  I '■lUTE/RANGE 
I'tSZ  261  156 


COMMAND 

ENT 

N/  HHlri/  AZM/  RNG/ 

DlRECTIuN/ 

SPEED/ 

STCOS/ 

2 1915  261,  159. 

27o.  1 

69.  7 

39,15 

3 1915  287.,  133, 

287. Id 

IID.  A 

59,  72 

2 (2  C • 13 « 

COMMAND 

DIS 


STCfM 

3.13A 

2,29? 


Bl  IME/C  ri  vE/RANGE/OVERLAV/ECHU  NUr.3ERS 

1915  2'315  2452  VIC  232v)fe>2«2 

command 

ACC 


ARfcA/I  UENSirV 
2512  A 

DAY/  TIME/  TILT 
16  7 1 92  2 « 

STUP 


j 


73 


1 VST  LCHOlS  found  KITH  -\R£A  UKcATER  THAN  25i1  SC  KM 


•VKEA/A/.IMUTH/RANGE 

2792  26A  l<tbceti  1 - (cfctfci  2 and  3 meAgzd,  zxt^potated  position 

fiOfL  old  ce££  I appzwu  ^&cL6ona.ble) 

2 StV  ECHOES  FOUND  KITH  ARcA  GREATER  THAN  250  SC  KM 


AKrA/A/ HUTH/RANGE 
5 7a  2 5A  169 

591  265  136 


COMMAND 

ent 

N/  HH.1M/  12M/  RNG/  UH 
1 192,3  26A,  1A6, 


TIO'i/  SPEED/  STCDS/  STCIH 


COMM  AN!) 
DIS 


Bf  I WE /v  TlMt./^  ANGC /OVERLAY /ECHO  NUMBER 
1922  2222  202  VIC  I 0 0 


COMMAND 

ACC 


AREA/1  TILNSITY 
252  A 


DaV/  TIMc/  tilt 

167  1 125  C 

STOP 


1 VST  ECHOLS  FOUND  KITH  ftKEA  GREATER  THAN  250  SC  KM 


aRFA/A/HUTH/RANOF 
2757  26?  lAacete  I 


2 SEV  ECHCtS  FOUND  KITH  AREA  GREmTER  THAN  250 


AREA/AZ IKUTH/RANGE 
510  25?  165 

A35  2 72  126 


-4 


CUMMANO 

EMT 

■4/  HH4M/  (UHf  RNG/  DlRECFIJ'S/  SPE  tU/  STCOS/  STUFK 
1 1'425  263.  IA0, 

313./  73.4  41.  til  0.72k) 

\&  J iE).  0. 

CGMMA>4'J 

DIS 

BT  Ir«L/cTI  ME/RAr4GE/UVcKLAY/6(.Ha  MUMRf^S 

1925  2025  200  VIC  1030000 

C-JMMANIO 

\CC 

AKfcA/rUL  4SIIY 
250  4 

OAY/  TIME/  IILT 
167  1 <33  0 

STOP 


I VST  fcCHLES  FOliND  >»lTH  -.RtA  GKEATE-'i  THAN  250  bC  KP 


I 


Ei 


ARcA/AZI  ILlTH/ RANGE 

2560  263  l32ce^  / - OACa  AtOAJUng  to  (feCACOAC 


2 Stv  ECHUFS  POUND  k*  I TH  ARLA  GRLATtk  THAN  25  » bC  KP 

AKEA/AZI  -EtTH/RANGfe 
555  252  155 

442  259  119 

CONmano 

E a 

.4/  HH1M/  M\M  R 4G/  OIRhCMUN/  3HEH)/  STCOS/  SfClM 
I 1933  26  3.  1 32. 

315.3  77.3  46.93  3.690 

0 0 0.  0. 

COMMAND 

DIS 

RTIME/CTIME/kANGE/UVEkLAY/ECHU  numbers 

1930  2030  200  VIC  1 0 0 0 0 0 0 0 


t 

Si., 


5' 


I 


I 

f Gi.. 

fv 


75 


C0MMA>40 

ACC 


AKEA/I>tTE>iSITY 
25i  4 

DAY/  TIME/  TUT 
167  1*35  0 

srup 


1 VST  cCMJtS  FOUND  WITH  /«EA  25/  SC  KM 

AKtA/A/1 4UTh/«AN3E 

225 A 2^C  IZtcijtt  1 - oAta  de.cAtiuing  fuxp-Ldty 

2 Stv  ECHOES  FOUND  ulTH  ^RtA  C'iEATfK  THAN  25>»  SC  KM 

AKt  A/A/ I lUri/AAN  jF 
522  251  ISa 

525  2^fc  lU 

CUMMuN  ) 

EiT 

4/  HMN.N/  \/l/  RnG/  jmtCTIli.N/  SPEIC/  srccs/  STCI  ' 

1 1435  26»>.  126. 

3l  3*i  7»i,  7 AA,d5  f1,66C 

V ■)  /•  P. 

CU AMA\  ) 

JlS 

^UMt/ETI  M--/J  \NGf /OVtKlAV/FuMU  NuMItPS 

14)5  2<4i^  2PC  VlClVVOldPut; 
ACC 

AKcA/l  nTcNSIT  Y 
250  A 

DAY/  TIME  / TUT 
16  7 I'HV  0 
STOP 


1 VST  ECHi'ES  FOUND  WITH  AREA  GKcATER  THAN  250  SC  KM 

AREA/A/IMUTH/PANGE 
2390  25E  IZBceU  1 


• 2 StV  cCHOtS  FOUND  V*ITH  ARhA  GkcATER  THAN  250  SO  KM 


AMEA/A2  IMUfh/RAUiiE 
A76  2’3l  I'te 

61A  2<S7  103 


CUMMANO 
E iJ 

N/  HH  1M/  iZH/  RNO/  UIRlSIIUN/  SPEED/  STCDS/ 
I 19A0  256.  120. 

311./  79.3  43,72 

0 .)  0.  0. 


STCFK 


I COMM  iN  ft 

')IS 

3 r I ME  / ■ T I ME  /«  ANGl  /OVF.  '<L  -.Y  /ECHO  NUMbt-iiS 

1940  2v)4l2  2J0  VIC  10  0 0 0 0 0 0 

j COMMAN  ) 

\CC 


ukEA/iMriNSirv 
2 30  4 


34 Y/  TIM.-/  TILT  I 

lt)7  1 145  0 I 

STOP  ( 


I VST  cCh..tS  FOUVD  .*ITh  ARtA  GKcATEk  THAN  250  SC  KM 


4KE  A/A/ MtlF/RANGt  3 

2246  2-15  I 11  ce££  J 1 


2 SFV  ECHCtS  FOUND  klTII  AREA  G>%tATE’<  TMA\  25/  SC  KM  h 

AREA  /AA  IM.UTH/RANf.L 
417  249  13H 

367  2tl  95 


^ CUMMAN.i 
* cNT 


N/ 

HH-'.M/  VZM/  R <G/ 

DIRtC  IIUN/ 

iPEEO/ 

STCDS/ 

STCIM 

1 

1945  259.  Ill, 

30  9.0 

79.  1 

50.48 

0.525 

0 

0 0 . 0 . 

77 


I 


-I 


CJMHAMl 

OIS 

riTIME/i.iTIML/RANGf  /OVE«L  AY/ECHO  ivUMBE*<S 
I'iAS  2<)A5  2ae  VIC  I i3  0 ; 

CUMMA.>tJ 

ACC 


AKEA/lMTtNSITY 

uee  3 

OAY/  Tl!3e/  n LT 

Its 7 I ziceZZ  1 bnjokz  -onto  htvzfiaZ.  c.oA.eA  ;tune; 

a tcLAgeA  cuiza  at  a toweA  tntzmtty  uxu  a&ed  a^teA  tkU 

ttm.) 

2 SIR  tChOLS  FOUND  WITH  AREA  GkEATE.R  FHAMUHO  SC  KM 

ARcA/AZIMUTE/HANGE 
*512  35A  2A9l 

AdAfi  2fc  3 1 1*31  data  takzn  at  Zong  Aonge.,  no  o6^tgnment  mcudo. 


1 VSr  6CHCES  FOUND  wl  TH  \RuA  GREATER  THAM(322  sc  KM 

ARcA/AZl  iUTH/RANGE 
3362  3A5  2A6 

DAY/  TIME/  TILT 
1;j7  1950  0 


1 STR  ECHCES  FOUND  wlTli  AREA  GKEATEk  THA\ie00  SC  KM 

akea/a/i.-iuth/range 

A7A5  26?  WZcztZ  4 - new 

DAY/  TIME/  TILT 
167  1955  e 

STOP 


1 STR  ECH^jES  found  wlTH  AREA  GRcA  TER  THAN1000  SC  KM 
ARCA/AZ IMUTF/RANGF 

A7A3  26A  10Rce££  4 - Indtcotexi  motion  unasuaZZt/  iZoui 


78 


CUMM^NH 

EO 


N/ 

HHMM/  \m/  RNG/  OIRuCTIuN/ 

SPEED/ 

STCCS/ 

STCfM 

4 

1963  263.  110. 

4 

196S  264.  108, 

219.9 

36.  2 

0.00 

0.00,5 

0 

Y , 0 . 

DIS 


BT  IME/rTIMC/'^ANGE/OVERLAY/eCHO  NUMbERS 

H5  5 26»5b  2'A9.  VIC  A 0 0 0 3 


CUYMAN'O 

ACC 


AREA/I  UENSIT  Y 
100  0 1 

uay/  time/  tilt 

167  ?00Z  0 

STOP 


I STR  ECHOES  FOUND  WITH  /^REA  GREATER  THAM00H  SC  KM 
AREA/AZ IMbTH/RANGE 

A8A'3  267  >77ce££  4 - [cznihoZd  poiX^Uan  change  veAy  ta/ige, 

echo  ieeiM  unitabte} 


command 
e JT 

> 

N/  HH  H/  iZM/  R JG/ 
4 2003  267,  97, 

D 1 R lC  II  (J''i / 

SPEED/ 

STCDS/ 

STCI 

0 0 0.  0. 

236.9 

71.  2 

28.38 

1.627 

COMMAND 

DIS 

BT  I ME/ET  I Mr /RANGE /UVERL  AY/ECHO  \UMt3ERS 

2000  2100  200  VIC  A 0 0 3 0 0 0 

COMMAND 

ACC 

area/I  jtensity 


1000 

3 

DAY/ 

TIME/ 

IILT 

167 

2005 

0 

2 STR  tCHGtS  FOUMD  WITH  AREA  GREATER  THAN14J0?'  SC  KM 


AREA/AZIMUTF/RANGE  c/uuzge 

ABA 5 263  c.iZt  5 - new)  {6iQYu.(tlcjoint  aJido.  and  crntnald  potAJUbonl 

1J68  29A  155  no  cni^ment  - tofign  patch  oi  6thjatlionm  JuUn 


DAY/  TIME/  TILT 
16  7 2!Uk3  Z 

STCP 


2 STR  ECHOlS  FUU  id  WITH  AREA  GREATER  THAMlfeiM'B  SC  KM 
4REA/A/IMUTH/RANGE 

263  7A  ceZl  5 - oA(ua  and  ccnt/vold  po6tXton  coM-utcnt 
1322  295  155  no  aiA-ignmznt  - tittle,  motion  p/icvlouA  PPJ 

COMMA Tb 
E7T 

9/  HrHM/  4ZM/  RNG/  DIRECTION/  SPEED/  STCUS/  STCTM 
5 2005  263.  82. 

5 2 010  26  3.  7A. 

263.0  Lv<5.0  0.20  0.000 

0 0 2.  0. 

COMMAND 

DIS 

ilTIME/ET  IVt/RANGE/OVERLAY/ECTC  NGMCERS 

2010  2110  200  VIC  50000002 

COMM  A D 
ACC 

AREA/I  TTENSITY 
1000  3 

DAY/  riMu/  TILT 
167  2015  0 

STOP 


2 STK  ECHOES  FOGND  WITH  ARcA  GREATER  THAN1000  SQ  KM 
AREA/AZIMLTH/RANGE 

A186  263  6A  cell  5 - Aapld  motlori.  Indicated 

1A3A  296  1 53  no  oAA^nmen^  - Aee  above 


80 


COMMAND 

fc'IT 

A/  AZM/  R,\iG/  DIRFCTIUM/  SPetO/  STCDS/  STCF'-I 

5 2315  263,  64. 

26  3.4)  IV3.  A ^,21  0,082 

kS  3 2.  0. 

command 

DIS 

rtriME/CTI.-!t/R  ANGE/OVCKL.iY/PCHU  JUMBFPS 

2015  2115  200  VIC  5 0 0 0 0 0 0 0 

COMMAND 


ACC 

AREA 

/I 

J T E N 

SIT  Y 

1000 

3 

DAY/ 

T 

IMf  / 

TI  LT 

16  7 

2 

02  0 

0 

2 SIR  tCHCEb  FOUND  WITH  aRUA  OklATER  THAN1000  SC  KM 
AREA/AZINUIF/RANGE 

516  2 2E5  =!1  aeZt  6 - neu{  no^k  e.nd  tinz  coiued 

15A3  2 46  1A5  wo  cu,6-iqm&fU  -cn  centfiold  and  oAea. 


I VST  ECHCes  FOUND  WITH  ARtA  GREATER  TFAN1000  SC  KM 

AREA/AZIMUTH/RANGE 
1250  2A0  R2 

OAY/  TIM'./  TILT 
167  2025  0 

STOP 


2 STR  ECHCES  FOUND  WITV<  ARlA  GREATER  THAN1000  SC  KM 
AREA/AZ IMUTh/RANGE 

5203  2PA  A6  ceZt  6 - oAta  and  ze.ntAold  po^Atcon  cxfm-utent 

IAA2  257  I AP  WO  OAA^gnmewi 


I VST  echoes  found  WITH  AREA  GREATER  TFAN1000  SC  KM 

AREA/AZ IMLTF/RANGE 
1247  23P  77 


81 


I 

1 


EiMT 

N/  HHMM/  AZM/  R.NG/  OIRtCTIUN/  bPEfcO/  STCGS/  STCT.l 
6 2020  2E5.  SI. 

6 2025  2«^.  <»6, 

294.1  66.6  0.00  0.00.? 

0 0 0,  0. 

COMP AMU 
DIS 

Br IME/t  ri ME/RAMCE/OVEKLAY/ECHQ  MUMBEKS 

2025  2125  202  VIC  6 0 0 0 0 0 0 

CUMMAMU 

ACC 

AKEA/IMfENSlfV 
1000  5 

DAY/  TIME/  TILT 
167  20J0  0 

STOP 


2 STK  tCHCtS  FOUMU  WITH  ARLA  GREATER  THANI000  SC  KM 
ARtA/AZ ImUTH/RANGE 

48  70  26  I 41  ce^  6 - Aome  decieoAe  -ot  oAea,  cent/io^  position 

17  78  3 00  l^5no  oM-igment  coMl6tent 


COMM AMU 
EMI 

M/  HHIfl/  WM/  RMG/  DIRECTIUM/  SPEED/  STCOS/  STCIM 
6 2030  287.  41. 

276. M 62.  1 20.79  0.097 

0 0.  0. 


CUMM4M0 

UlS 


BTIME/C  T I M«- /RANGE /overlay /ECHO  ^UMBERS 

2030  2130  200  VIC  60000000 

command 

BYE 


82 


APPENDIX  C 

Annotated  Command  Structure,  November  2,  1974 

rwc 

HHMM 

2400 

:UHHA'43 

SCO 


COMMAND 

ACC 


At^EA/MTEMSITY 
150  4 

DAY/  TIHc/  TILT 
306  1225  0 


2 VST  LCHUES  FOUND  WITH  AREA  GREATER  THAN  150  SO  KM 

area/azimuth/range 

154  251  IbczZt  1 - nejuj 

540  311  95  cM  2 - new 

DAY/  TIME/  TILT 
306  1228  0 


3 VST  ECHOES  FOUND  WITH  AREA  GREATER  THAN  150  SC  KM 


AREA/AZIMUTH/RANGE 


219 

250 

15  cell  1 

189 

304 

95  celt  3 - new 

322 

321 

94  no  066ignmeJtt 

DAY/  TIME/  TILT 
306  1230  0 
STOP 


Split  oi  celt  2,  6kip  4 in  sequence, 
may  uic  txteA 


4 VST  ECHOES  FOUND  WITH  AREA  GREATER  THAN  150  SC  KP 
AREA/AZIMUTH/RANGE 


225 

251 

It,  cell  1 

179 

306 

95  cell  3 

337 

321 

no  a^6i^nment  (no  motion  jjAom  pfievioub  time) 

222 

34S 

126  cell  5 - new 

83 


1 


CGMMAMD 

Ei'JT 


N/ 

HHMM/  12M/ 

RNG/ 

DIRECTION/ 

SPEED/ 

STCCS/ 

STCi  y 

1 

1225 

251. 

76. 

1 

1228 

25£. 

75. 

30  3.3 

33.  1 

0.00 

0.000 

1 

1230 

251. 

74. 

256.  0 

23.  8 

27.41 

0.80  i 

3 

122H 

3£A. 

<35. 

3 

1230 

3£6. 

95. 

21  5.  0 

96.  8 

0.00 

£ <00.5 

£ 

2 

2 * 

£. 

command 

DIS 

BTIME/^ T1 N -/« ANGC/QVtRLAY/ECHC  KUMHEWb 

123e  13345  2K£  VIC  1 3 0 fc)  45  45  w 2 e 

CuMMANH 

ACC 

Ar<EA/I  gfENSI  FY 
lb£  ^ 

OAY/  TlMt / TILT 
326  1232  e 


A VST  tCHOcS  FQLNO  4hITH  AREA  GRLAIER  THAN  1^0  SC  KM 


AKEA/A2IMLTH/RANC.E 


21£ 

249 

llC-dJUi  1 

194 

3£6 

94cete  3 

364 

321 

93no  oM-igme/it 

193 

35£ 

IZTceU  5 

DAY/  TiMu/  IILT 
3£fc  123A  £ 


3 VST  ECHOES  FOUND  WITH  AREA  GkLATER  THAN  150  iC  KM 
AREA/A2 IMLTH/RANGE  * 

OQ  1 

(meAqe  oi  3 iMXh  unasi^ned  ceM, 
centAO-cd  position  cotui^ttnt  uUjth  old  cett  2) 


192 

249 

TlceJU 

55£ 

315 

‘iZceM 

216 

349 

\27ceJU 

DAY/ 

TIME/  TIL 

T 

306 

1236 

e 

STOP 

84 


1 


2 VST  ECHCES  FOUND  WITH  AREA  GktATER  THAN  Ibid  SC  Kf« 

AKEA/AZIMLTH/BANGE 
ise  2AS  uceZt  1 

558  31  ? <i2CtU  2 


CCMMANO 

ENT 


N/ 

HHMM/  A2M/  RNG/ 

DIRECTION/ 

SPEED/ 

STCUS/ 

STCfM 

i 

1232 

245. 

73. 

26  3.1 

30.  8 

43.72 

0.755 

1 

123A 

245. 

71. 

278.3 

36.  6 

41.52 

0.72? 

1 

1236 

245. 

72. 

274.  1 

37.  3 

38.72 

0,675 

5 

1230 

34S, 

126. 

5 

1232 

352. 

127. 

235.  1 

72.  5 

0.20 

0.000 

5 

1234 

34S. 

127. 

16«.'7 

15.  0 

15. £4 

1.06  0 

2 

1225 

311. 

55, 

2 

1234 

3 15. 

52. 

24  7.  7 

47.  9 

0.20 

0 • 0 2 >i 

2 

1236 

317, 

52. 

24  2.6 

5 3.  3 

20.54 

2.50/ 

2 

0 

2. 

2. 

COMMAND 

DIS 

kllME/c  n ME/KANGE/OVERLAY/ECHO  r.UMuERS 

I23t  1336  20e  VIC  1 2 5 0 0 0 4) 

command 

ACC 

AkEA/I NTLNSirv 
150  A 

DAY/  TIME/  TILT 
326  123H  2 


3 VST  ECHOES  FOUND  WITH  AREA  GREATER  THAN  150  SC  KM 


AREA/A2IMUIH/RANGE 


156 

252 

69ce£e 

J 

551 

315 

9lceZe 

2 

152 

351 

132cete 

5 

\ 

85 


d 


DAY/  TIMt/  TILT 
1126  L2Ae  ?. 
STOP 


2 VST  fcCbUES  FOUND  talTh  tRtA  GKtATER  THAN  15P  SC  KP 

AKEA/AZIMUTF/RANGE 
216  2S£  68 

5A6  32£  Se 

CuMMANO 

EVT 


N/  HHMM/  MM/  RNG/ 
1 1238  25£.  69. 

DIRECTION/ 

SPEcD/ 

STCCS/ 

STCIM 

1 12A0  252.  68. 

265.2 

35.  6 

39.77 

0.67‘> 

2 1238  319.  91, 

259.9 

39.  6 

37.31 

8 ♦ 6 3 5 

2 12A0  322.  98. 

292.  1 

59,  6 

17.93 

0.61V 

5 1238  351.  138. 

29  3,2 

61,  7 

18.19 

0.55g 

21A.7 

0 0 2.  8. 

COMPAND 

DIS 

BTIME/HTIPH/RANGE/CVEKLAY/ECHC 

98.  9 
NUMflERS 

36.72 

1.32  V 

12A2  13A0  282  VIC 

COMMAND 

ACC 

AREA/INTEiMSIT  Y 
152  A 

1 2 

5 

0 0 

0 0 

IJAY/  TIME/  TILT 
34!6  12A2  £ 


3 VST  ECHOLS  FOUND  WITH  AREA  GREATER  THAN  150  SC  KP 
AREA/AZlMLfH/RANGE 


199 

252 

blcM  1 

159 

2F9 

86  80  oAA^wnent 

561 

322 

92cett  2 

86 


- 


04Y/  TIME/  lILT 
STOP 


2 VST  ECHOES  FOUND  WITH  AREA  GREATER  THAN  SC  KP 

AREA/AZlMtlH/RANGE 

lie  25fc  b*tceZt  1 tb&veAal  ceZt6  m&figed,  may  be  n.e^ult 

9AA  311  83ce££  6 - new  in  hiteju^Lty  oh.  Hadan.  powen 

( ^Zuctuatlon 

I SEV  ECHOES  FOUND  WITH  AREA  GKEAItR  THAN  lb'1  SC  KP 

AREA/AZIML  FH/RANGE 
157  3?7  92 

COMMAND 

ENT 


N/ 

HHMM/  AZM/ 

RNG/ 

CIREC  TION/ 

SPEED/ 

STCCS/ 

STCI 

1 

12A2 

2S2, 

67. 

256.7 

3A.  2 

35.21 

2 

.6Hi 

I 

12A5 

252, 

6A. 

25A.  1 

35.  b 

33. AB 

2 

.61') 

2 

12A2 

i22. 

92, 

2A2.5 

6A.  9 

18. A5 

2 

. 5Ao 

6 

1 LAS 

311. 

ou 

e 

2 

e. 

2. 

COMMAND 

uis 


HI IME/=T I MF /RANGE /OVERLAY /ECHO  NUMBERS 

12A5  13^5  222  VIC  122232^22 

Command 

ACC 

AREA/l  DENSITY 
15K  A 

DAY/  TIMT/  TILT 
3K6  12A7  2 


I 

\ 


AREA/AZIMUTH/RANGE 
153  252  6lC.ett  1 

85(2  313  80ce^e  6 

19?  353  130cete  7 - new 


I SEV  ECHOES  FOUND  KITH  AREA  GREATER  THAN  150  SC  KR 

AREA/AZIHUTH/RANGE 
225  325  8A 

DAY/  TIME/  TILT 
306  1251  0 

STOP 


3 VST  ECHOES  FOUND  HITH  AREA  GREATER  THAN  150  SC  KH 

area/azimuth/range 

19A  255  btceJU  1 

9A6  312  eecezi  6 

151  353  l32ceU  7 


I I SEV  ECHOES  FOUND  WITH  AREA  GREATER  THAN  150  SC  KH 

j AREA/AZ (MUTH/RANGE 

257  326  83 

I 


I gfLOWth  Jin.  both  c.eJZ6 


COMMAND 

ENT 


N/ 

HHMM/  A2M/  RNG/ 

DIRECTION/ 

SPEED/ 

STCDS/ 

STCTM 

1 

12A7 

251. 

63. 

250.5 

36.4 

33.63 

0.593 

I 

1249 

252. 

61. 

246.8 

37.8 

33.45 

0.599 

1 

1251 

255. 

60. 

240.6 

39.0 

40.84 

0.576 

6 

1247 

310. 

81. 

346.1 

74.1 

0.00 

0.000 

6 

1249 

313. 

60. 

268.5 

62.2 

69.93 

0.944 

6 

1251 

312. 

80. 

271.2 

39.4 

62.41 

1.099 

7 

1249 

353. 

132. 

7 

1251 

353. 

132. 

180.0 

0.0 

0.00 

0.0011 

0 

0 

0. 

0. 

COMMAND 

DIS 

BTIME/ETIME/RANGE/OVERLAY/ECHO  NUMBERS 

1251  1351  2BS  VIC  1 6 7 0 0 0 B 

COMMAND 

ACC 

AREA/INTENSITY 
150  A 

DAY/  TIME/  TILT 
1253  0 


0 0 


I 

j 

i 


2 VST  ECHOES  FOUND  KITH  AREA  GREATER  THAN  150  SO  KM 

AREA/AZIMLTF/RANGE 
ISS  256  b9ceZl  1 

971  31A  l^cdUi  6 gfiowth  cjOYitiming 


2 SEV  ECHOES  FOUND  WITH  AREA  GREATER  THAN  150  SC  KM 


AREA/AZIMUTH/RANGE 


177 

267 


300  77 

327  85 


DAY/  TIME/  TILT 
306  1255  0 
STOP 
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2 VST  ECHOES  FOUND  WITH  AREA  GREATER  THAN  150  SO  KH 


AREA/AZIMUTH/RANGE 
241  25S  STcetZ  ; 

1123  315  B§ceJtl6 


bo^  cetti  gftotacng,  ceZt  1 acceZeJuvUng 


1 SEV  ECHOES  FOUND  talTH  AREA  GREATER  THAN  150  SC  KH 

AREA/AZIMUTH/RANGE 
471  316  78 

COMMAND 

ENT 


N/ 

HHHM/  A2M/ 

RNG/ 

DIRECTION/ 

SPEED/ 

STCOS/ 

STCT 

1 

1253 

256. 

59. 

236.1 

39.9 

39.61 

0.555 

1 

1255 

259. 

57. 

231.4 

41.9 

41.75 

0.624 

6 

1253 

314. 

79. 

260.7 

43.3 

58.90 

1.076 

6 

1255 

315. 

80. 

248.2 

41.6 

57.34 

1.015 

0 

0 

8. 

0. 

COMMAND 

DIS 

BT I HE/ETl HE/RANGE/OVERLAY/ECHO  NUMBERS 
1255  1355  208  VIC  1 6 0 0 

COMMAND 

ACC 

AREA/INTENSITY 
150  4 


DAY/  TIME/  TILT 
306  1257  0 


2 VST  ECHOES  FOUND  WITH  AREA  GREATER  THAN  150  SC  KH 

AREA/AZ1MLTH/RAN6E 
205  260  54cete  1 

1019  319  Bfcete  6 


1 SEV  ECHOES  FOUND  tolTH  AREA  GREATER  THAN  150  SC  KH 

AREA/AIIMUTH/RANGE 
293  326  80 


1 


r 


OAY/  TIME/  TILT 
3i6  1259  0 


2 VST  ECHOES  FOUND  klTH  AREA  GREATER  THAN  150  SC  KM 

AREA/AZIMLTH/RANGE 
260  262  52ce^e  7 

1000  320  79cete  6 


1 SEV  ECHOES  FOUND  WITH  AREA  GREATER  THAN  150  SC  KM 

AREA/AZIMLTH/RANGE 
231  327  78 

DAY/  TIME/  TILT 
306  1300  0 

STOP 


2 VST  ECHOES  FOUND  WITH  AREA  GREATER  THAN  150  SC  KM 

AREA/AZIMLTH/RANGE 
278  262  Sleeve  1 

1003  322  68ce££  6 


1 SEV  ECHOES  FOUND  WITH  AREA  GREATER  THAN  150  SC  KM 

AREA/AZIMUTH/RANGE 
239  327  77 

COMMAND 

ENT 


N/ 

HHMM/  AZM/ 

RNG/ 

DIRECTION/ 

SPEED/ 

STCCS/ 

STCl  1 

1 

1257 

260. 

54, 

22  8.9 

44.4 

40,65 

0.677 

1 

1259 

262. 

52. 

22  7,1 

46.9 

39.60 

0.685 

1 

1300 

261. 

51. 

22  6.7 

48.  3 

42.70 

0.680 

6 

1257 

319. 

80. 

23  7.8 

54.6 

57,82 

1.402 

6 

1259 

320. 

79. 

236.7 

59,5 

54.73 

i.3n 

6 

1300 

322. 

80. 

234.2 

65.0 

57.53 

1,359 

0 

0 

0. 

0. 
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COMMAND 

OIS 


BTIME/ETIME/KANGE/OVERLAY/ECHQ  NUMBERS 

L3ee  1AB0  280  VIC  1 6 0 0 0 0 0 0 E 

COMMAND 

ACC 

AREA/INTENSITY 
150  4 

DAY/  TIME/  lILT 
306  1302  0 


1 


2 VST  ECHOES  FOUND  klTH  AREA  GREATER  THAN  150  SC  KM 

AREA/AZIMUTH/RANGE 
278  267  5\C.iJUi  1 

993  324  81 CZJU  6 


1 SEV  ECHOES  FOUND  ViITH  AREA  GREATER  THAN  150  SC  KM 


AREA/AZIMUTF/RANGE 
207  329  77 

DAY/  TIME/  TILT 
306  1304  0 


2 VST  ECHOES  FOUND  WITH  AREA  GREATER  THAN  150  SC  KM 

AREA/AZIMITF/RANGE 
266  269  5BCZlt  1 

1048  326  82  cett  6 


1 SEV  ECHOES  FOUND  WITH  AREA  GREATER  THAN  150  SC  KM 

AREA/AZIMLTH/RANGE 
202  331  78 

DAY/  TIME/  II LT 
306  1306  0 


1 SEV  ECHOES  FOUND  WITH  AREA  GREATER  THAN  150  SC  KH 

AREA/AZIMUTH/RANGE 
205  332  79 

COHHANO 

ENT 


N/ 

HHMM/  A2M/ 

RNG/ 

DIRECTION/ 

SPEED/ 

STCCS/ 

STCr 

1 

1302 

267. 

51. 

229.3 

50.2 

50.80 

0.701 

1 

1309 

269. 

50. 

22  2.3 

51.  8 

99.67 

0.683 

1 

1306 

270. 

99. 

221.0 

52.  8 

CD 

• 

00 

0,670 

6 

1302 

329. 

81. 

231.5 

69.2 

55.99 

1.292 

6 

1309 

326. 

82. 

229.2 

72.  6 

53.93 

1.235 

6 

1306 

327. 

61. 

22  9.0 

73.3 

52.15 

1.190 

0 

0 

0. 

0. 

COMMAND  i 

DtS  i 

BTIME/ETIME/RANGE/QVERLAY/ECHO  NUMBERS 

1306  1906  200  VIC  1 6 0 0 0 0 0 0 0 < 

COMMAND 
ACC 

j AREA/INTENSITY 

150  9 

DAY/  TIME/  TILT 
306  1308  e 


2 VST  ECHOES  FOUND  WITH  AREA  GREATER  THAN  150  SC  KM 

AREA/AZIMLTH/RANGE 
290  267  97cea  1 

1052  322  BlcM  6 


’I 


DAY/  TIME/  TILT 
386  1310  8 
STOP 


2 VST  ECHOES  FOUND  talTH  AREA  GREATER  THAN  150  SC  KH 

AREA/AZIMUTH/RANGE 
285  278  1 

108A  328  eacejU  6 

COMMAND 

ENT 


N/ 

HHMM/  AZM/ 

RNG/ 

DIRECTION/ 

SPEED/ 

STCDS/ 

STCTi 

1 

1308 

267. 

47. 

221.1 

52.9 

51.53 

0.710 

1 

1310 

270. 

46. 

220.9 

53.2 

50.72 

0.69J 

6 

1308 

328. 

81. 

229.0 

72.6 

50.14 

1.149 

6 

1310 

328. 

80. 

229.8 

69.8 

48.96 

1.143 

8 

0 

0. 

0. 

COMMAND 

DIS 

BTIME/ETIME/RANGE/OVERLAY/ECHC  NUMBERS 

1318  1410  208  VIC  1 6 0 0 0 0 0 8 8 

COMMAND 

BYE 


APPENDIX  D 

Computer  Program  Listing  for  ECHOPRED 

CJf'f'CN  X8(9.'-'!  ),Xfc(  SLI2Z  3)  ,PHI4I22  !|  ,PHIE  ^2a^■| 

♦,  ISVGT(4n,GTMlN(4:)  ,GTMAX(4U  ),SUEN6I  U).NH4M|  1 ? ) , lECA  Z I 1 2. 90  , 

*lECRNGnn;,9C.1),  I ECA0,IGTLENG,HAP(  120)»MEN.COFI  I Z. 4. 10)  .gate  (2201  , 
♦lOU,  Br.XTP.lhSKC  19) 

OIPENSICA  IPC(NT(2  e).LNSAV(4?) 

D1  PENSION  NCB(  l,r  ),  «»X  ( 12),  BX(  12  ),4Y(  1C),BV(  )>*  ) , I ENC  1 1 0 ) 

OIPENSIUN  STE2(U),SCE2(U) 

OIPENSION  SMT(l2),Swr2(  U’),SMX(  1J  ) , SWY(  1 ?)  . SMTX  ( 1 . ) , £bT  V(1  0 ) 
DIKENSION  SWIW  ) ,FT1ME(  1.-') 

DIPENSION  JCIRCS,;),  K Y M ■>(»)  . AL  INE  ( 24 ) . I EC  ( 11 ) , SERTP  ( 53  ) 

DI  PE  NS  ION  SNi  (52  ),  SNi(  3’I),S'J3(52),SPTP(  50)  ,S0IST(5'.)  ,CP(9) 

01  PE  NS  ION  JEC(13  ),  *;Rt4(  l.M.STORM(t) 

OIPENSION  lECINT(r) 

OIPENSICN  IKNT16) 

OIPENSION  ELPSI3) 

real *8  SW,SWT,Swr?,SbX,SWTX.SWY,ShTY,XT,  YT 
INTEGER  GTMIN.GTPAX, T0.T 1,T2, bOXTM 
INTEGER*?  lECAZ 
CHAR  lECRNG 
char  GATE.IHSK 

DATA  CP/*P  ’S  S'SES'E  '.‘NE'.'N  '.'NW'.'U  •/ 

DATA  CAIR.OBYE, CENT, CCISiQOEL.QTWC.QPOS.QwEE.UIGN 
♦/  3EflIR,3HeYE,3EE.\t,3HDIS,3H(.EL,  3ETWC,3HPOS,3HbHE.3HIi,N  / 

DATA  CKEY.QACC, CGCC, CBQC/ 3HKE Y, 3PACC, 3HGC0 .3HB0C / 

DATA  ST0RP/24E  LGT  MCC  STR  VS  I SEV  EXT/ 

FAD1R(X,Y)=RT D*ATAN2(Y,X )*180- 

HUUR(K)= FLOAT (K/1?  ) 4FL0 A T ( NO n( K,  100)  1/62. 

JHEPPIT)  = T*6?'.-»4,'*IF  1X(T  ) 

1107  FQKPAT(2I5,2F5.2 ) 

1108  F0RPAT(1X,I2,  I5,2FS  . ) 

1706  FORPATd  115) 

6706  FORPATdX,!  115) 

MXNCE=13 
MEN=  1 '2 
MXEP1=PXNCE*1 
KTC1=57. 2957795 
ALN2=^.693l 
DTR=3. 0174532925 
04GP=4.?/3. 14 15927 
C READ  STATION  NAMES  ONTC  DISK 
REUINC  2 

1731  REA0(5,17?2 )TN1,TN2, TN3,XA,YA 

1732  FUBPATl3Ai,,3X,2F6.2) 

1F(TN1.EC.4HTAF  )GC  to  1733 

CONVERT  NP  TC  KP  FOR  AIRFOPTS  AND  STORE  ON  DISK 
XA=XA#1,e52 
YA=YA*1.B52 

WRITE(2, 1732)TN1,TN2,TN3,XA,YA 
GO  TC  1731 

1733  END  FILE  2 

c... enter  overlays  CNTC  nsK 

RFNINC  3 
DC  1794)  1*1  ,2 

1706  REAC(5,17e7)XB(  1 ),Y8t1),XE( ll.YEI  1) 

1767  FORPAT(7F10.4) 

IF(XB(1) .EQ.999. IGC  TO  1TBS 
L=1 

CALL  SETLINIL) 

WRITEO,  1787)XB(  1),YP(  1 ),XE(  I l,Yt(  1 1 . PHIP  ( 1 ) .PHIE  ( 11 , <L(  1 1 
IE(L.EC.2)HRITEI 3,  17e7)XB(2), V6(2),XE(2).VE(2),PHI6(2I ,PH1E(2I.SL( 
♦2) 

GO  TC  1736 
1788  ENU  FILE  3 
1790  COM  INUE 

C DEFAULT  TIPE  WEICHT  CONSIANT  IS  32  MINUTES. 

TWC*-ALN2*0.5 

C DEFAULT  GROUNC  CLUlTER  DISTANCE  FOR  ECHO  CONTOURING  IS  29  KM 

IGC=2': 

C....IF  NC  TIME  IS  AVAILABLE  FROM  BFNDIX  DISPLAY  ROXTM  * *. 

BDXTP*(1, 

C REAC  CCPMANC 

DO  1L20  1*1,10 
JEC(  I )=2 
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00  1S2.0 

00  1f2(>  K«1.1C 
1022  CUFI 

, . 00  Un 

1010  M&P(I)x3 

1 = 1 

1035  00  lO^ig  I = I,MXNCE 
lENCm=3 
1040  NOB  I I )=3 
21f»P  MRITE(6i21(ri) 

2l0l  F0RP4H*  ^COPPANC*  I 
C***360  C6PFKDENT 

REAO  (bflt  21051CMC 
C**«0*** 

URITE(6.2ie6)CMC 
2100  FORPATUX.24A3) 

2105  FORPAH2AA3) 

IFCCPC  .EC.  CENT1GC  TO  ll^v' 
IFCCPC  .EC.  CGCOGC  TO  12;ir 

,EC  

.EC 
EC 
.EC 
.EC 


IFCCPC 
IFICPC 
IFCCPC 

ifIcfc 

IFICPC  

IFICPC  .EC. 
IFCCPC  .EC. 


CTWOGC  TO 
CBYE1GC  TO 
CIGMGC  TO  I5w. 
GCEL1GC  TC  163:^ 


CAIR1GC  TO  173. 

CCIS1  CO  TO  lai- 
CFCS  »Gt  TO  19.V 
IF1CPC  .EC.  CACCIGC  TO  2?;<* 

IFiCfC  ,EC.  CWFEIGt  TO  3'^0? 

IF(CPC  .EC.  CeCCIGt  TO  5iii:. 

UR1TEI6.2112)CMC 
2110  FORPAT»lX,A>,*  INV/LIC,') 

GO  TC  2H3 
C 

C TIFE  WEIGHT  CONSTANT 

C ENTER  ECHC  CfiSERVATlON 

1100  WRITE(6, 11P1) 

1131  FORMAT!*  N/  FHMM/  A2M/  RNG/  niRECTION/  SPEED/ 
00  1 Vi2  1 = 1 ,KXNCE 
1102  JECm=T 

C RE  AC  ECHO  CBSERVfTlON 

1105  JE=r‘ 

C***360  CEPENCENT 

READ  (53,  11071JE,  JTF,  £AZM,£R.^G 

WRITE16. 11P81JEt JTF, EAZM.ERNG 
J=1ABS(JE) 

IF(  J .EG.  P )GC  TC  21  '’3 
C PERFCRF  INPUT  VALICITY  CHECK. 

IFU  .GT.  F 

IF1FCCIJTH, 

IFIEAZF.GT,  _ 

CCFFUTE  TIME  IN  FOURS*  X AND  Y COORCINATES. 


SIOCS/  STDTM'I 


irLXr  INKUI  VALILIIT 

.GT.  MEMGC  TC  1160 
;cijth,19pT  .ge.  abigo  to  ik 
iZF.GT,360)GC  TO  1170 


1115 


ETIF=FCUBCJTH) 

AZF=EAZM*nTR 
XT=DCCSI AZM)*ERNG 
YT*OSIM  AZh)*ERNG 
CHECK  FCR  NEW  ECFO. 

L=FAP(J) 

IF(L  . GT  . 0)  GC  TO  1130 
IFIJE  .LT.  eiGO  TO  1133 
NEW  ECHC.  ALLCCATE. 

00  1115  L=1*PXNCt 

IFIIENCILI  . EQ.  CC  TO  1125 

CONTINUE 

ALL  hCRKlNG  ECHCS  USED. 


WRITEI6. 11201 lENC 
C***36i  CE PENDENT 
1120  F0RPAT1<UNABLE  TC 


}RPAT(<UNABLE  TC  ACCOMODATE  .<EW  ECHO.' 
*•  CELETE  CNE  Of  THE  F0LLUmIN6*/1(I5) 

GO  TC  2100 

; start  new  ECHO  WITH  Index  l 

1125  NO0(LI>I 
IENC(LI« J 
MAPUI-L 
SOE2IL)«0. 

|tE2|LI-(I. 

SMUI-I.A' 

swtTli=etih 
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SWT2  U)»ETlf'**2 

SMXIU>XT 

SWTX  «l»»XT*ET  IM 

SWVU»*YT 

SWrY(L»*ETIK*n 

FTIKE(L)»ET IK 

GO  TC  1135 

IFiFTlKElLl.GT.ETIM  CTn*ETM^24. 

CCKPUTE  WEIGHT  F4CTCR 
W=EXF(TWC*(  ET  IM-FT  IMf.d.  » ) ) 

TEST  FCR  CeSERVATICN  DEtbTinrn 
IF  ( JE.GT.«»  ) GC  TO  1 Hi? 
atLETE  PREVIOUS  OESFRVATION 
IFIETIK  ,LE.  FTlKElLllGO  TU  1135 
WRirE(6. 113AI 

FORKAK*  YOU  CAM  CELETE  SUMETHH4G  THATS  NOT  THERE*) 

GO  TC  210 
NUB(L)=  NCB(L)-1 
IF(NCP(L ) .LE.  i )GC  TO  1157 
SMIL)=SW(L)-W 
SWT  ( L ) = SWT(  L)-W*ET  IM 
SWT2m-SWT2CL)-W*tT  IM**2 
SWXU  )=SWX(  U-W*XT 
SWYUl^SWYI  L)  -W*YT 
SWTX(L)=SWTX(L)-W*IT  IM*XT 
SWTY (L)=SWTY( L)-W*IT  IM*YT 
GO  TC  1135 
iENcm=a 
MAPI J)=3 
GO  TC  1175 
AUL'  NEW  CflSERVATICN 
DT=Aes(EriM-FTlKc(  L)  I 
IF  (CT  .EC.  7.)GC  TO  1175 
IFINCP(L)  .LT.  2)GC  TO  IKS 
CCKPUTE  ERRCRS 

PTIKE=( AXIL )*(XT-BX( L) )♦  AY( L )*( YT-DYIL ) ) ) / 1 AX  I LI **2+ A Y I U **2 ) 
PX=AX  ID  *PT  IKE+eXI  L ) 

PV=AYILI*PTIKE+eY(L) 

ERDST=SUPTI IPX-XT ) PY-YT ) »*2) 

E.^TIK  = ABS(ETIM-PriNE  ) 

RERDS=ERCST/CT 
RERTK=ERTIK/CT 
STE2(L)=  ST  E2IL  ) ♦RIRTM*'*2 
SOE2  (D  = SCE2(L)*REFOS**2 
:iOb(L  l = NCB(  Ll-H 
SWIL)=W*SW(L)  ♦ 1. : 

SWT  I L ) = W*SWT(  L)  ♦ tT  |M 
SwT2  (L)»w*SWT2(  L ) ■*  FTIM  * ET  IM 
SWX  (D>W*SWXIL)  * XT 
SWYI L)=W*SWY( L)  ♦ YT 
SWTX  (D=W*SWTXI  L I + eTIM*XT 
SwTY  IL)»W*SWTYI  L ) •*  FTIM*YT 
FTIKEU)  = AKAXUFTlKEa)t  cTIMI 
CCKPUTE  NEW  APPROXIMATING  PCI  Y ANC  SPEEO/DIR. 

DEN=SWIL  )•SWT2(  L )-SW  TIL  )**2 
VV=I SWIL )*SWTX( L I-SWTIL )*SWXI  I ) l/OEN 
UU=ISWlL) *SWTYI L )-SWT(L )*SWYI  I I l/OEN 
AXIL  )=WV 
AYIL  l=UU 

BXILKXT-VV*ETIM 
flYILI=YT-UU*ETIM 
AOIR^FAOIRIVV.UU  I 
SPC=SCRT IUU**2>VV**2 I 
ANCE^NCeiL) 

STCCS»SCRTISCE2ILI/Af,OEI 

STDTK=SOHTISTE2IL)/ANOE) 

E0PH  = AKAX1ISTCCS,1'.  I 

BT2A»HCUR|3CXTM I 

IFIET  IK.GT.BT  2A  IBT  2A=BT2*»*2A. 

IFI£TIK,NF,eT24IGO  TC  1162 
DO  1 ISa  JK  l.MXNCE 
A=CCFI1,  1,J1) 

B=CCFI1,3,J1  I 

ASU0  = SCRT[I  A-YT  )**2*IB-XTK*2I 
IFI  J1.EC.1IGC  TC  1 ITS 
IFIaSUB-akIM  1179,  llEI",!!"'. 

AMlN=ASUe 
JECILKJ1 
CONT INUE 
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IF(A^*I^.LE.ECPH)GO  TC  1162 
111116(6,1181) 

1181  FUBPATjfaCHECK  CENTRCID  OF  THIS  6CFO.  THIS  OBSERVAIICN  DbLETEO'l 
GQ  T L 1135 

1182  WRITEC6,  ns^lACIR.SPC.STCDS.STQTM 
1150  FQRKATI19X, F) 2. 1 , F7 . 1. F9 .2t F9 .3 1 

IF(SPC.GT.120.)GC  TO  1185 
GO  TC  11 ?5 
lAPUT  ERRORS 
HKITE(6, 116 1) 

FQRRATC'  echo  NC.  IOL  LARGE*) 

GO  TC  2W* 

WRITE(6,1166) 

FORPAT(*  SCPETHING  SCREWY  AdOUT  THAT  TIME* I 
GU  TC  21 ’0 
WRlTEie, 1171) 

FOHPATC*  AZIPUTH  CCES  <^OT  COMPUTE*) 

GO  TC  212? 

WRITEI6, 1176) 

FORPATI*  PREVIOUS  tCFO  AT  SAME  THE.  THIS  OBS  IGNCREO*) 

GO  TC  1105 
WR1TEI6, 1186) 

FURPATf*  ECHO  SPEEt  GREATER  THAH  125  KM/HR,  CHECK  THIS  CO.'l 
GO  TC  11 V 5 

C...SET  GRCUNO  CLUTTER  DISTAMCE  FOR  CONTOURING 
1200  REACI50,  1205MGC 
1205  F0RPAT(I3) 

WRITE(6| 1205) ICC 
GO  TC  210®' 

1300  WRITE(6, 13051 
1335  FORPATC*  HHPP') 

C***36C  CEPENCENT 

REAC(50, 1107)ITWC 
WRITE«6.17P6)1TWC 
TWC=-ALN2/HCUR(  ITWC) 

GO  TC  210.' 


'mi 

1)65 

1166 

1173 

1)71 

1175 

1176 


11> 

IK 


: BYE.  STOP. 

U00  call  TPKEWI  10  ) 
1FWTP«11) 
TPREWI 1) ) 


call 

call 

STOP 


C 
C 

C*A*360 
1500  REACI 


IGNCRE  TEXT 


CEPEAOENT 
ACtr-0,5l«l5)ALlNE 
WRirE(6.2106)ALiNE 


IFCALINcI 1) 
GO  TC  1530 


.EC.  QKEYiGu  TO  2 nr 


c delete  an  echo. 

1600  00  U'S  Hl.MXEPI 

1605  IEC(i)*i3 
WR1TE(6, 1606) 

1606  FCRPATC  WHICH  ECHCS*) 

C***36B  CEPENCENT 

REAC  (50.  1706) lEC 
C**<360  DEPENDENT 

WRITE(6.67e6)IEC 
DO  16le  Hl.PXEP) 

K=IEC(I) 

IF(K.LE.3.0R.K.GT.NEN)G0  TO  l6l0 
L=>PAP(K) 

IF(L  .LE.  2)GC  TC  16  10 
1ENC(L)»0 
HAP(K)>0 
1610  CONTINUE 
1615  K>0 

00  1620  H1.PXNCE 
IFdENCd)  .EC.  0)CO  TO  1629 
K=K+  1 

IEC(K)>|EN0(I ) 

1620  CONTINUE 

IF  IK  .EC.  0IGC  TC  16  30 
WRITE(6,  1625)(IEC(  D.I-I.K) 

1625  FORPATI***  ACTIVE  ECHOS  •.10IA) 
GO  TC  210A 
1630  WR1TE(6, 1635) 

1635  FORPATI***  NC  ACTIVE  ECHOS') 


98 


GOTC  21^-? 


C.  PRECICT  ECPC  PJTE  CVER  AIRPORTS 

J7ii3  MR1TE(6,  17^13) 

1703  FORKATI*  BTIKE/  ETIME/  STC  OEW/  ECU)  NO. S' I 
DC  17.:5  I=1,PXNCE 
1705  IECm»«! 

C***360  CEFF.KCENT 

READ  I5P,  17P7)  JTP.J  lE.SDf),  ItC 
C***«**« 

C**  FORPAT  17^6  FAS  PUL T I PL E REFS. 

1737  FOBPATI2I5.F5.?,11 15) 

6707  FORP4T(lX,2l5,F5.'Sni5) 

WRlTEl6.67B7)JTe,JTE,SUC. lEC 
BT24=FCUR(JTB) 

ET2<»  = FCUB(JTE) 

DO  1730  JE=1,MXNCE 
J=1EC(JE) 

L=LEGEM  J) 

IFIL  .LE.  B)GC  TC  21 
«RITE(6,n0«»)J 
1708  FOKPATt‘CECFC'I<i) 

C CCPPUTE  BCUNCARY  OF  bOX. 

NoeL=AceiL) 

ance=ncbl 

RTP=FTIPE«L ) 

IF(BT2A.GE.RTM)CC  TO  AtT? 

BTIP=ET2A423. 

ETIP  = ET2A^2<i. 

GO  TC  AiMfl 
BT1P  = BT2A 
CT1P=ET2<( 

1F(ET1P.LT.BTIP)ET  IM=ETIM-'24. 

^010  EOPH=APAX1(SCO*SCRT(  SCE2(U/AiWE).  1C.) 
ETPH=«PAX1(SCRT (ST £2(L )/ANUE). 5.1 ) 

AXL=AX(L ) 

BXL=eX(L) 

AYL=AY(L) 

BVL=EY(L ) 

RSSC=  1.0/  ( AXL*AXL+  'YL*AYL  ) 

EC=SCfiT( RSSC)*ECPH 

Wd=AeS((  eTIP-BTP)*/XL*cC) 

WE=ABS((  ETIP-RTP)*AXL*EC) 

PB= AXL*8T IP+BXL 
PE=AXL*ETIK-fBXL 
PXfe=APIMlPe-we,PE-WE) 

PX)  = APAx1(PE  + WB.PE+WE) 

WB=AESM  BTI  P-RTP  )*<YL*cC  ) 
WE=AES((ETIP-RTP)*iYL*£C) 

PB=AYL«8T IM+BYL 
PE=AYL*ET IM4BYL 
PY0=AMNlIPB-WEtPE-WE) 

PY1  = APiX1(Pli4W8,PE4ME) 

NT=C 

NCK=C 


REAC  TABLE 


TOWNS  ANC  REPORT  TFOSE  TOWNS  IN  PAIH  CF  ECHO 


1710  WRITE16, 1713) 

1713  FORPATt*  TIPE  AIRPORT  DI  ST 

REWINC  2 

1733  REAC(2,1732,ENC=  174.5  )Tn1,TN2,  TN3,  XA,  YA 
E CHECK  FCR  LCCATICN  IN  BOXED  AREA. 


FTIM  ITIP') 


: CHECK  FCR  LCCATICN  IN  BOXED  AREA. 

NCK=NCK+ 1 

IF(XA  .LT.  PX0  .LR.  XA  .GT.  PXD30  TO  I7AB 
IF(YA  .LT.  PY3  .CR.  YA  .GT.  PYDGO  TO  17A* 

: CCPFUTE  ERRORS 

PTIPE=(AXL*(XA-eXL  1+ayL*«YA-BYL11*RSSC 

PX»AXLAPT IME48XL 

PY=AYL*PTIPE4BYL 

OIST  = SCRT((  PX-X  A )**2  4«PY-YA)**2) 

Or»ABS(RTP-PT  I)4£  ) 

1737  IFICIST  .GT.  ECPF*CT  )GO  TO  17A3 

IF(PTIPE  .LT.  BTIH  .CR.  PT IME  .GT.  ETIHIGO  TO  ITAl 
1730  NT«NT4l 

S'lUMl^TM 

SN2(NT)=TN2 
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1 


SM3(M»«TN3 
SOISTIKT  »>CIST 

Jl}IR|KT)>C4PP*AT  PX-XA,PY-VA  i*9,9 
SPTP(Mi«PTlPE 
SfcRTPIKT »»Cf*ETPP 
GO  TC  173*^ 

C.  , ^ ALL  ST  AT  I CNS/ A IMPORTS  CHECKED. 

W45  IFINT.GT.OGC  TC  KS-  . 


..  WRITEI6.  IU6I 
1146  FOKHATI*  NC  E 
60  TC  17  3*1 


NCCUNIERS  PnEOICTEO* I 
NS  EV  ENCOUNTER  TIME. 


C slrt  static  . _ . .. 

1753  CALL  lSCPTiSP7K,KV,NT) 

00  1775  I»1,NT 
J«KYI  I) 

JPTH*JHHPPISPTHU»  ) 

JT1»JFHHPISPT  M|J)>SERTM(J)  I 
JT2=  JFHHPISPTMI  J KSEHTMI  J » I 

i??f 

WRITE  16. 1777INT 

1777  FORPATIlFP, 14.*  ENCOUNTERS  LISTED* I 
1783  WRITE(6t 1785) 

... I 

C 

1785  FORHATI//// ) 

GU  TC  2131 

C...FINO  eCUNCARV  POINTS  OF  CONE. 

COMPUTE  LCCATICN  CF  POINTS  ALONG  PATH  UF  ECHO  AT  RTIM  aNO  EIIN 


.DISPLAY  CCNTCUREC  ECHCES  WITH  aARNING  AREAS  ON  6EN01X  SYSTEM 


1 


WRITEI 


1905) 


16.1! 

1855  FORMAri»’eT IME/ETIHE/RANGE/OVERLAV/ECHO  NUMBERS*) 

irug.OYLY,  lEC 

IRNG.UVLY, lEC 


REACIfJ,  ie);i)JT8*JTE,l 
1810  F0RMAT(3I5|2X.A3,1 1151 
WRITEI6.  13K)JTe.JTEtl 
1GTLENG>  IRNG/2i;<' 

00  18.?')  LINC»1,4: 
ISVGTILINC)*-) 
LNSAVILINC)«0 
GTMAXUINC)-'* 

1839  0TMIMLINC)«3 

OC  len  LINC»I,2H2 
YB(L  INC)»3. 

XEILINC)-^. 

YE(LIAC)«3. 
PH18UINC)«'^. 

PHIl  (LINO-?. 

SUL  I NO -I. 

1611  XB|LINC)«0. 

00  1650  JE-1.MXNCE 
J'lEClJE ) 

L«LEGEM  J) 

IF|L.LE.0)GC  TC  1850 


10  1652 
R 

1852  lECHijECU) 


WRiTE(6| 1853)  J 
1853  FORMATI*  YOU  CANT 
GO  TC  1853 


PRECICT  from  one  OB.  FCHU*.I3.'  IGNCKEO. •) 


iFIIECH.GT.eiGC  TO  1848 
OU  1651  I»1|1B 
IFtCCFn.1tI).NE.t;..CR.C0Fllt3.I).NE.e.)Ga  TO  1851 
i ECH« I 

PX«AXIL)*FTIMEIL)^6XILI 
PY»AYILMFT  IMEiLHEVtL) 

COM  1.3,  lECHl-PX 
C0FI1,1,1ECH)»PY 
C0F(2,1. IECHI-10. 

C0Ft2,4l IEc^)«1e. 

NHAP(!eCH)«2 

SLENG(IECH)>45. 

JECILl-lECH 
CO  TC  18AB 
1851  CONTINUE 


1849  ?8UW^’g^‘’ii!NTCUR  FOR  ECMO*. I?l 


GO  TC  185P 
1848  CALL  ECHCIIECH) 
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1650  CUM  IhUE 
LINE*  » 

I NT* 56 

DO  I£I6  1*4.9 

. GATEII4|3)*1>-SK(  I«10> 

1816  GATEm«lHSK(  I) 

GATEn^l«JTe/1Cr)i; 

jTe/1  n,ij» 

CATtI  J2»*l»ccl  JTe,  1 c»/n 
GATE(13I*KCC(  JTE,1  I 
GATEnei*I6TLENG 
Br24*£CCIRUTe  ) 

ET24*£CU«UTEI 
00  1E22  JE>  IfMXNCE 
J*1EC(JE ) 

L*LECEM  J) 

IE(L.LE.9)GC  TC  19?<? 

N06L*NC8«L1 

IHNCeL.LT.PIGC  TO  1^23 
lECH* JECIL1 
ANCE*NC0L 

EUPH»i>'AX1(SCRT(S0E2(L)/AN0E).  le.  ) 
RTP=FTIHE(L ) 

IFIBT24.GE.I»TM)GC  TO  1814 
BTIP*eT24*24. 

ETIP*ET24*24. 

GC  TC  181S 

1814  BTir'*PT24 
ET1P*ET24 

IFIETIP.LT.BTIMIET  IM*ETIM«24. 

1815  BT*AeSleT IH-RTM) 

ET=AeS(ETIR-RTM) 

ETPH*«PAX1ISCPT  (S»T  C2(L  l/ANUtI.il.1  I 
IFjtTFF.GT.ABSI ET-ET  nETPH*ABSI ET-BT  »/2. 
OX|  = AV(L  >*(  BT  1H-ETFHI4BYIL  ) 

0V1*Axa  )*I  PT  IM-ETPHMBXIL  » 
ax2* AVIL  »♦( ETIN4ETPHI4BYIL I 
QY2*AXIl  >*(ETIM«ETPH»4bXfL) 
xy*(CY2-CY1 )/(CX2-CX II 
CALL  ELLlPS(XP.TANb,XK.ELP6t IFCHI 
IF(ECFF.GT,XK)ECPH*XR 

RAC*SCRT  MECP»-*eT4XK  >••2/1  TAN»*TANW4  1 . I ) 


RXl*CXl4RA0 

RX2=CX1-RAC 

RYI*T4NW*(RX1-QX 1 1 4QY1 
RY2=TANW*(RX2-0X 1I4QY1 


Xl-RAC 

4NM«(RX1-0X1)40Y1 


R4D*SCRT  ( (ECPH4ET4XK  1**2/ ( TA9«i*TANh4  1 . )| 

SX1  = C)«24R»C 

SX2*CX2-RAC 

SY1*TANm*(SX1-(3X2I4QY2 

SY2* TAKW *15x2-0x2  I +0Y2 

LINE*lINE4l 

XB(L INEI*RX1 

Y0IL1NEI»RY1 

XE(L1KE)*RX2 

YE»LINE)*RY2 

CALL  SETLIN(LINE) 

LINE*LINE*1 
XB(LINE)«RX1 
Y3ILINEI=RY1 
XE(LINE)«SX1 
YE(LINEI*SYI 
CALL  SETLINILINEI 
LINE*LINE4l 
XaiLINE)*SXl 
YBILINEI»SYl 
XEIL  INEI*SX2 
YEIL  1NEI*SY2 
CALL  SETLINILINEI 
LINE-LINE41 
XBIL1NEI-RX2 
YB(LINEI*RY2 
XEILINEI«SX2 
YEILIKEI-SY2 
CALL  SETLINILINEI 
1820  CONTINUE 

IWARN«LINE 

CHECK  IF  CVERLAY  IS  TC  BE  CISPLAYEU 
REWIND  3 
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MXNCLN=LIKE 

IF«CVIY.EC.3»-  »GC  TO 

IQVLV=1 

LINE  = II.NE*1 

IF(CVLV.EC.3HST4  MCVLY*2 
GO  TC(ld25. ia35).ILVLY 
1»25  REAC(:,lB33tENC»U  '5  )XttlLlNE) 

183e  FOPH4T(7F10.^ J 
GO  TC  182? 

U35  READ  n,l  330, 6ND«  IX0(  L INE  U YB(  L INE  1 tXE  ( LINE  ) t YE  ILI  hE  ) »PHI  3 ILINF 

*1 ,PHIE(LINE)»SL(LINE  » 

MX^CL^*L INE 
LINE«LINE-»1 
GC  TC  1335 

134J  CALL  LASCRTILNCNT,  IPClNT.MXNOLNI 
ILNE=” 

DO  3 L*1»LNCM 

CALL  GETLINI ILNE.L, IPOINT ) 

LNSAva>  = lLNE 
3 CCMIAUE 

C... START  CC  LCCP  CF  AZIMUTH 
DO  bC.  Il«l,36e 
IAZ=II-1 

CCAVERT  TC  PATFMATICAL  CEGREESP 
IF(  I AZ.LE.9  ;»  ANGLE*9>IAZ 
IF( I AZ.Gr.90)ANGLE-AE0-IAZ 
A'4GLE»ANGLE*CTR 
DO  Ub  lQT*20,22<t, 

GATEfIGT)«P 
146  CONTIAUE 

GATE  m>IAZ/10J 
GATE<2)»RCCnAZ,  V-  )/l4 
GATE(3l>KCO(  lAZ,  IC  ) 

DO  59  L»1,LNCNT 
58  LINE*lFCINT(LNSAV(Ln 

IF(xe(LIAE).EC.999.)CO  TO  59 
INT»8 

IFILIAE.LE. IWARN)IAT=56 
IPHie=FHlPILINEl4.£ 

IPHI  = FHE(LINE»*.5 
JAZ»iaz 

IF( JAZ.LT.1PH1B1GO  TC  59 
PH1DFF«AES( PhIBILlAE  »-PHl£(L iNE 1 ) 

IFCPHICFF.GT.l.ANO.PUDFF.LT.  179.  1GU  TO  99 
IRNG=GTMIN(L) 

75  IFIIAT.GT.GATEI  IRNG)  )GATE(  IRNGUINT 


1RNG=IRNG*1 
IF( IRAG.GT.GTPAXIL »)GC  TO  57 
GC  TC  75 

IF(aBS(XE(LIKE)-XB(L  INE) 1.GE..ep*1)GO  TO  126 

SUL  IAE)>999. 9999 

RANG  E«XE I L I NE )/ CCS  I ANGL  E ) 

GC  TC  127 

RANGED (YE(L INE) >SL(L  INE) *X9(L INF) )/( S INI  ANGLE )-SL (LINE )*COS (ANGLE  I 

*) 

RANGE=>ABS(RAnGE  I 
1RNG=PANGE/FLCAT ( IGTLENG)^23. 

IFUAZ.GT. IPFIB.ANC.JAZ.lt. IPHI  )GO  TO  IBB 
IF(IAPS(IRNG-GTKIn(L  n.LT.IABSI  IRNG-GTM4X(L)OGO  TC  129 
IRNG»GTMAX(L) 

GO  TC  128 


IRNG=GTMIN(L) 

1f( IRNG.GT.GTNAXIL I) IRNG«GTHAX(L) 
IF( InT.GT.GATEI IRNC)  )GATE( IRNG)>INT 
IF(ISV6T(L).EC.0)GC  TO  53 
JRNG>ISVGT(L) 

IFdAESI  IRNG-JRNG)  .LL.DGO  TO  53 
IF(IRNG-JRNG)54. 59.56 
JRNG» JRNG«1 

IF  ( INT.GT.GATEI  JRNO  )G  AT  E(  JRNG  )*INT 

GO  TC  55 

JRNG»JRNC— 1 

GO  TC  52 

ISVGT(U>IRNG 


ISVGT(U>IRNG 

IF(GATE( (RNG-1) .EQ.0  .and. GATE! IRMG« 1 ) .EQ.0 )GATE ( IRNG- 1 ) >I NT 


IF(  JAZ.LT.IPFDGC  TO  59 
CALL  GETLINI ILNE.L. IPOINT) 
LNSAV(L)>ILNE 
GO  TC  58 


I 


COM  UUE 
INf»lt 

00  \6tS  JE>ltHXNCE 
J>IEC(JE) 

L’cLECEM  J) 

IF(L.LE.»»GC  TC  18t5 
I£CH*JEC(lT 

IFUECF.EC.  ■'»GC  TO  1 e65 
IE«SLE^G^ECF  ) 

Og  1360  I>1  . IE 

IFIIECAZ  I lECF.n  .Nt.  lAZ  )G0  TO  1»5? 
KANG«IECRNG(lkF.l  I 
IF«KRNG.GT.223)STaF 
IFIIM.GT.GSTEIKKNC)  )GATE(KRNG)*IUT 
IFIGATECKRNG-n.LT  .IM  IGATECKMG- 1 >»l 
IFIKRNG,GE.219)GC  TO  186/ 
IFCGATECKRNG+n.LT  .INT  )GATEC<RNG+ 1 » 


1360  CONTINUE 
1865  CONTINUE 

CALL  TP3B  ITCH.  GATE,  22^f  $5'»,  $5.M  ) 

63  CONTINUE 

GO  TC  212'' 

500  MRlTE(6|1fl4) 

U4  FQRWATClF''.23FENC  CF  TAPE  CNCuUNTEREO//// » 

GO  TC  2li1 

501  0R1TEI6|1''5) 

105  F0BPAT<1f,1,28HERR0P  IN  WRITING  ENCOUNTEREO////I 
GO  TC  21or 

C PRECICT  PCSITICN  AT  GIVEN  TIME 

1900  WRITE  16, 1902) 

1902  FORHATC*  ECfC  NC/  FHPM* ) 

CPAA36P  CEPENCENT 

REACCEk’,  1737)JE,  JT  IM 

WRITEC6.6707) JE. JT IM 
L>LcGEN(JE) 

IF  CL  .LE.  0)GC  TC  27' 

GTIF«ECUR(JTIM) 

IFCGTIP.LT.FTIMECL  )l  GTIM»GI  IM*2A, 

PX»AXCL)AGT IP4BXIL1 
PY»AYCLI *GT IP^EYIL I 
OT*AeS(GTIP-FTIP£lL) ) 

RAC1»CT*SCRTISCE2I  L) /I  NOBIL  1*1  AXIL  1**24AYIL  )**2m 

RAC3>3.*RAD1 

RNG=SCRTCPXA*24PY**2 I 

AZPaFACIR(-PX,-PY) 

WRITE(6f  19C'5)AZHtRNG.RA0l.KA03 

1905  FORPATC*  AZP,RNG  »*2F6.1,'  RAD  1SC»*F5.I,*  RAO  ?SD»*F3,1| 

GO  TC  213<^ 

3090  CONTINUE 
270  WRITEl6,271l 
271  FOBPATl*  NC  ECHCM 
GCT  C 2103 

8.. .ACCEPT  PPI  FRCK  BENOIS  DISPLAY 

2000  HRITEI6t230 1) 

2001  FORMAT  1 1 F2, • are  A/ INTENSITY • ) 

REACC53.2005) lABEA.INTSW 

2005  FOWMATC2I5) 

WRIT EI6. 2005) lABEA.INTSW 
INTEN^SAINTSW 
C...REAO  PPI  DATA  FRCP  LX 
2010  CALL  CCNTURIINTEN,  lABEA, AREA,  lECINT,  IGL) 

REACC5O,2015) INTRPT 

2015  FORMAT(AA) 

WRITE(6|2316) INTRPT 

2016  FOBMATClX.AA) 

IECN*IECNC-1 
1F(IECN.EC.9)GC  TO  21B9 
00  202A  I»I.\TSW.7 
ICT*0 

00  2u25  IECF>1, lECN 
IFClECINTIIECHi.EQ.l )ICT>ICT4l 
2025  CONTINUE 

IFCICT.EC.eiGO  TC  202A 
WRITE(6,2'27) 

2027  format  (IJ-P) 

WRIT E(6| 201 2 nCT. STORM!  1 )*  UREA 

2012  FORMATllB.AAt*  ECHLES  FOUND  WITH  AREA  GREATER  THAN'.IA,*  SC  KM'/ 


)SC»'F5.I,'  RAO  ?SD»'F3,I| 
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^ARE4/AZIKUTH/R4NGE*) 
nc  2323  lECh=l, 1 ECK 
IF(IECINr«IECh).NE.I  )G0  TO  2.' 

X=CCFn,  1.IECM 
Y=CCFM,3,IECF) 

R=SCPt(X*X*Y*Y»*FLCAT(  IGU 
CALL  ARCT*N(X,Yt  Ai  ) 

1AR=AREA( IECF)*.3 
M=R>.E 
I A2» 

WUIT  £«6,  P”"!  » 1AM,  I ZZ  , IK 

2313  FURFATl tX, lA, I6i 16 1 
2323  CGMIKUE 
232A  CUMUUE 

GO  TL  2-^TA 

2339  WRITE(6i2?n)STCRMlIMSW  ),  lARFA 

2311  FQRFATl*  ’NCtAA,'  ^CFCES  FOUM'.  WlTF  AREA  GREATER  THAN*, 16, • SC  KM* 

2316  IF(  IATRPT.EC.6ESTOF1 CO  TO 
GC  TC  2n  13 

C PRECICT  WFEN  ECHC  ULL  3E  NEAREST  GIVEN  POINT, 

3333  WRITE  (6,  33'*  1 ) 

3331  FORPATl*  ECFC  NC/  tlV/  RNG'l 
C*6*36'J  CEFEKCEAT 

REAC  (fkl,3P’t?2>JE,  AZF,  RNG 

3332  FOWPATn5,2F5.;') 

WRITE(6,3*03)JE, AZ', PNG 

3, S3  FORPATnd,2F5.P) 

C**«««* 

L=LEGLM  JE) 

IFIL  .Lt.  *)GC  TC 
AZ=AZF*OTR 
XA=RKG*CCSI AZ ) 

YA=RNG6S  IM  AZ  ) 

0EN=AX(L)**2+AY(LI«*7 

PTIP£=(AX(L»A1XA-BX1 L» )+AY(L )*<YA-BY1L) ) l/OEN 
DT=AeSlPT IKE-FT IPEIL I) 

SO=STE2( L )/ACB( L ) 

IFUCeiL)  ,LT.  3)St=  .5 
ERTP=CT*SC 

T ’ = JFFFY(PT  IKE-ERTH 
T1=JFFFP«PT IRE) 

T2=  JFFRMl  PT  IRE-f  ERTM 
PX=AX  (L) *PT IRE+eXl I ) 

PY=AY  (L) ♦PT IRE+EYl L ) 

DS=SCRT( (PX-XA)**2+( pY-YA)**2» 


SDC=STE2«L)/(NCe(L )*CS) 

WRIT  EI6,3CS5)CS, SCL, T1.T7, T 
FORRAT(lx,F5.  1,  * MlL^sU/-' 


GO  TC  2133 


T2 

•FA.l,' > 


I '16,' 


•14, 'I') 


2JP  WRITE(6,231) 

201  FORRATl*  IN  ERR' ) 

GGTC  210'* 

C..,CUALITY  CCNTRCL  FCR  EENCIX 
5300  WR1TE(6,5305) 

5005  FORRATl*  PPI  CFECK') 

REAC(60,5P13) ALFA 
5010  F0RRATIA6) 

WRITE(6,5R1P) ALFA 
WRITEI6,E'*15) 

5015  FORRATl*  TEST  PATTERN') 

REACI5S,5P1?)BETA 
WRITE  (6,5ei3)BETA 
IFIALFA.NF.AF  YESItO  TO  51’.5 
11*1 

5020  CALL  TFREACnC,GATf,  220,LENR,»5I  3P,»*4SCI 
DU  3fc25  K=1.19 
5025  IHSK (K)>GATE|K) 

CALL  CECCCEI  IRSK,  I AZ  , HLT,  ISTC,  JUL  , I T IMb , lOL Y,  IGL  , I TC  ) 

ISAVTR*ITIRE 

JULSAV>JUL 


IFIII.GT.DGC  TC  5.3  5 


WRITE(6.5P3P)IAZ, ITLT. ISTCtJUL, ITIHE,  lOLY.IGL.ITC 
5030  FORRATl*  lAZ  TILT  STC  JUL  TIME  CLY  GL  TC • /2I 5 ,2 14, 1 7, 14 ,21 3) 


C.,..SUM  UP  eiTs. 

DO  5031  1*1, 7- 
5031  IKNT(1)»'7 
5035  DO  5SA5  K»23,21S 
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)GU  TO  SC5* 

LAST  racial  =•, 


I A) 


K\T=GATE (K)/e 
IKM  IKNT  )=IKNT(  KNT  )♦  1 
53A5  COMIAUE 

IF( 1 AZ.N6.3S9.ANC. n .Nt.36 
WRirE(6,S-’5'')IAZ,I  I 
5050  FORNATI^  •'LAST  AZI»1LTF  =',U, 

WRITE(6,5."'6'') 

5«6T  FQRRAkJ’INTENSITY  BIT  COUNT  • /9X,  IM  , 9X,  1H  2 , 9X  , 1H3  t9)i , 1HA,9  X , TH5  , 
*9X,1F6,9X,1F7 ) 

WRITE  (6.5T6  5M  IKNT  ( n,  1=  1,  71 
5365  FOPRATltllT) 

GL  TC  513C 

5355  CALL  T IMCFF (JUL , IT  PI  E, JULSAV,  ISAVTM,  IDFF , I IFF ) 

IF( I AES(  ITFF) .GT .1  .OR. lABSI IDFFI.GT. 1 )WR1 T£( 6,5 ’7Z» JLL.I  TIRE 


5373 

5375 


5030 

5385 

539? 

5395 

5130 

C.  • . . 

5105 


FQRRAT (•  3ERBCB  IN  CATE/TIMt,  CATc  =',IA, 
IF(  I AZ.NF. I I- nWRI TE(6, 5375)  IttZ, I I 
FURRATI' AZIRUT)-  =',I4,»  RAiJIAI.  =»,IA) 
11=1  l + l 
GO  TC  5)2" 

WKITE(6,5C85) 

FOPR  AT  (*:»EQF  ENCCUNTCREO  ON  I/P'////) 

GO  TC  21?'’ 

WRITMeiS^gS) 

FOKRATI' ZERRCB  IN  t-EiCING  1/P'////) 

GO  TC  5Z2? 

IFIBEIA.NE.AF  YES)C0  to  21»;' 

SET  intensity  SWirCFES 
WRITE (6,5105) 

FOBRATI* ’SET  INTENSITY  SWITCHES  TO  1 2 3 
IP=1 

11=1,351,). 


' time  =',m 


1 2 3 2') 


't  )/!■ 
H ) 


5126 

5110 

5112 

5125 

5130 

5115 

5116 

llil 


DO  5 n?  _ _ 

00  yl?5  IJ=1,10 
IAZ=ll-'IJ-2 
GATE  (1)=  IAZ/lpi» 

GATE  (2)=RCC(IAZ,  1 
GATEI^)=RCC(IAZ, 

DO  5 126  1K  = A, 19 
GATE  (IK)  =.' 

DC  3 11k?  1=2'',  159 
GATE(I)=(I*B)/2‘3 
00  5112  1=160,219 
GATE(I)= IP*« 

CALL  TFWRIT 1 1 1, GAT  b, 223, *5 1 15.$512£ ) 

CONT  INUE 
IP=IP4l 

IF( IP. EC. 5)  IP=1 
CONT INUE 
GO  TC  2W 
WRIT  E (6,5116) 

FOpRAT(‘0ENC  CF  C/F  TAPE'////) 

GO  TC  21'^ 

WR1TE(6,5121) 

FOR-VATI*  0ERPCR  IN  WRITING  O'!  n/P  TAPE'////) 

GO  TC  21 
ENC 

FUNCTICn  LEGEN(K) 

CORRCN  XP(900  ),XE(  ^7*'),YB(  9aC  ),YE(2C?),SL(  2C3)  ,PHId(20  ')  ,PHIE(2<!i5l 

♦ , ISVGT(4’),GTRIN(A.  ),GTMAX(',7  l,SLENG(  It  ) ,NHAM  ( K I , I EC  A Z ( 1 0 , 900)  , 

♦ lECKNG(10,9i3 1, I6CN0,IGTLEnG,MAP( 1 20  1 ,MEN,CQF I 1 Z ,h ,U 1 ,G ATE (220) , 
♦IGC.eCXTR,  IHSk(19) 

INTEGEPA2  lECAZ 
CHAR  lECPNG 
CHAR  GATE,IFSK 
LEGEN=0 

IF(K  .GT.  0 .ANC.  ^ .LE.  M C.M  )l  EGEN=MAP  ( K ) 

RETURN 

ENC 

SUeRCUTINE  CCNTUP(  IN  TEN,  I AREA. AREA, lEC INT, 16L ) 

CORRCN  X8(930  ),XE(  2 P J 1 ,Y B(  900  I , YE(  3i( 0 1 ,SL ( 20<1 1 , PHI B(  20 1 .PHI E (200 1 

♦ , ISVGT(43  ),GTMIMA5),GTMAXI  Ai'  ),SLENG(  1 P ) ,NhAM(  1 0 ),  1 EC  AZ  ( 1 0,903)  , 
*IECRNGn8,9?0),  1ECN0,IGTLENG,MAP(  100  ) .MEN.COF  ( 1 0 ,,,  ,1  0 1 .GATE  (220  I , 
♦IGC,eCXTR,IFSK( 19) 

OIRENSICN  AREA( 13) 

DIRENSICN  IECINT( 1 . 1 
INTEGER  PCXTR 
INTEGER*2  lECAZ 
CHAR  lECRNG 
CHAR  IARY(36<7,2£0  I 
00  13  1=1, 1j 


, GATE,  IFSK 


00  IC  ((*1,4 
00  1 ».  J=  1 , 1 
\i  CCfU.K,  J)  = (’. 

4«tA^X=5’*IAREA 

ARtA*(-. 

5 I»1 

lcCNC=1 

I0L*1 

jeGATfc=l99 

S CALL  TPREACd  v'.GAT  t,  r2c,LLNR,  45",$5?l) 

CALL  TPWRITm,GATc,52-,ii.2,*5i?A^ 

DC  2..’  K=1,19 
2jJ  1HSK|K)=GATE(K) 

'STC.JUL,  ITIME,IULY,IGL,ITC» 

BOXTR^IT  IRE/lr®' 

JTIR=PCXTV 

WR1T£«6,2^4» 

2i34  FURPATI* -'CAY/  TIKE/  TILT*) 

, WRITE  (6,275 )JJL,JT IM.fTLT 

235  FCRNAT(Ia,2I6) 

233  IF(I.FC.1.ANC.IA2,,vF.7>GO  TU  5 
IF(IAZ..YF.I-1)GC  tc  u 
, „ 00  2-72  J=1,IGC 

232  lARY(I,J)=fl 
1ST*  IGC>2Y 

00  201  J=IST,213 

„ IFiGATfcl  Jl.LT.lNTtMCATElJ)*'’ 

231  I ARYn,J-19)  = GATE(  J) 

1 ARY  (1,2  3-')  =? 

IF(IA2.EC.359)GC  Tl  !5 
1*1+1 

. GO  TC  5 

16  WRITE(6,17) 

17  FCRRATI* '•A**  - IWiLIC  OAfA.  AZIMLTF  EKOQR  - »♦*  • ) 

. RtTUP^ 

13  IF( I .AE.362 1GC  TC  16 
J=JECA7E 
ISAV=>6») 

C...ZERO  CUT  lECAZ  ANC  ILCHNG 

00  9iS  IK=1,6»3 
leCRNGd  ,IK)*a 
It'CRAG(2,IK)  = ’ 
lECAzn,  IK)=3 

935  IECAZI2,  IK)=<1 

1 = 1 

23  K=3 

AR£A« IECKCl*e . 

MIaTEN*56 

C...Lci4TE  NCRTHEAST  CCRNFR  OF  ECHO 

25  J=J-1 
JSAV=J 

IFU.GT.  IGOGC  tc  26 

I*i+1 

J*JEGATE 

IF(I  .EC. ISAV1KETURN 
GO  TC  25 

26  J1*J 
11  = 1 
JHI\»J 
JHmX* J 
1?<IK*I 
I(«AX*I 

CALL  eNCRVl  lARV,  I,  J,  M,Jl,tl2/.) 

GO  TC  25 

C...SCAN  FCP  NEXT  PCINT 


rLK  mcxi  rLiiNi 

fi.  99  (155,35,65,5:,  <.5, 75, c5, 951,  IP 

155  11*1 


J1*J+1 
GO  TC  13  1 
11*1-1 

IFII  KEC.")  M*36J 

J1-J+I 

GO  TC  13  1 

J1*J 

n«i-i 

IFIIl.^^^ri 11*362 
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•. . 


11  = 1-1 

IF(I  1 .EC.  ’)  n=3t  ■ 

GC  T C l <i  1 

n=i 

GO  TC  1-11 
11=1+1 

im  i.Ec.?6nii  = i 
ji=j-i 
G^:  TC  in 
ji  = j 
11=1+1 

IFU  1. EC. •’61)  11  = 1 
GO  TC  m 
J1=J+1 
11=1+1 

IF(  I 1 .EC.'>61  ) 11=  1 

call  "^ORY(  I4RY,  I, J,  I1,J  1,S12^I 

IP=IF+1 

ICCU6T=I CCUNT+1 
IFdF.GT.Pl  IP  = 1 
IF( ICCUNT.GT.eiGC  TO  125 
GO  TC  29 
K=K+  1 

IF(K.FC.9'»’I  ISTCP  CCNfUK 
J=J1 
1 = 11 

IFU -IMN.GT.  lenct  TC  122 
IF(  I .IT.  I»'IM  IM  Ii\=  1 
irU.GT.  I^AX  1 IYAX=  1 
GO  TC  123 
IH=I-?tJ 

IFUB.LT.IMiMIf'lN  = lF 
IFU.GT.  Jl'AX)JMax=J 
IFU.LT.  J»'IMJMIi\=J 
IECAZ(2.K)= I 
IECR^G(5,K)sJ 
ICCA2n,K)=I-1 
I tCR6Gn  ,K)=J 

IFU  ARY(  1,  J 1 .LT  .P  INT  rNl!Ul+TEf4  = lAAYl  1 ,J) 

I Ai^Y  1 1,J  1 = 1 ABY(  I , J 1+  ia>. 

ICCOM  = l 

IFUF-3171,71,72 

IP=I F+6 

GO  TC  29 

IP=IF-3 

GO  TC  29 

CCNT INUE 

DO  L=1,K 

I4=IA“SnECAZ(2,K)-ieCAZ(2,L)  1 
IFUA.LE.l.CR.n.GE.;59)GO  TO  165 
GO  TC  139 

IFU  AESl  lECPNGI  2.K  )- lECR.NGl  2t  L 1 l.Lb,  1 IGO  TO  17? 
l4RYUECaZ(2,L)  , lECR-MGl  2,L  11=  lARYl  lECAZl  2,l1  tIECRNG(2  ,L1 1-lP? 

COi\T  lAUE 

WRITE  (6,  151  1 lUN,  I^AX,JMn,JMAX 

FORMAT  P^’ECfC  0ETWfEN*,2IA,  • AZ  . ANC»,2IA,  • RANGE  NOT  CCNTCURABLE* 

♦ 1 

RETURN 
K=K+  1 

IFCAZM(K)=IECAZ(1,L  1 
I£LRNGU,K)  = IECRNG11,L) 

IFUMN.EC.  1)JMIN=2 
MAXI NT=1NT 
IHCLC=IMIN 
MINTEN=KINTEN+A 
I=IKIN 

IFUPIN.LT,1)I=36»'  + UI\ 

DO  223  J=JHIN,JMAX 

IFU  ARYI  ItJ-1  )>GT.  ISIC.ANC.IARYI  I.  Jl.GT.C.AND.IAKYIl  f JI.LT.6A) 

♦ lARY  U,J)  = I ARYl  l,J)+l-»3 
I M*  I - 1 

IFUP.EQ.'’)  IM  = 36i 

IFU  ARY(  IM.Jl.GT.l  3.Ai^U.UkY(  I,J).GT.0.AND.IAKYU,J).LT.1P?) 

• 1 ARYU,J)  = IARY(  I,J)+1?J 
IFUARYI  I.Jl.LT.K  )GC  TO  221 

IFU  I ARY  II,  Jl-l  ’ )1  .GT.MAXI.NT  1 MAXI  NT  = I AR  Y (I  , J 1- 1 . k.’ 

CONTINUE 
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DC  221  J = JH  ji*/SX 

1F(  l4fJY(  I ,J  ).GT  . I"'.  I 'RcA(  Ii:C><U»=AReA«  ICCMO  )*FLOAT«  K;L*IGL*(2*J-1)  » 
♦*3.7266<»6E-i 
CUAT lAUt 

Irdf'IN.Lc.  tf'flX  )CO  TC'HV 
I r(AHCi(  IFCaC ).LT.  iRfciNX >GU  TO  191 
l'^lN=t1-CLn 
I = If'IA 

IF(  I tN,  LT.  1 ) l=3t . -UMij 
OU  22:  J=JK1N,JKAX 
IFUAtiVI  1,J).LT.P  )rc  TO  223 
I ARY  ( I,J  )=I  ARY(  I , J )-  '>!  . 

IFIIAwYI  I,J).LE.KIivTfc.Y)IflRY(  = 

CCNT  lAUE 
lMli\=IMA*1 

IF(  I MA.  LF.  l^AX  ICO  TC  2?2 
I=1FCLC 

1FNFCLC.lt  .1)1=1 
J=  JSAV 
GO  TC  195 
IWIA=1FCLC 
1 = IMN 

IFN  .LT.  1 ) 1=362'*  IM  IN 

00  19*  J = JK1N,JI'AX 

IFN  6RYI  I,J  I.GT  . r I lAKYI  1,J  ) = ’ 

COAT lAUt 
I.'11A=1MN'*1 

IFN  MA,  LF.  I^'AX  ICO  TC  192 
I=1HCLC 

ifnfclc.lt  . 1 n = I 

J=JSAV 

IFI AREA!  leCAC  I.LT.rLCATI  ItrtcA  ) ICO  TO  23 
IFIK-L.LT. A IGC  TC7' 

CALL  CSTLACIL.K) 

IFCIATNECAC)=WIATLN/E 
ieCAC=IECAC+1 
IFNECAC.lt. 11)CC  10  2i 
IFN  .LT.36C  IWRITEI  -f  19o) 

FOHRATI^TEN  ECHCtS  FCUin  before  PPI  CCMPLtlED,  RLStT  tCHO  LIPITSM 
RtTURA 

WR1TE(6,5?6 I 

FOHFATMTECF  ENCCUATLREO  UN  1_V///I 
STCP  CCATUR 
WRITF(6,5''7  I 

FORPAN*  ■'ERRCR  I(\  rF  ClNG  ON  10*////| 

RE  fURA 

kRITt (6,5’3  I 

FORRATiIf" , 23FEAC  LF  TAPE  ENCOUNTERED////) 

STCP  CCATUR 
WRN  C(6,5''5  ) 

FURRATI1F1,26FERR0R  in  writing  F"/C0UNTERE0////I 
ftClORA 
£ NC 

SUHRCUTl AE  3E1L  lAI LINE  ) 

CCRRLA  XPI9i  2 ).x  t(  2'’  ) , Y B ( 9^  ^ ) , YE  ( 2 ' 3 ) , SL  ( 2 ’ 2 I . PH  I CN  2 2 :)  .PHI  E 1 2 ) 

♦ , ISVGT  (Ai  ),GTM1  AU.  ) ,GTMAX(*.J  I.SL  ENG(  1 f ) ,NF  AMM  : I , lEC  A 2 I 1 0,923)  , 
*!ECRAGN  ■’,9.3  ),  I tC\0,  IGTLcNG.HAPI  IgP  I.MEN.COFI  1 2 , , 1 2 ) .GATE  (220)  , 
♦IGC, ELXTR, IFSK{ 19I 
CHAR  lECRAG 
lNTLCtPA2  lECAZ 
CHAR  GATE.IFSK 
A=X6(LINEI 
0=YB(LlAE) 

C = XE  (LINE  ) 

0=YE  (LINE ) 

SL(LIAe)  = (0-B)/IC-N 
Call  arctania.b.e) 

CALL  APCTAMC.C.F) 

IF! A.EC.  1..AAC.P.EC.  .)E=F 
IFIC.EC. AAC.C.CC.  .)F  = £ 

IFIE.EC.FIGC  TC  5 
IFIE-F.GE.179.5)G0  TC  2 
IFIF-E.GE. 179.5 )G0  TC  4 
GO  TC  5 
XC(lIAE)=3. 

YE(LIAE)=e-SLIL IACI*A 
PH19 (LINE )=£ 

PHIE(LINE)=E6'i. 


2'* 


25 


.jpHieii  r^.i=E 

• ’ ) 

■ 1 ) 


IF(Yt(UKE)  .LT.  1 
LlNfc  = U.»r4l 
X M L I^E»  ='^ 

Y>^(LI^E)»VE(L^\f• 

SLILINE)=SULINC- 
PHl:5(Lir,£»= 

PHl£^lI^£l=F 
I'=(YC(LINE)  .LT.  } 

XE«LUE»=C 
YEIL  INF)=r 
fti:rUR^ 

XH(LUE)='’. 

YUlI.'E)  = R-SL(UIi\E»*J 
XE(LINE»  = /' 

Yc<  L INE) =e 
PH10(LINE»  = . . 

PHlt (LINE )=F 

IFIYEdINE)  .LT.  1 .)FH  IR(L  I.IC)  = E 

LINE=LINE+1 

X UL  UE»=C 

YlMLINE)  = C 

Xf  IL  1NE)  = ’<. 

Yt(L  INEI=Ye(L INE-1  ) 

P-(IP(LII\^l=f 
PNIE  tU\E)=?fc'. 

1F(YF(LINF).LT. 1.)FHIE(L IN£)=F 
SL(LdE)=SL(LlNc-l  ) 

R^TUPN 

IF(E.LT.F)GC  TC  2 

PHia(LINE»=F 

PHIE(LIN£)=t 

RETURN 

PHIU(LIi>iE»=C 

PHIE (IINE)=F 

RETURN 

E\L 

SUtRCLTINE  iRCT /N(L, V, (HETa) 

RTL=57. 29577951 

IFIU.FC.  ■ ..«NC,V  .G5.  . )TE6U=.  . 

IKU.EC.'  ..iNC.V.LI.  . ITEETAs  137, 

IFIU.GT,  .)TEcT^=9  ,-ATAM(  V/U  I’CkTD 
If  (U.LT.3.1TEET  A = 21''  .-MT6'i(V/U)*fiTC 
RETURN 
ENG 

SUdRCUTINF  ISCRT (V  ,KY,id 
REAL  V(2» 

INTEGER  KY(2) 

iMTAILIZE  key  and  PulNTERS 
DO  2'  I»1,N 
KY(I  1 = 1 
K-»l 
N =N-  I 
NK1  = N ’ 

TEST  FCR  ENC  CF  SORT 
IF(NK1  .LE.  "IRETURN 

NCT  FIMSfEC.  STAitT  NEXT  PASS. 

K1  = NK  1 


>K1 


NKI  = !. 

00  J=K 
M1=KY(J) 

M2  = KY(J-»n 
IFIVIFl)  .L£, 

KV(J)=R2 
KY(J4l)=M 
IF1NK1  ,FC.  ? 

NKl  = J-»1 
CONT INUt 
K.;  =R  AX^H  1 ,K'  ) 

NKI  = MN3  (NP.NKI  ) 

GO  TC  25 
ENL 

SUERCU1 INC  ECHC( lEtH  » 
CORPCN  xeC9k..1),XE(2p-' 

♦ , ISVGTIRO  1.GTPIM4.  ), 

♦ lECRNGII  ’,9'’f’  ),  lECNO, 
♦IGCjECXTP.IFSKI 19) 

1 irEGEP*2  lECAZ 
char  iecrng 

CH«R  GATE.lfSK 


V(R2  ))tC  TO 


' )K/=J-1 


- ”)fYB(9.;,''),Yi.(2*JC),SL(^"2),PHltt(2Ki  1,PHIE(2t!4)l 

♦ .ISVGr(R0),GTPIM4.  I.GTMAXI  A"*  ),SLENG(  1C),NfAMn:),IECAZ(la,9i  J)  , 

■ 0,IGTLFnG,MAP(  lk'),MfcN,C0F(H,«,,l?),GATE(22P), 
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•»^5 


197 

19? 

199 

2i2 

2^3 

9^\ 

993 

9i;3 

2oi> 


Pl  = : .H13926 
Iii=SLENG(  leC*-  ) 

= I CCM 

DO  9 1=1. IE 

Y = CCF(1,  >,ICCH 
X = CCF  ( 1,  I , I FCF) 

ou  \=?,  iF4y 

U=2.♦t^-l)♦PI*I/FH.4IC  IE) 

CCSU=CCi «UI 
S1.\U  = SIN<IJ) 

X=X+CCF( 4, 1 , IFCF  )*f03U+CUF(  2,  I ECF  ) * S IN U 

Y = Y + CCF(  \ ,3,  lECh  l^CO'iU^COFI..,  4,  IHCF  )*Sn.U 
CONTINUE 

I tC9N''(I  ECH,I  )=SC9T(  X^X^Y'CYI/FLOiTI  IGILENG  )-»2 
C4LL  JRCTiMX.Y,  /!£  ) 

IEC4Z ( lECF, I )=AZ 
CCNT  INUc 
REIUFN 
E )L 

MXNCLN) 

W P tf  n /t  Ta  ^ 


YEI  2!^'^'),SLl  2?7)  ,PHiai2.:  ’)  ,PHIE  (2!-;i‘) 
SLENCI  1C  l.NHAMIl  1 , I EC  A Z ( 1 2 , 

P(  l;  " ),HCN,COF(  lv,<.,1i)  .GATE  (229)  , 


♦ IGC,bCXTf',ll-SKt  19) 

OIPENSICN  IFCINT(2  7 ), PH  IFULUI 2.  )) 
lNT£GtS*2  lECaZ 

CHAR  lECRNG 
CHAR  CATE.n-'iK 
0(J  19’’  J=1,22''' 

PHIHCLC(J)=  . 

I PLl M I J ) = . 

0(’  ?:■  J=l,f'XNCLN 

lCcuNr= j 

UJ  2.  I = l,PXNCLi\ 

IFlPMECtn  I).EC.9  9.)G0  TU  2*2 

ICCUM=ICCUM4l 

IF(  ICCUNT.EC.  DGC  TO  19E 

IFlPFlei  1)>PHM\)  )9  -♦  199,;  :2 

1PL1M(J)  = I 

PHlFIN=PFie(  I ) 

PH£NC  = F)-  IE(  I ) 

GU  TC  2 <2 

IFtPH^NC-PHlEII))?  2,2-2,193 
CUNT  INUE 
L INt= IFC INT (J  ) 

PHIhCLCl  L INE)=99'». 

CONT INUt 

IF(PXNCLN.GE.23- ISTOP  I.NSOPT2 
00  9t  I=1,36G 
I NZ= 1-1 
L'1CM=f 

00  9« 1 J= 1, WXNCLN 
K=IPC1NT (J) 

IFIPMEI  K ).L£  .1  az.  i.NC.PHI  E(K)  .GF.  lAZ  )LNCN T =LNCN T+  1 
CONT INUE 

IF( I .tC.  1 )GC  TC  9 
IFUNCM.LE.LNMAXItO  Tu  9.:  ;i 
L9I*AX  = LNCNT 

IP(LN‘'AX.GT  .A'’)SIDP  LNSOi^IA 
COnT  I NUE 

IF  I I LNNAX  + MXNCLN  ).  t T .2.,;  )STOP  LNSGRT- 

00  l=i,lnpax 

IPCIM  (L  + pxNCLM=L-**4XNw»LN 
xai  L + f-XNCLN  ) = 999, 

m<ncln=mxncln+lncni 

RETURN 

e ac 

S'JBPCUTINF  GETL  IN(  ILNE.L.IPOINT) 

COFf'CN  xe(v'3),  XE(  2<-,  ),YB(  5»fi  ) , Y * ( 22  C ) , SL  ( 2 ” ) ,PMl  >3  ( 2 v'  1 , PH  IE  (2 
7.  ISVGT(A  ’ ),GTMIN(9  ) , GTM  AX  ( 9 J ) , SL  ENG(  l-’J  ..XHAM  ( 1 1. ) , I EC  A Z ( 1 - ,9  • M 

• fecRNcn  ),9:u I,  I ecno.igtleng, mapi  irc i.men.cof ii<.,a,t“ 

♦IGC.etXT^.IFSKI 19) 

OIKENSICN  IPCINT(2;?) 

1 ITCGEB  BCXTP 
INTEGER  GTHIN.GTFAX 
INTEGcR*2  lECAZ 
CHAR  lECRNG 
CHAR  GATE.IFSK 
ISYGr(L)=* 


lECAZd-  ,9  • M , 
,1«  ) ,GATC(22'’)  , 


110 


I1WW-* 


ILN6=ILNE^1 

IFIILNE.GT.24  •'’)ST0‘  t£Tl  IN 
J=IPCIKT IILNE) 

IFIxeU)  .cC.999.  IRtTtRN 
A=Xt  Ul 
B=VL>  I J) 

C=XE  U) 

0=YE ( J) 

RI^SCRTI «*4+E*e ) 

R>=SCPTI C*C+C*C ) 

R->=2P‘  .*FLO«T  I IGTL  LNG) 

1FIR1  ,GT  .R3.4NC.B2.GI.R3  IGO  T(j  57 
IF(R1-R2J76,76,  7-i 
r,TMML)  = R1/FLCm  lOUENGI+2:. 
GTRAXID^RP/FLCm  IGTL  tNG  » *2  ' . 

GC  TC  75 

GTKIML)=R2/FLC«  n IGTLcNG  )^2<’-. 
GTFNX(L»aR1/FLC«T(  IGTLtMG»-»2«  • 
IFIGTNAxm  .GT.EJ’  IGTMAXIL 
RETURN 
E'JC 

SUdRCUTINE  TINCFM  IC  I,  IT  I , IC2-  IT 2, 
IUFF=‘ 

ITFF  = 

I OAY1=ICl 

IOAY2=IC2 

IHR1  = IT1/1R;,?2 

IM1M  = (IT1-IFRI*  J’  :•  I/l". 

ISEC  I = IT  I-I  IFRl  ♦ i"'  r ♦ IM  IN  1*  1 » .' ) 
IUR2=  IT2/1I1  ?'i 
lNiN2  = IIT2-IFR2*r  7 
ISEC2*IT2-(  IFR2*1-  ♦ MINE*  I ► 5 ) 

IF(IC1-IC2) 1,A,5 


TINCFM  ICl,  ITl,  IC2-  IT2,  lOFF,  I TFF  I 


I/l".  . 

♦ lNtIN1*1»  .') 

)/^r^ 


IF(IC1-IC2) 1.A.5 
IF|ITI.EC.IT2)RETURN 
IF(IT1.GT.IT2)GC  TC  i 
IFIISEC2.GE.ISEC  mo  TU  2 
ISEC2=1SEC2+6(J 
IMU2=IMIN2-1 
ISCCF*ISEC2-ISeC 1 
IFIIMN2.GE.IKINDCO  TU  J 
IKIN2*IRIN:2+6.; 

IHR2«IFR2-1 
IM.\CF  = IRIN2-IRIM 
IFIIFB2.GE. IFRl )GO  T l 9 
IHRi;aI»-H2  + 29 
I JAY2=IC4Y2-I 
MRGF*IHR2-IFR1 
IuFF=  IC4Y?- ICAY 1 

ITFF»-(ISCCF+liii:  *I  IMNCF+K  ♦IFRDFII 
RETURN 

IFIISECl .GE.ISEC2ICO  TU  6 

I>£Cl=ISECl4ft3 

IMIM  = IRIM-1 

ISCCF-ISEC1-ISEC2 

IFIIRIM.GE.II'INEICO  TO  7 

IHIMaIHINU6? 

IHRl  = lt-«1-I 

If'NCF=IMINl-IMIN2 

IF(IF«I.GE.IFR2)GO  TC  a 

IHH1«IFR  W2A 

lOAYUIOAYl-l 

lH«Cf*IHR1-n-H2 

I0FF»ICAY1-IC4Y2 

ITFF«ISCCF^  IMN0F4K'5*IHRI  F» 

RETURN 

ENC 

SUbRCtTINE  8NCRYI  I JRY.  I,  J,  IlfH#** 

CURRCN  Xei9i'?),XEI  , YPl  9C  3 ) , YE  I 2 3?  J t SL  I 2<»?I  . PHI  hl2  7E  ) ,PHI  E 12  12  I 

♦ ,I5VGTI-.'»l,GTFIM4  ) ,GTH  AXU.7  ).  SL  cNG  I IJJ.NHAMI  1 I)  , IECA21 12.9121  , 
*IECRNGn'*t9lp  ).  lECNO.IGTLENG.HAPI  12 1 ) .HEN.COF  (1  ' ,9 .1  PI  .GATEIPEtf)  . 
♦iGL.eCXTR.IFSKI 19) 

InTEGER*2  IECAZ 
CHAR  lECRNG 

CHAR  lARYI36i^,2?>.  I .GATEilFSK 
IFlIARYIlIt Jl I.GT. 1?;|RETUKN 


IFIIAPYI  11.  Jtl.LE.M  RETURN 

IFIIAPYI  I.J1).EC.'*.ANC.IARYI  I 1.  J ) .EC  .URETtRN 
IFU1.GT.2R?  ISTtP  tNCRY 


1.  Jtl.LE.MPETURN 
.J1).EC.'*.ANC.TaryI  1 


K’«l  1-1 

IFU2.fes;.'l  12=3c.' 
irin. ft. ■•61)13=  I 

i^«|45yi  II,  ji-fT ) 

in«uuvn7,  ji ) 

IG*UfiVl  11,J1-1  ) 

M»I«FV(  P.Jll 
IH  If  .EC.MTcnT-  ICi\T*l 
IFIIC.EC.t  ) ICNT  = ICM«1 
IF(U.tt.'')ICNT*IC\T-*l 
1F(  1 l.tt.'^l  ICNT=  1C  »T  ♦! 

IFIICM.CC.'  )RETUR^ 

IC=ia5V« I2,J1*1 ) 

I£=I 4PVI 12, Jl-1  ) 

n,Jl-l) 

IJ=  I HTVl  (3,  jl  + U 

IF! I e.6C.-, 4NC. It.^E  .? )G0  lu  2 
IM  IC.Nt,  •.  JNC.  IC.  \E  .?  )GC  fiJ  2 
^F^L.^t.6.4^l..It..\F.^l)^o  ru  ✓ 

IFdC.KE,  ‘.aNr.lG.f.E.'loO  TO  2 
IF(IO.NE.’.*»NC.lF.i\E."  )GU  TO  2 
IFllF.^C.*’#^NC.II.tNE./)GO  TO  ✓ 

If m,^E.'’,A^r. IJ.^E.^Muo  tj  > 
lFU4.^E.••,a^C.le.^E.^)Go  To  2 
I A3t  ( 11, J1  )s" 

r:iurn 

IFjll.£C.''.iNC.IL.(-0.«'.lND.IC.>.T.EC.21G0  TO  1 
''^£^J'^*‘^Cr**’'^'**C''iT.r0.2IGO  TO  1 
I F (I  C^T  , LT  IBETURN  ' 

I-RY(Il,Jl)*r 

m-ruR6 

E a 

SOHRtUTlKE  CSTLNtlL.I'  ) 

CCRfC6  XE(9.  llfXtd  ' )fVC19.r),YLl2i!  ’ItSLli'*^),  PHI  C 1 2 '’'•).  PHI  F 12  ‘"G  1 
• ,IiVaT(4’«l,GTKIM4  ) ,GI  •UX(*.r  list  ENG  1 1 l.NHAMl  I > ) , IFUZ  I'.t  - ')  ; 

MGC  ECxTplfhSM  1*  f »»HFN,COf}  1 ,4  , H ) ,GA 122.  f , 

isrdGLP^z’iccai 

char  IEC'<A6 
CHAR  GATE.IFSK 
OTH*.  1VA532 
OF»l 

P|*5.  U1^9?6 
lECH.  lECNC 
MCCUM»K-L'»  ' 

IHAP=9 

I f 1 PCCUM  .CT  . 18  1 1H.'M=PC'0UN1  /2*1 
NM»M  ItCF)»  IF  AM 
OG  1 J* 1,1^ 

00  1 I’^l.A 

COFl  J.I,  leCF)=k., 

A=IECAZ( l.L  ) 

R*IECRFG(1.U 

Xsi»P*SIMi*CTR) 

YS1»R*CCS1 A*CTH  1 

XLEKG=‘ . 

L1*L4  I 
PO  6 I»L  1,K 
A»IEOZl  1,1  ) 

R*IECPNG« 1, I) 

X8II  )=RASIMA*CTR| 

YBU  I»«*CCS1A*CTP) 

xs«xeii) 

YS-YHII) 

XLENG=XL£NG-»SCRT  M XS-XSl  )**2+l  YS-VS  1 )**2) 

XS 1 > X3 
YS1»YS 
CCAT lAUE 

SLENGIIECHUXLEAC 
DO  2J  NF»2, IFAM 

Sl«i  . 

OG  li'  I»L1,K 

xs-xem 


PX»IKS-X«>1)/IS2-sT  ) 
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pv.|ts-rsn/«S2-si ) 

QX«XS-PX*S? 

QV»YS-FV*S2 
C'4«2*N*PI/XLENG 
U1«C^•SI 
U?«C(V*S2 

sikT«siniui ) 

icsi'»*'  — 

T|*CJ/£N!S!.N?*PXM  CN**|-?!f*COS 

C CP ( I » I 
CO^CCEPF  CF 
T’«-CX/CN* 

T »«-CX/CX* 


IN2*SIN(U2I 
CSJ»CCS(UTl 
0S?»CCS(y2i 

L CF  ■v»'  A vwunud 

,CN*SIN2*PX*(  CN**I  -2)*C0S2'»S2/CNPSIN2  I 

CN*siNUPx*(CNMU2i*cest4sl/CN*snil 

j’ >ECH|42,/XLENG*It1-, 
* CCOKLS 

*CCS24PX*«C\**T-2»*Sli'#2-S5/Ch*CQS?l 


CC 


LUi;i^^.3,][CP>«CCF(NF,  i,  ItCH 
Cb  ClEFF  CF  SJNE  FOH  Y cioRtS 
T2«-CY/CK*CCS 24PV*ICN*<»(-2I«'$  IN2- S2/CK*C0S2  ) 
FI*~CY/CR*CCSl4PY*(C*.**l-2l*Slr<l~  s1/CN*CtJS  t ) 
CnF^^^,4.1ECP^-CCF^^^^,H:IlcH|i2^XL^NC•^ 


S1«S2 
XSI-XS 
..  YSI*VS 

U*2  CCMlKUe 

240  CURTINUE 

c pfbwcnic  cqeff 

CUf < 1.2.  1ECF)*0. 

CCFI  1,4.  lECP>«iF. 

S1»^  . 

00  3 •'  I»L1 ,K 

x$«xe(i! 

ys«veii» 

«:!?|:{l!!^!ll:j  I 

QX>XS*FX«S2 
OY»Y5-FV*S2 
f?».;*FX*S2**24CX*S2 
T>», i*PX*S1 **24CX*S1 

5^I’0;{ii:JiSEUsV’ 

T1».5,*FY*S1  ••24CV*S1 

ICCP>«CCFI  I.  :,  lECPUC  T2-rn/XLENG 

XSImXS 

VS1«YS 

30-J  CCKTINOE 
RETUPR 


rn 


• I5C  eCXTpIlPSxI  11^''°*’®^*-^'"'*^^'’*  ><.’PI.MEft,C0F(1?,4,ia»  .GAT^I^IpJ, 
DIp5r5ICR*ElIs<3» 

INTEG 


MrEG£R*2  lECAZ 
bHAR  lECRNG 
CHAR  GATE.IPSK 

TAM(«-1./XP 

ihnvi 

A2>A«« 

B2«e«e 


V 

? 


AC»A*C 
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ec»B*c/ 
io>e«c 
co«c*c 
OEB>«”*(eC*C 


• l*C2  )-'•€*(  82*CC-AB*02)4C2*  I32«  IA04qC)-AB«e0) 
)-AC*l  XK*CC )4C2* IXK*|A04«C) ) 


XMSCATI  A24e24C2402) 

FNI-  XK*I  ec*cc-c  AC4?  C fC2  . . „„ 

Fftt>A2*|  XK*CC )-XK*IBiALD-AP*oi  )4C3*(-AB*XK  J 
El^stf* * FM/cIn  * ***‘  * -AC*«-AB*XK  »4XK*  ( B2*  ( Al? 


AL'4HC)-AB*BC) 


;LFSI2)>  . 
:tFSn)«FN3 


|/CEN 

N«Af  AAiELPf  l2)/ I ELPSm-ELPSI  i<m/2. 

TAAw-TAMm) 

SiKM«SlftlW) 


CQSM>CCSiW) 

^o2!!I“I!'E5!  I!!59§^!*i*^>-Ps«2»*siNw*cosw4ELPS(3i*stuh**2 

TANkt>-1./XM 

PETUPK 

E\C 


